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ABSTRACT 

We present the results of an imaging observation campaign conducted with the Subaru 
Telescope adaptive optics system (IRCS+A0188) on 28 gravitationally lensed quasars 
and candidates (23 doubles, 1 quad, 1 possible triple and 3 candidates) from the SDSS 
Quasar Lens Search. We develop a novel modelling technique that fits analytical and 
hybrid point spread functions (PSFs), while simultaneously measuring the relative 
astrometry, photometry, as well as the lens galaxy morphology. We account for sys- 
tematics by simulating the observed systems using separately observed PSF stars. 
The measured relative astrometry is comparable with that typically achieved with the 
Hubble Space Telescope, even after marginalizing over the PSF uncertainty. We model 
for the first time the quasar host galaxies in 5 systems, without a-priory knowledge 
of the PSF, and show that their luminosities follow the known correlation with the 
mass of the supermassive black hole. For each system, we obtain mass models far 
more accurate than those previously published from low-resolution data, and we show 
that in our sample of lensing galaxies the observed light profile is more elliptical than 
the mass, for ellipticity > 0.25. We also identify eight doubles for which the sources 
of external and internal shear are more reliably separated, and should therefore be 
prioritized in monitoring campaigns aimed at measuring time-delays in order to infer 
the Hubble constant. 
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1 INTRODUCTION 


* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan, 
f E-mail: cerusu@ucdavis.edu 
f E-mail: masamune.oguri@ipmu.jp 


The_ first _ gravitation ally_ lensed quasar 

(|Walsh, Carswell, fe Wevmannl Il979h has been discov¬ 
ered more than 30 years ago, turning gravitational lensing 
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from an obscure theoretical field into a mainstream ob¬ 
servational one. More than 100 strongly lensed quasars 
have been discovered to date, and it has convincingly 
been demonstrated that these objects provide insights 
into various topics in astrophysics and cosmology, as well 
as being a unique tool for studying the dark universe. 
Applications include the stu dy of the q uasar host galaxies 
at high redshift (e.g., iPeng et al.ll2006l ), dark matter sub¬ 
structures and lumin o us satellites (e. g., lChiba_et^_alJ [200^; 
McKean et al.l 120071 : iMinezaki et al.l 120091 ; Vegetti et al.l 


201^ . the_ structure and ev ol ution of massive g a laxies 


fe.g.. IKochanek et al.l l2000l; iRusin fe Kochanek 120051 : 
iKoopmans et alT ~ 20091 : lOguri. Rusu. fe Falcol 20141) . and 
microlensing applied in t he study of the st r ucture of 
quasa r accretion disk (e.g., lAgol fe K rolikl [1 9991: iDai et al.l 


2010fl. broad line regions(e.g., Richards et al. 2003 : 


Sluse et al. i mi) , as well as to measure the stellar mass 


Mediavilla et al.l 2009); 

Schechter et al. 2014). Follow- 

ing early work by e.g. 

Turner, Ostriker. & Gott 

1983); 

Fukugita, Futamase, & Kasai (19900; Kochanek 

1993), 

analyses of statistically 

well-defined samples of lensed 


quasars (i.e. samples in which lens candidates are selected 
by a homogeneous method whose completeness is known) 
can now constrain the cosmological constant/dark energy 
by comparing the number and distribution of image sep¬ 
aration of le nsed quasars with theoretical models (e.g., 
lOguri et alJ l2012il . Time delay measurements between 
quasar images constrain the Hubble constan t free of the 
calibration in the cosmic distance ladder (e.g., IOgurill20Q7l) . 
Finally, the distribution of lensed image separations, from 
galaxy to cluster mass scales, reflects the hierarchical 
structure formation and the effects of baryon cooling (e.g., 
IKochanek fe Whitell200ll l. 


Th e Sloan Digital Sky Survey Quas ar Lens Search 
(SQLS: lOguri et ahl 20061 : Ilnada et al.ll2012l ) is a systematic 
survey for lensed quasars, aiming to construct a large sample 
of gravitationally lensed quasars at optical wavelengths. It 
relies on the large homogeneous sample of spectroscopically- 
confirmed quasa rs from the Sloan Digital Sky Survey (SDSS; 
I York et al.ll2000l) . The techniques employed by the SQLS to 
identify lensed quasar candidates are described in the ref¬ 
erences above. We followed up all candidates with deeper 
imaging (usually with the University of Hawaii 2.2m Tele¬ 
scope; UH88) to detect the lensing galaxy. We then per¬ 
formed follow-up spectroscopy of the most promising can¬ 
didates, to confirm their lensing nature. SQLS is at present 
the prominent search for lensed quasars in the optical, com¬ 
prising of 62 lensed quasars t o date (December 2014), ~2/3 
of wh ich are new discoveries dKavo et alJ 2010: Ina da et al.l 
I2012lf1 . It has also produced the largest statis tically com- 
plete sample of lensed quasars (26 objects; Ilnada et all 
120121) . 

A disadvantage of SQLS, like other ground-based opti¬ 
cal strong lens surveys, is its poor detail in imaging lensed 
quasars. Even when performing follow-up observations with 
the UH88 telescope, the pixel scale ~ 0.22” is large, and the 


1 A publicly available list of all lensed 
quasars in the SQLS is available at 
http://www-utap.phys.s.u-tokyo.ac.jp/~ sdss/sqls/ 


seeing ~ 0.8” is similar to the image separation of a galaxy- 
scale strong lens (~ 1” — 2”). Therefore, high-resolution 
imaging of these quasar lenses is the key to turning each lens 
into a highly useful astrophysical and cosmological probe. 
This is necessary for obtaining accurate relative astrometry 
and point/surface photometry for the quasar images, lens¬ 
ing galaxy, and quasar host galaxy (in case the latter is 
detected), which are used to constrain both the light and 
the mass distribution in these systems. 

In the following, we highlight three of the applications 
enabled by the high-resolution images of a large sample of 
objects, such as the sample provided by our work. 

Estimating the Hubble constant from high resolution 
and time delay measurements: Although monitoring obser¬ 
vations, for which only relative fluxes are relevant, can be 
performed with small telescopes aided by image deconvolu¬ 
tion, to determine time del ays measurements between mul¬ 
tiple images, fe.g.. lEulaers et ahll2013l : iRathna Kumar et al.l 
120131) . high resolution single epoch observations are sti ll re¬ 
quired to construct an accurate lens model (e.g. JSuvu et al.l 
l2012ll . As time delays are currently being measured by the 
COSmologic al MOnitoring of GRA vItational Lenses (COS- 
MQGRAIL: lEigenbrod et al.l [20051 ) for many of the SQLS 
l enses, high r e soluti on data is in demand. For example, 
ICourbin et al] (1 19971 ) obtained high resolution images of a 
four-image lensed quasar with an early adaptive optics sys¬ 
tems, resulting in a relative lens galaxy position three times 
more precise than before, which allowed to measure the Hub¬ 
ble constant two times more precisely than in previous stud¬ 
ies. 

Quasar galaxy hosts and the correlation with Mbh.' 
The tight correlations found between the masses of su- 
permassive black holes Mbh and overall properties of 
the host galaxy bulges, such as vel ocity dispersion , lu¬ 


minosity and stellar mass 

e.g.. Kormendv & Richstone 

19951: 

Ferrarese & Merritt 

20001; Haring & Rix 2004; 

Gultekin et al. 201l|) suggest 

that the black hole growth 


is coupled to the galaxy evolution. These correlations 
have been established for local galaxies, based on spatially 
resolved kinematics. Since spatially resolved kinematics are 
very difficult to obtain at z > 1, the most straightforward 
correlation that can be studied at large redshifts is the one 
between Mbh and the bulge galaxy luminosity. However, 
in this case AGNs or quasars must be employed, as the 
nuclear activity allo ws the estimation of Mbh using the 
virial technique (e.g., lHolll999l : IVestergaardll20C)d ) . 

The difficulty then lies in accurately decomposing the 
faint quasar host, which is subject to the cosmological sur¬ 
face brightness dimming , from the brigh t nuclear source. As 
was demonstrated by IPeng et al.l fl2006l ). this task is facili¬ 
tated for gravitationally lensed quasars by two effects. First, 
strong gravitational lensing produces natural magnification 
typically by several factors, by increasing the total flux and 
the apparent size of the host galaxy, while conserving sur¬ 
face brightness. Second, lensing distorts the host into an arc 
and gives it a dramatically different morphology from that 
of the Point Spread Function (PSF), making it significantly 
less susceptible to systematic problems in the PSF model 
(an inaccurately known PSF). High-resolution observations 
are paramount in measuring the surface brightness of the 
lensed host galaxy. 

The ellipticity and orientation of mass and light in 
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the tensing galaxies: An ongoing question is to what ex¬ 
tent there are correlations between the azimuthal luminous 
profile of the lens galaxies and their overall (dark + lumi¬ 
nous) mass profile. Strong lensing provides a unique tool 
for testing these correlations. There exists a known corre¬ 
lation bet ween the mass and lig ht orientation, withi n ~ 10 
deg (e.g.. iFerreras. Saha, fe Buries! 1200^ : iTreu et al.ll2009lh 
It is not presently known however if the scatter in the re¬ 
lation is intrinsic, or due to measurement errors. On the 
other hand, there is less consensus on a correlation be¬ 
tween the mass and light ellip ticity. For example, su ch a 
tight correlation was found by iKoopma ns et all (120061) for 
the Sloan Lens ACS Survey (SLACS: iBolton et alJ 20051 ') 
sample of galaxy-galaxy lenses (i.e. foreground galaxies act- 
ing as strong lenses for background galaxies), although 
iFerreras, Saha, fc Burled (120081) did not find a correlation for 
a selected subsample of the SLACS lenses. A c orrela tion was 
also f ound for the SL2S galaxy-galaxy lenses llGavazzi et al.l 
l2012h . but has significantly more scatter. There are also re¬ 
cent results from different samp les (e.g., iDve et ak 1 12014k 
iKostrzewa-Rutkowska et al.l 1201^ 1 . The difference between 
the SLACS and SL2S lenses is that the former lie in an 
environment with significantly less average external shear 
( 7 ) 0.035, and also that they have smaller Einstein radius 

f?Ein relative to the characteristic scale (effective radius) of 
the galaxy A e fr, resulting in the more relevant role played 
by stellar mass in defining the potential within the critical 
curve (REin/Reft ~ 1.1, 0.5 for SL2S and SLACS, respec¬ 
tively) . For much l arger v alues o f REin/R e s, studies b ased on 
weak lensing (e.g. Ihloekstra. Yee. fc Gladdersi |2004) found 
that the dark matter haloes, which dominate the mass pro¬ 
file at these radii, are flatter than the luminous profile, in 
agreement with theoretical expectations. 

For the lenses of strongly lensed quasars, initial Hub¬ 
ble Space Telescope ( HST ) observations have failed to 
find a correlation between the m ass and light ellipticity 
llKeeton. Kochanek. fc Falcdl 1998 0. More recent results em¬ 
ploying a uniform analysis te chnique and decon volution of 
HST data are presented in ISluse et al.l ()2012h and show 
that there is a correlation for their sample of four-image 
lensed quasars (quads). The SLACS lenses and the lenses of 
strongly lensed quasars are known to probe different popu¬ 
lations of lenses, with lensed quasars typically lying in richer 
environments fe.g.. lOguri et al1l2005i 'l. In addition, for most 
lensed quasars f?Ein is on the order of a few R e g, therefore 
pro bing an int e rmed iate range. The new correlation found 
by ISluse et al.l d2012i ') . after eliminating the outliers, holds 
only for ellipticities < 0.25. For these small ellipticities, as 
the authors note, one possible bias resulting in the correla¬ 
tion is that the effect of a rounder dark matter halo on the 
total mass distribution may be harder to detect. 


The goal of this work is to enhance the value of the 
SQLS gravitationally lensed quasar sample through the 
use of high-resolution observations. Among the 62 lensed 
quasars in the SQLS sample, only ~ 22 have available high- 
resolution observations, mainly supplied by HST. However, 
recent progress of adaptive optics (AO), especially the use 
of laser guide stars (LGS), makes it possible to obtain high- 
resolution images using ground-based telescopes. Due to the 
smaller diffraction limit, near-infrared imaging with Subaru 


Telescope AO can potentially achieve three times the spatial 
resolution of HST, while also employing a finer pixel scale. 

AO has previously been used in the literature in 
the st u dy of gravitationa lly le nsed quasars. |Ledoux et al.l 
ill 9981) , lAuger et al.l ll2008lf and ISluse et all d2008l'l used AO 
to study ga la xy sc al e two-image len s ed qu as ars (doubles). 
iRusin et al.l (l2000l ~): ICourbin et al.l (l2002l ~): iMevlan et al.l 
1 20051) discovered additio nal lensing g alaxie s in known 
lensed quasars, using AO. iFaure et al.l (1201 ill observed a 
quad in order to look for substructure, which McKean et alJ 
(120071) was successful in dete cting for another quad. In ad- 


:r q uad. 1 

1. feoid) 


have 


dition, iMarshall et al.l (l2007l j and ISuvu et al.l 
used AO to study the arcs of strongly lensed galaxies, recon¬ 
struct the source, or study in high detail the lensing galaxy. 
However, these have been isolated studies focused on single 
systems. The observational campaign described in this pa¬ 
per is the largest dedicated study of a large sample of lensed 
quasars with AO. 

In Section [5] we present our observing strategy and de¬ 
scribe acquired data including data reduction. In Section [3] 
we present the morphological modelling technique we em¬ 
ployed, and in Section [|] the mass modelling technique. We 
continue with a comparison of our results with those found 
for the same systems in the literature (Section (5j, and a 
concluding discussion of the technique we used, in Section [5] 
Section[T]shows results obtained on the quasar host galaxies. 
Finally, Section [8] studies the lens environment and the rela¬ 
tion between mass and light for the systems in our sample, 
and Section [9] concludes this work. In the main appendix, 
we describe in detail the analysis of each individual system. 


Throughout this work, the concordance cosmology with 
H 0 = 100 h km" 1 s” 1 Mpc" 1 , h = 0.70, = 0.27 

and Q.a = 0.73 is assumed. For SDSS J1001+5027 and 
SDSS J1206+4 332, where time delay measur ements are 
available, the lPlanck Collaboration et all ll2013l ) result, h = 
0.673 ±0.012, is assumed. Magnitudes are given in the Vega 
system. All object coordinates assume J2000. All observed 
images and residual plots are in logarithmic scale. 


2 OBSERVATIONAL STRATEGY AND 

ACQUIRED DATA 

2.1 Observational strategy 

Sky coverage is a limiting factor when using AO, as all tar¬ 
gets must be located close to a tip-tilt (TT) star of suitable 
brightness in order to perform AO correction, even in the 
case that LGS is used. Subaru Telescope however presently 
has the most relaxed constraints on TT star, among com¬ 
parable systems on 8-10 m telescopes (Section 12. 2. 111 . This 
makes it an ideal tool to use in the observational campaign 
of the SQLS quasars. Indeed, 54 of the 62 SQLS quasars are 
accessible to the Subaru Telescope AO system. 

The AO correct ion functions better at longer wave¬ 
lengths (e.g. iDavie IHH), making A'-band the natural 
choice among the widely used JHK near-infrared bands. An 
additional reason to use A'-band in the study of quasars and 
their host galaxies is that for redshifts < 2.5, A'-band falls in 
a wavelength region where the quasar host galaxy is bright¬ 
est compared to the nuclear source (Figure[T|). Moreover, for 
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Figure 1. Typical galaxy (E-type) and radio quiet quasar spec¬ 
trum, with the location of the K band at redshifts 0 and 2.5 
marke d. The composite quasar spectrum is from IShang et al.l 
^Olll h For t he galaxy template, the m e an sp ectrum of local 
galaxies from IColeman, Wu, Sz W eedman] (|l980| ) is extrapolat ed 
using the evolutionary models of lBruzual A. &; Chariot! (|l993l ). 


the purpose of modelling gravitationally lensed quasars, ex¬ 
tinction, intrinsic variability and microlensing are all weaker 
at progressively longer wavelen gths, contaminating the true 
imag e flux ratio less (e.g. lYonehara , Hi rashita, fc Richter 
12008 ). However, exposure times in the iv-band are limited by 
the strong sky background emission. On the other hand, K'- 
band provides a good compromise between background level, 
sensitivity, and AO performance, and we therefore used this 
band throughout the observation campaign. 

As one of the goals of the campaign is to measure ac¬ 
curate lens galaxy morphologies and potentially resolve the 
quasar host galaxies from the point-like nuclear source, a 
good knowledge of the PSF is essential. For the majority 
of the systems, we used the SDSS CasJobt@ search query 
to identify bright PSF and guide star pairs. These were se¬ 
lected such that the guide stars of the PSF and the target 
have similar brightness, and the separations between the tar¬ 
get and its guide star, as well as the PSF and its respective 
guide star, are similar. The PSFs are close on the sky to 
the target (within ~ 30'), and they were observed immedi¬ 
ately before or after the target, in the same natural guide 
star (NGS) or LGS mode. This was done to ensure as much 
as possible that the atmospherical turbulence characteris¬ 
tics and therefore the AO correction is relatively constant 
between the PSF and the target. However, as we will discuss 
below, these separately observed PSF stars turn out not to 
be suitable for the accurate characterisation of the target 
image PSFs. SDSS J1515+1511 and SDSS J1002+4449 are 
the only observed systems that have a bright star in the de¬ 
tector field of view (FOV), which can potentially be used 
as PSF. Here we centred the FOV, and therefore the LGS, 
between the star and the target. This was done, again, so 
that the AO correction at the location of the target and the 
PSF star is similar. 


2.2 Acquired data and data reduction 

2.2.1 IRCS+A0188 

The AO imaging observations were perfor med with The In ¬ 
frare d Camera and Spectrograph (IRCS; iKobavashi et al.1 
l2OO0ll at the Subaru telescope, along with th e Laser Guide 
Star Adap t ive Optics system ( LGS+AOI 88 ; lHavano et al.l 
I2008L l201(ll: iMinowa et al.ll201(tl . IRCS uses two 1024 x 1024 
ALLADIN III arrays for imaging and spectroscopy, and 
in imaging mode provides two plate scales, of 52 and 20 
mas per pixel. After the distortion correction described in 
Section 13.71 the pixel scales are 0.05284" ± 0.000017" and 
0.02053" ± 0.000005", respectively. The available FOVs are 
54" and 21", respectively. 

AO 188 uses a curvature sensor with 188 control ele¬ 
ments, operating at 2000 Hz, and a 188 element bimorph 
deformable mirror (BIM188). It operates at the Nasmyth fo¬ 
cus of the Subaru Telescope. Both NGS and LGS modes are 
available. The recommended guide star constraints in NGS 
mode are brightness larger than R = 16.5 mag, within 30" 
from the target. The LGS mode use s an artificial sodium 
laser guide star llHavano et ali l2006l l for high-order wave- 
front sensing. The LGS is typically less bright than ~ 10.5 
mag in R-band. The recommended tip-tilt star constraints 
in LGS mode are brightness larger than R = 18 mag, within 
60" (acceptable up to 90") from the target. More informa¬ 
tion about AO 188 is provided on the instrument web pagc0. 


2.2.2 Acquired data 

The observations were performed between 2011-2014. In to¬ 
tal, ~ 7 nights were assigned for the present campaign, 
about half of which were lost due to telescope/instrument 
trouble or cloud coverage. Observation priority was given 
based on the scientific interest of each target and the de¬ 
gree to which the targets are suitable for AO observations 
in terms of TT star brightness and separation, and dis¬ 
tance to zenith/airmass. The 28 objects that were success¬ 
fully observed, the observation modes, exposure times, fil¬ 
ters, pixel scales, typical airmass, photometric stars and ob¬ 
servation dates are given in Table |T] Two of the objects, 
SDSS J0820+0812 and SDSS J1206+4332, were observed as 
back-up, without AO. 

All observations were performed with 5 or 9-point 
dithering, in order to remove bad pixels and cosmic rays, 
and to allow for flat frame and sky frame creation from the 
data. Although dome fiats were obtained for some observa¬ 
tions, these have proven to be of inferior quality. Although 
the targets were chosen to be accessible to AO, many of these 
border the recommended limit in terms of the TT bright¬ 
ness and separation. As such, low Strehl ratios < 10% (Table 
IA2I) were typically obtained. The 52 mas pixel scale mode 
was typically used, in order to increase the signal-to-noise 
(S/N) and not to avoid significantly oversampled PSFs. The 
list of observed TT and PSF stars is given in Table IaD 
For 12 objects, PSF star observations were skipped because 
of rapidly changing atmospheric conditions. Typical seeing 


2 http://skyserver.sdss3.org/casjobs/ 


3 http://www.naoj.org/Observing/Instruments/AO/system.html 
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during the observations can be inferred from Table IA2I and 
was generally ~ 0.5” — 1". 


2.2.3 Data reduction and calibration 

Data reduction was performed with IRAF0, using the IRCS 
IMGREE0 package designed to reduce data obtained with 
IRCS. The reduction consisted of the following steps: 

(i) Each frame was checked for dark patterns by divid¬ 
ing it to the next observed frame. Also, the number of 
counts at the location of the target was checked for lin¬ 
earity. All targets are well inside linearity limit, except for 
the core of the star inside the FOV for SDSS J0743+2457 
(within 5%), and cores of the bright quasar images in 
SDSS J0904+1512 (within 3%), SDSS J1322+1052 (within 
2%) and SDSS J1353+1138 (within 4%). 

(ii) As persistent fringe artifacts were found in the dome 
flats, sky flats were created by masking bad pixels and bright 
objects, then median-combining the dithered frames which 
were divided by the mean count value. Each raw frame was 
corrected for different pixel sensitivities by division to the 
sky flat. 

(iii) The geometric distortion correction map described in 
Section nm was applied. 

(iv) Sky-background frames were generated by median¬ 
combining the flat-fielded frames. As the near-infrared sky 
is highly time-variable, the sky-background frames were gen¬ 
erated for each dither sequence individually. Pixel values for 
the masked pixels were interpolated from the combined ones. 
Before subtracting the sky frame, each frame was normalized 
so that their median values match each other. 

(v) Position offsets between the images were obtained via 
cross-correlations between regions that contain bright ob¬ 
jects, typically the targets. Cosmic rays falling in the re¬ 
gions used to perform the cross-correlations were manually 
masked. The full-width at half-maximum (FWHM) of star- 
like objects in each frame was measured by fitting a Moffat 
function to the radial profile with the IRAF IMEXAM task, 
and the frames with significantly different values were dis¬ 
carded. The frames were average-combined using rejection 
algorithms (average sigma clipping or clipping based on the 
CCD parameters). Pixel binning was performed for a few 
frames with significantly oversampled PSF, as mentioned in 
Appendix m 

(vi) Photometric zero-point calibration was performed, 
typically using standa rd stars. All o bjects were corre c ted fo r 
Galactic extinction dSchlegel, Finkbeiner, fc Davis! 11998[) . 
and atmospheric extinction relative to the standard staiQ 

It must be noted that the standard star catalogue magni¬ 
tudes at A'—band were used to calibrate the photometry at 


4 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 

5 The package is available at 

http://uv#.naoj.org/Observing/DataReduction/index.html 

B According to the Mauna Kea summit extinction values 
(http://###.jach.hawaii.edu/UKIRT/astronomy/utils/exts.html ) 
the atmospheric extinction is of ~ 0.1 mag/airmass for A-band. 


A''—band. However, the expected differences based on inter¬ 
polation are small (0.01 — 0.03 mag). 


3 MORPHOLOGICAL MODELLING 

TECHNIQUE 

After reducing the data, it became clear that none of the 
separately-observed PSF stars represent suitable PSFs to 
model the corresponding systems. This was also concluded 
by other stud ies employing separate ly observed PSF cal¬ 
ibrators (e.g., iKuhlbrodt et alll2005i ). The reason is PSF 
variability, both in the science targets and the separately- 
observed stars. This is caused by the rapid change of the 
atmospheric turbulence characteristics, which induce dif¬ 
ferent responses from the AO system. Figure [2] shows the 
time variability of the core component FWHM in the PSF 
for both the target and the separately-observed star in 
SDSS J0946+1835, as well as the outstanding residuals ob¬ 
tained if the star is used as a PSF to model the system. 

It is imperative that we construct a well-characterised 
PSF for each system, in order to reach the science goal, i.e. 
obtain an accurate astrometrical, photometrical and mor¬ 
phological characterisation of the multiple components of 
each observed system. In order to guide the reader, here 
we provide a short summary of the following subsections, 
in which we illustrate in detail our PSF reconstruction and 
modelling technique. In section [3711 we give general consider¬ 
ations on the AO PSF. In Section [3. 2 1 we briefly mention the 
AO PSF modelling techniques encountered in the literature. 
In Section we introduce our technique of modelling an¬ 
alytically all components of a gravitationally lensed quasar 
simultaneously, as their light profiles overlap, and derive the 
most suitable analytical profile. We describe the implemen¬ 
tation of this technique, where we take special care to prop¬ 
erly explore the parameter space. In Section 13.41 we show 
that for selected systems we can go beyond an analytical 
PSF, and we derive a hybrid PSF that specifically accounts 
for the non-analytical components. In the next Sections (1331 
and 13.61 respectively), we describe how we need to modify 
out technique to model systems where we detect a quasar 
host galaxy, as its profile fitting requires a mass model. We 
then insure that the error estimates for each derived param¬ 
eter accounts for systematics introduced by our modelling 
technique. Finally, in Section f3.7l we address the issue of how 
the instrument held distortion impacts our derived astrom¬ 
etry. 

3.1 The adaptive optics PSF 

For a seeing-limited image, spatial resolution is usually 
characterised by the FWHM of a point-like stellar profile. 
This characterisation becomes unreliable near the diffrac¬ 
tion limit, as the FWHM measurement is complicated due 
to the diffraction rings. In the diffraction limit, the shape 
of the PSF is described by a two-dimensional Airy function 
(ignoring the telescope spider), with the first diffraction ring 
at an angular distance of 1.22X/D from the centre; here A 
denotes wavelength, and D is the diameter of the telescope 

' http://www2.keck.hawaii.edu/inst/nirc/UKIRTstds.html 
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Table 1. Summary of observations 
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Object 

Guide star 

Exposure 

[B] 

Filter 

Pixel scale 
[mas] 

Airmass 

Photometric 

calibration 

Observation date 
(UTC) 

SDSS J0743+2457 

NGS 

20 x 60 

K' 

52 

1.13 

yes“ 

2011 November 29 

SDSS J0819+5356 

NGS 

14 x 60 

K' 

52 

1.21 

yes c 

2012 February 19 

SDSS J0820+0812 

w/o AO 

14 x 60 

K' 

52 

1.38-1.75 

yes 6 

2012 February 19 

SDSS J0832+0404 

LGS 

16 x 60 

K' 

52 

1.07 

FS127 

2014 January 14 

SDSS J0904+1512 

LGS 

17 x 60 

K' 

52 

1.05 

FS126 

2013 April 28 

SDSS J0926+3059 / 

LGS 

11 x 60 

K' 

52 

1.06 

no 

2014 January 15 

SDSS J0946+1835 

NGS 

6 x 60 

K' 

52 

1.08 

yes“ 

2011 May 17 

SDSS J0946+1835 

NGS 

33 x 240 

K' 

20 

1.03-1.27 

P550C 

2012 February 2 

SDSS J1001+5027 

LGS 

27 x 40 

K' 

52 

1.25-1.38 

FS127 

2012 February 20 

SDSS J1002+4449 

LGS 

16 x 50 

K' 

52 

1.10 

yes“ 

2011 November 29 

SDSS J1054+2733 

LGS 

16 x 60 

K' 

52 

1.02 

FS127 

2014 January 14 

SDSS J1055+4628 

LGS 

8 x 60 

K' 

52 

1.12 

FS126 

2013 April 27 

SDSS J1128+2402 

LGS 

15 x 60 

K' 

52 

1.02 

FS127 

2014 January 14 

SDSS J1131+1915 

LGS 

7 x 60 

K' 

52 

1.00 

FS127 

2012 February 20 

SDSS J1206+4332 

w/o AO 

14 x 60 

K' 

52 

1.09 

FS126 

2013 April 28 

SDSS J1216+3529 

LGS 

9 x 60 

K' 

52 

1.05 

FS126 

2013 April 28 

SDSS J1254+2235 

LGS 

30 x 60 

K' 

52 

1.04-1.12 

FS23 

2012 February 21 

SDSS J1313+5151 

LGS 

48 x 25 

K' 

52 

1.19 

FS23 

2012 February 21 

SDSS J1320+1644 

LGS 

31 x 60 

K' 

52 

1.00 

no 

2013 April 28 

SDSS J1322+1052 

LGS 

24 x 60 

K' 

52 

1.26-1.40 

FS126 

2013 April 27 

SDSS J1330+1810 

LGS 

12 x 60 

K' 

52 

1.09-1.13 

FS127 

2012 February 20 

SDSS J1330+1810 

LGS 

13 x 240 

K' 

20 

1.05 

FS126 

2013 February 27 

SDSS J1334+3315 

LGS 

10 x 60 

J 

52 

1.07 

P272-D d 

2011 February 18 

SDSS J1334+3315 

LGS 

5 x 60 

H 

52 

1.09 

P272-D d 

2011 February 18 

SDSS J1334+3315 

LGS 

10 x 60 

K' 

52 

1.05 

P272-D d 

2011 February 18 

SDSS J1353+1138 

LGS 

23 x 60 

K' 

52 

1.31-1.47 

FS126 

2013 April 27 

SDSS J1400+3134 

LGS 

15 x 60 

K' 

52 

1.48-1.67 

FS126 

2013 April 27 

SDSS J1405+0959 

LGS 

19 x 60 

J 

52 

1.02 

FS126 

2013 April 27 

SDSS J1405+0959 

LGS 

26 x 60 

H 

52 

1.02-1.07 

FS126 

2013 April 27 

SDSS J1405+0959 

LGS 

14 x 60 

K' 

52 

1.10-1.14 

FS23 

2012 February 21 

SDSS J1406+6126 

LGS 

21 x 60 

K' 

52 

1.46 

no 

2014 January 15 

SDSS J1455+1447 

LGS 

23 x 60 

K' 

52 

1.08-1.18 

no 

2014 January 15 

SDSS J1515+1511 

LGS 

22 x 40 

K' 

52 

1.01 

FS23 

2012 February 21 

SDSS J1620+1203 

LGS 

2 x 60 

K' 

52 

1.13 

no 

2011 August 8 

SDSS J1620+1203 

LGS 

16 x 60 

K' 

52 

1.16-1.22 

2MASS e 

2013 April 27 


Exposure times show the number of frames that were combined in the final science frame. a Based on 1 star in the 
FOV; b based on 1 star and 1 galaxy in the FOV; c based on 3 stars in the FOV; d dark sky pattern may affect model 
photometry; e compared the total magnitude with the measurement in 2MASS; f refers to SDSS J092634.56+305945.9; 
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Figure 2. Top: PSF FWHM for each individual frame in 
the SDSS J0946+1835, as well as for the separate PSF star. 
Frames are numbered chronologically. Bottom left: original 
SDSS J0946+1835 observed frame; bottom right : residuals after 
modelling the system with the separate PSF star, coadded from 
its first four frames. 


use of the typical structure of a lensed quasar: two or four 
point-like images with a lensing galaxy in between. While 
the quasar point-source is safely unresolved, since it is typ¬ 
ically located at high redshift, the quasar host galaxy may 
or may not be resolved. In all cases, only doubles where the 
host galaxy is unresolved were modelled in the literature, 
and in the case of quads, the host was actively removed 
during modelling, without estimating any of its physical pa¬ 
rameters;_ _ 

ISluse et all (|2008l ) modelled the double SDSS 
J0806+2006 with the PSF estimated from the more 
isolated, brighter image, and employing image deconvolu¬ 
tion. They do not mention h owever if any system atics are 
introduced by this technique. I Davies et ali ( 2004 1 however 
argued for the use of convolution against deconvolution, 
since the PSF is poorly defined, the S/N is limited, and 
smooth extended emission (in the case of hosts) makes 
the results uncertain. When there is a physical model to 
test, they advise that the best approach is to convolve 
that model with the best estimate of the P SF an d then 
compare the result to the observations. iFaure et alj (1201 lh 
took advantage of the cross configuration of the quad SDSS 
J0924+0219, where the host galaxy has a different mor¬ 
phology at the location of each image, and estimated the 
PSF by combining the super posed and normalised i mages, 
in an iterative procedure. ILagattuta. Auger, fc Fassnachtl 
rt20ldf analysed the quad B0128+437 by iteratively fitting 
a PSF modelled as three Gaussian components, together 


with the lensin g galax y modelled with a Sersic profile. 
iKoptelova et al.l ( 2014) 4 modelled the PSF for a lensed 
double by using a hybrid composed of the observed bright¬ 
ness distribution of the bright image at the centre, and 
analytical wings. 

3.3 Modelling with an analytical PSF 

There is enough information (e.g., [Marshall et al.l l2007h 
in the multiple images of a lensed system to reconstruct 
the PSF. In this paper, we combi ned t he approach of 
ILagattuta. Auger, fc Fassnachtl d2010l 4 and IKoptelova et al.1 
( 2014 ) into a new technique, by adopting an analytical PSF 
model and simultaneously optimising for the global PSF pa¬ 
rameters as well as the individual objects comprising the 
system: multiple point sources, and the lensing galaxy mod¬ 
elled as a Sersic profile convolved with the PSF. Where pos¬ 
sible, we take an additional step where we refine this fit by 
using a hybrid PSF (Section 13.411 . The o ptimisation is per - 
formed using a downhill simplex method dPress et al.lll993) 
to find the minimum goodness-of-fit y 2 in parameter space, 
and is implemented in a new software, Hostlens, which was 
developed by one of us (M. Oguri). Host lens co mbines the 
functionality of Galfit with that of Glafic (Ogur Like 

Galfit, it allows the simultaneous fitting of both point and 
extended sources, but in addition can simultaneously op¬ 
timise the parameters of an analytical PSF. Like Glafic, it 
can fit an observed lensing configuration with a mass model, 
while simultaneously fitting analytically the morphology of 
an extended lensed source, such as a quasar host galaxy 
(Section 13.51 . We note that most of the functionalities of 
Hostlens are currently implemented in Glafic. 

We considered several analytical PSF models to use 
with Hostlens. According to Section 13.11 the AO PSF 
consists of two components, one describing the compact 
PSF core, and the o ther one the extended seeing wings. 
lAkivama et al.l (|2008l 4 used a single Moffat profile, whose 
wings are more extended than that o f a G a ussian , to ap¬ 
proximate the full AO PSF. I Davies et al.l d20ld) stated 
that the PSF is generally matched quite well by the sum 
of a narrow Gaussian for the core, and either a Mof - 
fat or another Gaussian for the halo. iFalomo et al.l (120081 ) 
used a Gaussi an for the core and an exponential fu nction 
for the wings. ILagattuta, Auger, fc Fassnachtl J20ld ). men¬ 
tioned above, used three Gaussians: one component rep¬ 
resents the diffraction-limited core, another represents the 
seeing-limited diffuse PSF, and the third encodes structure 
in the PSF. In order to find the most appropriate model, we 
fitted the majority of bright stars observed during the obser¬ 
vation runs (i.e. the PSF stars) with analytical profiles. Two 
of these fits are shown in Figures 0 and [I] The first figure 
shows a two-dimensional fit and the subsequent residuals af¬ 
ter subtracting elliptical but concentric 1 Moffat, 1 Gaussian 
+ 1 Moffat, 2 Moffat, and 3 Gaussian profiles, respectively. 
The second shows the fitted radial plot corresponding to the 
same models. As can be seen, both in terms of the formal 
X 2 /d.o.f., as well as in terms of the appearance of residuals, 
and radial profile matches, the 2 Moffat profile provides the 
best fit. This is found in the overwhelming majority of fitted 
stars, although in some cases, such as the second example 
shown, the 1 Gaussian + 1 Moffat profile provides a simi¬ 
larly good fit. The residuals in the other models all show a 
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concentric halo, and do not match the overall radial profiles. 
On the other hand, none of the models (including 2 Moffat) 
is able to model the non-analytical components close to the 
centre. This entirely empirical analysis resulted in the de¬ 
cision to use the 2 Moffat profile as the analytical PSF for 
Hostlens. The better fit of a Moffat core, which is more ex¬ 
tended than a Gaussian, may be due to the very small Strehl 
ratios < 10% typically obtained in the observations. Indeed, 
the typical core FWHM is ~ 0.2", about three times that of 

the diffraction-limited c ore. _ 

The Moffat profile (lMoffatlll969f i has surface brightness 
of the form 


I(r) = [l + (2 1/p - l) (2r/FWHM) 2 ] P , r = J x 2 + (^ 

( 1 )' 

where r is the elliptical radial distance, x ( y ) is the distance 
to the centre along the semimajor (semiminor) axis and q 
is the ratio of the semiminor to the semimajor axis; ft is 
a parameter specifying the shape of the profile, related to 
the core-to-halo flux ratio. The Moffat profile results into a 
Gaussian profile as the f) = oo case. In Hostlens, the PSF is 
modelled with 9 analytical parameters: two for each of the 
four components of an individual Moffat profile: FWHM, el- 
lipticity, position angle and /3, as well as one parameter spec¬ 
ifying the relative distribution of flux, fluxl/(fluxl + flux2). 
The two Moffat profiles are concentric. Each point source 
is modelled with three parameters: the positions along X 
and Y, and the total flux. In addition, every Sersic profile 
has ellipticity, orientation, effective radius and Sersic index 
parameters. Finally, there is a parameter for the sky value. 
Therefore, a typical two-image lensed system requires 24 
parameters to model analytically. We modelled each sys¬ 
tem with Hostlens, at least at an initial stage. The derived 
analytical parameters of all PSFs are shown in Table IA2I 
We selected suitably large cut regions (typically ~ 7" x 7") 
around the targets, and masked all objects that were not 
simultaneously fitted by Hostlens, in order to allow reliable 
sky value estimates and model fitting. 

Due to the large number of parameters and the possi¬ 
bility that the parameters might get trapped in a local \ 2 
minimum, we chose 500 sets of different initial parameter 
values from a range of plausible values, e.g. from a flat dis¬ 
tribution in the case of position angles, and a gaussian dis¬ 
tribution around the IRAF IMEXAM-derived value for the 
core FWHM component. We chose the model with the low¬ 
est value of x 2 (we also checked the second-best model to be 
similar), and in addition we performed a parameter search 
around it using Markov-chain Monte Carlo (MCMC), for 
the purpose of both refining the parameter values, and esti¬ 
mating confidence intervals (error bars). We ran ten MCMC 
chains, and removed a fraction of ~ 10% — 30% of the points 
from the head of each chain (the “burn-in stage” points; 
i.e sensitive to the start values). To test the convergence 
of these chains, we compared the variance of the distribu¬ 
tion of points in each chain to the va r iance of the combined 
distribution, following iGelman et al.1 (Il995l i. After applying 
a smoothing technique, the combined chains provided 68% 
confidence intervals on the analytical parameters (lcr in the 
Gaussian approximation). 

Throughout this work, unless otherwise stated, for all 


the systems where we did not detect a quasar host galaxy, 
and which we fitted with a Hostlens-produced analytical 
PSF, we also checked the modelling results for consistency, 
with Galfit. In order to produce with Galfit error bars 
that account for the uncertainties of the analytical PSF, 
we produced 500 realisations of the Hostlens-produced PSF, 
by drawing at random from the probability distribution of 
the analytical parameters of the PSF, obtained during the 
MCMC runs above; we then ran Galfit with each of the 
500 PSFs, and took the standard deviations of the result¬ 
ing parameters as an error estimate. We did not use Galfit 
for the systems with host galaxy detections, as it lacks the 
functionality to fit this component (see Section [3.51) . 

We note that Hostlens and Glafic measure the effective 
radius of a Sersic profile using a circularised radius, while 
Galfit measures the radius along the semi-major axis. The 
effective radii measured in this paper are reported following 
the Galfit definition. For the lensing galaxies, we considered 
both an unconstrained Sersic index n and n = 4 (Tabled , in 
order to ease comparison with other samples of lenses where 
th e light prof il e is modelled as a de Vaucouleurs profile (see 
IQguri. Rusu. fe Falcoll2014l . and references therein). 

3.4 Hybrid PSF 

The observed systems can be classified into two categories, 
according to the degree in which the bright image A is 
separated from the other objects, and the visual appear¬ 
ance of the modelling residuals. For systems with small 
separation where the non-analytical residuals at the cen¬ 
tre clearly overlap with a large fraction of the galaxy light 
(e.g. SDSS J1254+2235 and SDSS ,11334+3315), as well as 
for objects where, due to the low S/N, no conspicuous core 
residuals are seen, data modelling stops as described above. 
For the other systems, it is possible to model away the faint 
quasar image and the lensing galaxy using the analytical 
Hostlens profiles, and obtain a clean image A to be used as 
the “actual PSF” (panels a-g in Figure(5j. The advantage in 
this approach is that the estimated PSF is the “true” PSF 
(accounting for all non-analytical features), albeit a noisy 
one. Next, we tested whether this PSF produces better sta¬ 
tistical results. A simple comparison with the analytical fit 
in terms of y 2 cannot be used, as image A is fitted perfectly 
by design. Therefore, we compared the standard deviation 
of the residual pixel values in a region containing the rest 
of the objects (typically B and G). We found that in al¬ 
most all cases, the analytical PSF produced better results. 
This is due to large noise in the wing of the “actual PSF”. 
Th erefore, we used anothe r approach, previously employed 
by iKoptelova et al.l (120141 ') . This consists in using the ob¬ 
served brightness distribution of the bright image only in 
its central regions, and replacing the wings after a certain 
cut radius with analytical wings from the 2 Moffat fit (pan¬ 
els h-i in Figure 0. This “hybrid PSF” approach succeeded 
in producing fits at least as good as the analytical PSF in 
terms of the standard deviation of the residuals. 

We considered two sources of systematics for the hybrid 
PSF approach. First, because the contributions from B and 
G are initially modelled using the analytical PSF, they may 
be improperly subtracted at the location of A, and there¬ 
fore affect the hybrid PSF estimation. This can be checked 
by building the hybrid PSF in an iterative approach, where 
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Figure 3. Residuals after analytically modelling with Galfit the PSF stars for SDSS J0946+1835 (January 2012) (top) and 
SDSS J1515+1511 (down). From left to right: original star, residuals after subtracting 1 Moffat, Gaussian + Moffat, 2 Moffat, and 
3 Gaussians models, respectively. The d.o.f. and y 2 /d.o.f. corresponding to each fit are shown. 
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Figure 4. Radial profiles and residuals after analytically modelling the SDSS J0946+1835 (January 2012) (left) and SDSS J1515+1511 
(right) PSF stars. The radial profiles were generated with the IRAF task ELLIPSE, and the residuals with the task PRADPROF. A 
colour version of this figure is available in the online version of the paper. 


in each subsequent step the hybrid PSF is used to model 
away B and G, resulting in a better hybrid PSF estimate. 
We did this for SDSS J1405+0959 and obtained that subse¬ 
quent PSFs are virtually identical. This is expected, as we 
applied the hybrid PSF technique only to systems with iso¬ 
lated bright images. Secondly, the final models may depend 
on the size of the hybrid cut region. We chose the original 
cut radius as the radius at which all non-analytical compo¬ 
nents disappear, of ~ 10 pixels. For all hybrid PSF systems 
(excluding the ones with a quasar host galaxy detection, as 
described below, or otherwise specified), we included into 
the error budget the scatter between the results obtained 
from this cut size, and another one 50% — 100% larger. 

We used the analytical and/or hybrid PSF approaches 
described above for all systems where a quasar host galaxy 


was not detected. Table [2] shows whether the analytical or 
the hybrid PSF was used to model each system. 


3.5 Fitting the quasar host galaxy 

A novel undertaking in this paper is the fitting of quasar 
host galaxies, for several systems, without an a-priori known 
PSF. This is only possible for gravitationally lensed quasars, 
because the host galaxies are tangentially stretched around 
the critical curves by the lensing effect, resulting in arcs 
that are distinct from the PSF shape, and can be modelled 
analytically (Figure [Gj. 

In order to model the spatially extended flux distribu¬ 
tion of the lensed images, it is necessary to model both the 
lens mass distribution and the extended source distribution. 
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Figure 5. Example of a hybrid PSF creation, (a) Original frame. ( b ) The section of the frame containing most of the light with the 
exception of the bright quasar image is cut. (c) The cut section is replaced with nearby blank sky. ( d ) The best-fit analytical models of 
all objects except the bright quasar image, which were created with an analytical PSF, are removed from the rest of the frame, (e) The 
analytical wings of the bright quasar image are added to the section where the sky was replaced. (/) The resulting frame is now clear 
from all objects except the bright quasar image, (g) The resulting PSF frame is centred, and the sky pedestal is removed. The resulting 
image is fitted analytically. ( h ) Only a circular region around the PSF centre is kept, large enough to contain all visible non-analytical 
residuals, (i) The portion of the frame that was removed is replaced with the best-fit analytical wings. 


This is done using the approach described in iPeng et al.l 
(l2006fi . and makes use of the fact that gravitational lensing 
conserves surface brightness. 

Hostlens was designed with the functionality to fit host 
galaxies modelled with a Sersic profile. When fitting a 
host galaxy with Hostlens, the parameters typically used 
to model the individual images are replaced by a single 
source profile. This profile consists of the parameters of 
a Sersic model, its flux as well as the flux of the central 
point source, and the lensing parameters: the Einstein ra¬ 
dius, shear/ellipticity and their orientation angles, where a 
singular isothermal mass model is assumed (Section [4j. The 
luminous lens galaxy is fitted simultaneously, as before, with 
a Sersic profile. 

There are five cases where the quasar host galaxy was 
clearly detected and fitted. We define a clear detection as 
a system for which there are clear visible arcs similar in 
orientation to the critical curves of the lensing model, and 
which are convincingly removed by fitting a host galaxy. 

In the case of two systems with host detections, 
SDSS J0904+1512 and SDSS J1322+1052, there were out¬ 
standing residuals when using an analytical PSF, and there¬ 
fore we employed a hybrid PSF. This was possible through 
an iterative process, where an analytical PSF is initially used 


to obtain a first rough approximation of the host profile. 
The host galaxy is then subtracted from the original ob¬ 
served frame, and subsequently a hybrid PSF is created as 
described in Section[3]4] This hybrid PSF is used in the next 
step of the iteration to refit the original system, and refine 
the host model. Subsequent steps lead to new hybrid PSFs. 
The iteration stops once the y 2 stops decreasing. 

To summarise, the general approach to the morpholog¬ 
ical modelling of the imaging data in this paper, whether 
the quasar host galaxy is detected or not, is described in the 
flowchart of Figure [7] 

3.6 Simulations for testing the reliability of the 
technique 

To test the reliability of the morphological fitting with 
Hostlens using the analytical or hybrid PSF methods, we 
performed a series of simulations for each system. In these 
simulations, we used a separately observed PSF star, or a 
star in the FOV that is the most similar in terms of analyt¬ 
ical parameters to the PSF of the respective system (pri¬ 
oritising for the core FWHM parameter). Table 0 shows 
which star was used as PSF template for each system. With 
Hostlens and/or Galfit, we simulated each system using 
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Figure 6. A simulated typical two-image gravitationally lensed 
quasar with a detected quasar host galaxy (color-inverted image). 


this corresponding PSF template, and the best-fit astrom¬ 
etry/morphology/photometry parameters we derived when 
modelling that particular observed system. Such a simulated 
system is shown in Figure [6] From the final science frame of 
each observed (i.e. real) system, we identified 100 blank sky 
cuts, and added them as noise realisations to the simulated 
system. We then remodelled each of the 100 simulated noisy 
systems in the same way as we did the observed system. 
The scatter in the resulting 100 values for each fitted pa¬ 
rameter, with respect to the known input value, represents 
an error estimate. The only exception is for the case when 
we simulated the system using the best-fit parameters, but 
remodelled it using Sersic index n = 4. In that case, we 
considered the scatter around the mean fitted value (Table 
0 . 

We used two time-saving approximations in remodelling 
the simulated systems. First, we did not perform a search 
starting from random values in the parameter space for each 
individual simulated system, as we did for the original sys¬ 
tem, but instead we started the fit from the known (correct) 
values. Second, in cases where we used a hybrid PSF, we did 
not create it in each of the 100 simulations from that very 
same image, but created it from only one of the simulations, 
and then used it as the hybrid PSF to model all simulated 
frames. Also, for the few cases when a hybrid PSF is used to 
fit a detected host galaxy, we only performed a small number 
of simulations, as these are time consuming. 

These simulations can account for the systematics intro¬ 
duced by the analytical modelling of a non-analytical PSF, 
or the use of hybrid PSFs. The orientation of the detec¬ 
tor was typically the same for the science objects and the 
PSF stars, therefore diffraction patterns due to the telescope 
spider are similar. The noise c uts from the actu al sc ience 
frames, previously employed bv lGabor et all fe.g. l200it ). ac¬ 
count for how the model parameters would depend on the 
existence of background objects, too faint to be visible, and 
are more realistic than artificially-generated noise. 

As final errors on astrometry/morphology/photometry, 


we used the largest error bars among those resulting from 
these simulations and all other methods described in the 
previous subsections. This is because the methods are not 
necessarily independent, and one typically dominates the 
error budget. 

3.7 Accounting for the IRCS+A0188 field 
distortion 

The most important observational constraint for gravita¬ 
tional lens mass models of lensed quasars is provided by 
relative astrometry. In order to obtain accurate astrometry, 
it is necessary to account for the IRCS+A0188 field distor¬ 
tion, which results mainly from the pair of off-axis parabola 
mirrors. 

3. 7 .1 Creation of a distortion map 

The existing estimate of the geometrical distortion in the 
IRCS+A0188 was produced by matching stars in the Octo¬ 
ber 2010 NGS AO observations of M15 with the high preci¬ 
sion (residual distortion error ~ 0.5 mas) H ST /ACS globu- 
lar cluster survey observation on May 2006 (1 Anderson et ali 
l2008h . A catalogue of the stars (S/N>10) in the NGS AO, 
with positions measured using the IRAF IMEXAM task, 
was matched to the HST catalogue, by using a third or¬ 
der polynomial transformation in X and Y (the IRAF GE¬ 
OMAP task), along with pixel rescaling and rotation; a 
distortion map was obtained, with residual distortion es¬ 
timated to be 7 mas/10 mas (4 mas/6 mas for the 20 
mas detector mode) across t he X and Y axes, respectively. 
iMcNamara, Harrison, fc Andersonl (120031) measured proper 
motions of stars in M15 that would amount to offsets of 
~ 1.3 mas during the time between observations, and can 
therefore be neglected; in addition, these are expected to 
have random vectors, without introducing systematics. A 
larger effect that has not been accounted for is the at¬ 
mospherical refraction, which is not present in the HST 
data and would stretch the image along the direction to the 
zenith. The effect would amount to a total of ~ 20 mas or 
more across the FOV. 

3.7.2 Estimating a more realistic residual distortion 

A more accurate distortion correction was not available at 
the time the bulk of the analysis presented in this paper 
was performed. While the measured residual distortion is 
an overall statistic measured over all pixels across a column 
of the detector, what matters for the purpose of this work 
is a measurement of the distortion that affects relative sep¬ 
arations of objects located only a few arcseconds apart. We 
used two methods to estimate this. 

First, we used a distinct observation of M15, also taken 
in October 2010, with the same detector configuration, in 
the 52 mas mode. We measured the relative separation of 13 
bright star pairs as they are observed at different positions 
on the detector, with a 5-point dither of size 2>". The se¬ 
lected stars do not have nearby companions that may affect 
the measurement, and are not located close to cosmic rays 
or bad detector pixels. The held was distortion-corrected 
beforehand, using the distortion map described above. The 


© 2014 RAS, MNRAS OOO.mi58l 















12 


C.E. Rusu et al. 



remodel using a 
lensing model and 
include the host galaxy 


remodel using a lensing model 


remove B, G; 

and include the host galaxy; 


build a hybrid 

build a hybrid PSF from A and 


PSF from A 

separate PSF & host iteratively 


and remodel 


Figure 7. Flowchart of the overall approach in the morphological modelling of the imaging data. 


Table 2. Technique of morphology fitting 


Object 


Technique 


SDSS J0743+2457 

SDSS J0819+5356 

SDSS J0820+0812 

SDSS J0832+0404 

SDSS J0904+1512 

SDSS J0926+3100 

SDSS J0946+1835 (2012) 

SDSS J0946+1835 (2011) 

SDSS J1001+5027 

SDSS J1002+4449 

SDSS J1054+2733 

SDSS J1055+4628 

SDSS J1131+1915 

SDSS J1206+4332 

SDSS J1216+3529 

SDSS J1254+2235 

SDSS J1313+5151 

SDSS J1320+1644 

SDSS J1322+1052 

SDSS J1330+1810 (2012) 

SDSS J1330+1810 (2013) 

SDSS J1334+3315 ( J/H/K') 

SDSS J1353+1138 

SDSS J1400+3134 

SDSS J1405+0959 {J/H/K') 

SDSS J1406+6126 

SDSS J1455+1447 

SDSS J1515+1511 

SDSS J1620+1203 (2011) 

SDSS J1620+1203 (2012) 


analytical fitting / simulation with star in FOV 

analytical fitting with host galaxy / simulation with star in FOV 

non-AO; fitting with star in FOV 

analytical fitting with host galaxy / simulation with star in the SDSS J0819+5356 FOV 
hybrid PSF created iteratively with host galaxy / simulation with self PSF 
analytical fitting 

hybrid PSF / simulation with the separately observed PSF star 
analytical fitting 

analytical fitting / simulation with self PSF 
analytical fitting / simulation with star in FOV 
analytical fitting 

analytical fitting / simulation with the separate PSF star from SDSS J1405+0959 

hybrid PSF / simulation produced with a cut around image A as PSF, refitted with the hybrid PSF 

analytical fitting / simulation with self PSF 

analytical fitting / simulation with the SDSS J1001+5027 separate PSF star a 
analytical fitting / simulation with the separately observed PSF star 
hybrid PSF / simulation with self PSF 

non-AO image fitting with star in FOV / AO image using image A cut as PSF 
hybrid PSF created iteratively with host galaxy / simulation with self PSF 
analytical fitting / simulation with the separate PSF star from SDSS J1405+0959 
analytical fitting / simulation with the separate PSF star from SDSS J0946+1835 
analytical fitting / simulation with the SDSS J1131+1915 separate PSF star 
hybrid PSF / simulation with self PSF 

analytical fitting / simulation with the SDSS J0743+2457 bright star 

hybrid PSF / simulation with the separately observed PSF star 

non-AO image fitting with analytical PSF / simulation with the analytical PSF 

hybrid PSF / simulation with the SDSS J0832+0404 separate PSF 

fitting with star in FOV 

analytical fitting 

hybrid PSF / simulation with the separate PSF star from SDSS J1405+0959 


a Gaussian convolution was performed on the PSF star to more closely match the PSF of the system. 
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measurements were performed with Galfit, by fitting two- 
dimensional Moffat profiles to the observed stars. The stan¬ 
dard deviations of their measured separations across the X 
and Y axes are typically less than ~ 1 mas, indicating that 
the residual distortion estimated above vastly overestimates 
the errors on the relative separation of close pairs. 

Second, we measured the separations between quasar 
images A and B of the actual AO data reported in this pa¬ 
per, in different dithered frames. This was more difficult, 
as the objects are intrinsically fainter, and there is also con¬ 
tamination from the lensing galaxy, usually located closer to 
the fainter image. To increase S/N, all dither frames where 
the objects fall on the same location on the detector were 
combined, by first calculating, via cross-correlation, neces¬ 
sary sub-pixel offsets due to telescope pointing errors. The 
observed systems were fitted with Hostlens, with the mor¬ 
phological parameters of the galaxies fixed at the values de¬ 
termined from the final, higher S/N science frame. The only 
free parameters were the positions and magnitudes of the 
objects, as well as the analytical 2 Moffat PSF parameters 
(since the PSF changes between frames). This method gave 
usable results (i.e. not affected by components poorly ac¬ 
counted for, or exceedingly large error bars) for 13 objects 
observed in the 52 mas mode, and one in the 20 mas mode. 
The resulting standard deviations are reported in Table [A~3l 

The results of these two methods are compared in Fig¬ 
ure [S] We found that the first method (the M15 star pairs) 
produces overall smaller standard deviations. This could be 
due either to the smaller measurement error bars (high S/N), 
or to the fact that the dither size was smaller, and there¬ 
fore the locations of the pairs on the detector are closer to 
each other. In the latter case, it would be expected that ob¬ 
jects with star pairs/quasar images with larger separations 
would also show comparatively larger standard deviations, 
but this is not obvious from the data. In view of the large 
error bars, the only statistical result that can be drawn is 
the average standard deviation, of ~ 1.5 mas on both axes, 
for the quasar images. Here, two objects with much larger 
values were omitted. This value may be a reflection of the 
large measurement errors, or it may be due to actual resid¬ 
ual distortion which is not corrected by the distortion map. 
As a result, this technique serves as an upper estimate of 
the mean residual distortion. 

For the systems for which the second method gave re¬ 
liable results, we used the measured standard deviations as 
estimates of the residual distortion (its upper estimates). 
The standard deviations are much smaller than the esti¬ 
mated residual distortion across the whole field, from Sec- 
tion l3.7.T1 For the other objects, we used the average value of 
~ 1.5 mas. There are two additional factors that can bias a 
true measurement of the object separation (relative astrom¬ 
etry). First is the uncertainty in the pixel size. This affects 
the object separation as \/~N x cr p i X ei, where N is the number 
of pixels separating the two objects. This effect is very small 
however and can be ignored (Table 1/01) . The second factor is 
the atmospherical refraction. We estimated typi cal numeri- 
cal values for its effect on each object based on lHelminial3 
(l2009ll . which studied the effects of differential atmospheri¬ 
cal refraction on precise astrometry with AO in the infrared. 
We considered the resulting values, typically much less than 
1 mas, as systematics, and reported them in Table IX3l 

To conclude, we found that the residual field distortion 
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Figure 8. Star or quasar image separation on the X and Y 
axis, as a function of its measured standard deviation. The error 
bars are the typical measurement errors in one frame (in the case 
of the observed objects, the average of frames in the same dither 
position). A colour version of this figure is available in the online 
version of the paper. 


affects the relative astrometry at a small level of ~ 1.5 mas. 
We took three different systematics into account, added in 
quadrature and used them as additional errors on the rel¬ 
ative astrometry of each system (including for the lensing 
galaxies located in between the quasar images): scatter in 
the positions at different dither positions, pixel size error, 
and atmospheric refraction. We present further considera¬ 
tions in Section [ 6 ] 


4 LENSING MASS MODELS 

All objects (two-image lenses, in the great majority), unless 
otherwise required, are modelled with the singular isother¬ 
mal sphere with external shear (SIS+ 7 ), singular isother¬ 
mal ellipsoid (SIE), and singular isothermal ellipsoid with 
shear SIE+ 7 ) mass models for the lensing galaxy. For 
the latter model, the lens ellipticity and orientation were 
fixed at the observed values. The error on the flux ratio 
measured in A'—band was increased to 10 %, in order to 
account for microlensing, dust extinction, and time-delay 
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dYonehara. Hirashita. fe Richtedl200i4) . The model s are fit¬ 
ted with the publicly available code Glahc (lOeur 1 120101 ) . 

In cases where error bars are calculated on the best- 
fit model values, the error bars on the magnification factor 
and the time delay, which are not parameters but instead are 
results of lens modelling, are obtained from the histogram of 
magnification/time delay when fixing the model parameters 
to those obtained in the MCMC chains. 

For several objects, the influence of nearby galaxies was 
considered by modelling them with SIS profiles, scaled to 
the mass of the mai n lens from their lumino sity, through the 
Faber-Jackson law dFaber fe Jackson! Il976lj L oc a 4 . This is 
in fact an upper limit to their influence, as it assumes that 
the galaxies and the main lens are at the same redshift, and 
they are early-type. 

A description of the analysis and science results for each 
individual observed system is presented in Appendix iBl and 
0 The astrometry/photometry, lens morphology and mass 
models for each system are provided in Tables [3] [4] and [5j 
respectively. 


5 COMPARISON WITH OTHER AVAILABLE 

DATA 

Here we discuss how the derived astrometry and morphol¬ 
ogy results compare with previous data of the same objects, 
archival or available in the literature. The comparison is de¬ 
tailed for each object in Appendix [B] 

5.1 Comparison with the previously available 
low-resolution data 

The astrometric results obtained from the current AO sam¬ 
ple are generally not in agreement with the previously pub¬ 
lished low-resolution data available, within the formal er¬ 
ror bars. The low-resolution data typically consists of ob¬ 
servations in the optical and I band, performed with the 
UH88/Tek2k (0.22 pixel -1 ) camera. In addition to the 
IRCS+AOI 88 pixel size being ~ 4 times finer (in the 52 
mas mode), the seeing was also ~ 4 times larger when com¬ 
pared to the AO observations, where the typical PSF core 
has FWHM ~ 0.2”. The difference in astrometry is on the 
order of tens of milliarcseconds and larger for the lensing 
galaxies. There are several systems for which the position 
of the lensing galaxy is offset by 1-2 pixels (referring here 
to the UH88/Tek2k pixel size), significantly affecting the 
published mass models (for example, the lens lying clos¬ 
est to the bright quasar image, which is more unusual). In 
all of these cases (SDSS J0832+0404, SDSS J0904+1512, 
SDSS J1001+5027, SDSS J1131+1915, SDSS J1216+3529 
and SDSS J1400+3134), the problem arises from the small 
separation between image B and the lensing galaxy G, of ~1 
pixel. As in all of these cases the peaks of the two objects are 
clearly separated in the AO data, the problem lies unequiv¬ 
ocally with the inherent difficulties in the two-dimensional 
profile fitting multiple components in low-resolution data. 

I 11 the cases where the lens galaxy morphology based 
on the low-resolution data is published, the results are typ¬ 
ically in disagreement within the formal error bars. Again, 
this is most likely attributable to the poor resolution. The 
data do not however require a distinct interpretation, with 


a few exceptions. The main lenses in SDSS J0743+2457 
and SDSS J1405+0959 are fitted in the discovery paper 
with a low and intermediate value Sersic index, respectively, 
whereas the new values in the AO data are larger, indica¬ 
tive of elliptical galaxies. Also, the lens in SDSS J1330+1810 
is originally fitted with Sersic index ~ 3.4, but ~ 2.2 
in the new data, closer to the morphology of a disk-like 
galaxy. In addition, SDSS J1254+2235, SDSS J1353+1138 
and SDSS J1406+6126, originally modelled with a single 
Sersic component, require two concentric components in or¬ 
der to properly fit the AO data. 

5.2 Comparison with HST data 

There is a single object in the current AO sample, 
SDSS J0832+0404, for which archival high-resolution data 
obtained with HST is available (proposal ID 11289, PI. 
Jean-Paul Kneib). Unfortunately, the morphology of the lens 
galaxy is unreliable for a direct comparison with the AO 
data (ignoring also the fact that a galaxy may exhibit dif¬ 
ferent morphology at visible and near-infrared wavelengths), 
due to the galaxy being very faint in the HST data, and in 
close proximity to the much brighter image B. Indeed, the 
morphological model diverges in terms of the Sersic index, 
effective radius and flux. Fixing the Sersic index at n = 4, 
similar to the AO data, leaves residuals, and results in an 
effective radius that is ~ 10 times larger than measured with 
AO. While a relative total flux comparison is problematic, 
the ratio of the peak pixel count at the core of B and G1 is 
~ 15, and it is B/Gl ~ 1 in AO (figure lB6ll . 

On the other hand, the relative astrometry is more 
straightforward to compare. The relative astrometry of A 
and B match at < 1 mas along the Y axis, and ~ 5 mas 
along the X axis in the HST and AO data. This is in good 
agreement considering the large error bars in Table [3] Even 
in case that the AO error bars are overestimated, the fact 
that the discrepancy along the X axis is larger can be ex¬ 
pected due to the quasar host galaxy arc being more elon¬ 
gated along that axis, and affecting the measurement. 
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Table 3: Relative Astrometry and Photometry of the Lens Systems 


15 


Object 

AA' ["] 

AT ["] 

A" [mag] 

K' 0 

n = 

= 4) 


SDSS J0743+2457 (0 

= 1.092" ± 0.020" 

, z s = 2.165, zi = 

= 0.381) 



A 

0.000 ± 0.001 

0.000 ± 0.001 

16.88 ± 0.05 




B 

0.807 ± 0.017 

0.736 ± 0.011 

18.93 ± 0.16 




G 

0.676 ± 0.004 

0.593 ± 0.004 

17.32 ± 0.08 

17.32 

± 

0.06 


SDSS J0819+5356 (0 

= 4.037" ± 0.003" 

, = 2.237, zi = 

= 0.294) 



A 

0.000 ± 0.006 

0.000 ± 0.003 

16.39 ± 0.14“ 




B 

-3.394 ± 0.005 

-2.185 ± 0.006 

18.13 ± 0.14“ 




G1 

-1.991 ± 0.002 

-1.330 ± 0.005 

13.65 ± 0.06 

14.00 

± 

0.01 

G2 

-23.18 

-2.00 

15.07 ± 0.10 





SDSS J0820+0812 (0 

= 2.090" ± 0.030”) 

", z s = 2.024, zi 

= 0.803 



A 

0.000 ± 0.004 

0.000 ± 0.004 

17.21 ± 0.03 




B 

0.918 ± 0.015 

1.878 ± 0.025 

17.98 ± 0.06 




G 

0.690 ± 0.036 

1.252 ± 0.063 

16.81 ± 0.22 





SDSS J0832+0404 (0 

= 2.182" ±0.006” 

, z s = 1.115, zi = 

= 0.659) 



A 

0.000 ± 0.0015 

0.000 ± 0.002 

16.50 ± 0.06“ 




B 

-1.774 ± 0.004 

-1.270 ± 0.003 

17.58 ± 0.06“ 




G 

-1.327 ± 0.004 

-0.947 ± 0.003 

16.33 ± 0.04 





SDSS J0904+1512 (0 = 

= 1.064" ±0.003”, ; 

= 1.826, zi — 

0.2-0.5) 



A 

0.000 ± 0.000 

0.000 ± 0.000 

15.73 ± 0.03“ 




B 

-0.928 ± 0.001 

-0.521 ± 0.002 

17.29 ± 0.03“ 




G1 

-0.570 ± 0.001 

-0.423 ± 0.005 

17.39 ± 0.19 

17.45 

± 

0.03 

G2 

0.912 ± 0.017 

-4.967 ± 0.031 

18.97 ± 0.09 





SDSS J0946+1835 (0 

= 3.055 ±0.002”, 

= 4.799, zi = 

0.388) c 



A 

0.000 ± 0.000 

0.000 ± 0.000 

17.10 ± 0.00 




B 

1.754 ± 0.000 

-2.502 ± 0.000 

17.57 ± 0.00 




G1 

1.459 ± 0.000 

-1.927 ± 0.001 

15.00 ± 0.03 

15.29 

± 

0.02 

G2 d 

2.307 ± 0.008 

2.751 ± 0.006 

20.03 ± 0.05 




G3 

-1.270 ± 0.007 

0.080 ± 0.013 

21.16 ± 0.20 





SDSS J1001+5027 (0 = 

= 2.928" ± 0.0015" 

', z s = 1.839, zi 

= 0.415) 



A 

0.000 ± 0.000 

0.000 ± 0.001 

15.84 ± 0.03 




B 

2.489 ± 0.001 

1.543 ± 0.001 

16.10 ± 0.04 




G1 

2.054 ± 0.012 

1.122 ± 0.009 

15.87 ± 0.11 

16.12 

± 

0.03 

G2 

1.799 ± 0.003 

-0.865 ± 0.005 

17.55 ± 0.22 





SDSS J1055+4628 (0 

= 1.145" ±0.003" 

, z s = 1.249, zi = 

= 0.388) 



A 

0.000 ± 0.001 

0.000 ± 0.001 

17.10 ± 0.03 




B 

-0.043 ± 0.002 

1.145 ± 0.002 

18.31 ± 0.06 




G1 

-0.048 ± 0.002 

0.850 ± 0.002 

16.68 ± 0.11 

16.60 

± 

0.15 

G2 d 

-12.03 

-6.05 

16.53 ± 0.02 





SDSS J1131+1915 (0 = 1.403" ± 0.002' 

", z s = 2.915, zi 

- 0.3) 



A 

0.000 ± 0.000 

0.000 ± 0.000 

15.97 ± 0.00 




B 

-1.396 ± 0.000 

-0.134 ± 0.000 

17.55 ± 0.00 




G1 

-0.946 ± 0.003 

-0.142 ± 0.003 

17.90 ± 0.09 

18.10± 

0.03 

G2 d 

-8.72 

-3.85 

17.94 ± 0.03 





SDSS J1206+4332 (0 

= 3.026" ±0.005", 

Zs = 1.789, zi = 

= 0.748)® 



A 

0.000 ± 0.003 

0.000 ± 0.002 

16.56 ± 0.02 




B 

-0.035 ± 0.003 

3.027 ± 0.003 

16.60 ± 0.03 




G1 

-0.508 ± 0.014 

1.813 ± 0.012 

15.85 ± 0.45 




G2 

1.120 ± 0.029 

6.070 ± 0.034 

16.80 ± 0.76 




G3 

-1.853 ± 0.014 

2.470 ± 0.060 

17.44 ± 0.06 





© 2014 RAS, MNRAS OOO.mt58l 























16 C.E. Rusu et al. 


G4 

1.70 ± 0.11 

4.75 ± 0.12 

19.48 ± 0.81 


G5 

-7.23 ± 0.10 

-0.27 ± 0.11 

18.38 ± 0.05 



SDSS J1216+3529 (i 9 

= 1.487" ± 0.013" 

, 2 , = 2.012, zi ~ 

0.55) 

A 

0.000 ± 0.002 

0.000 ± 0.002 

16.82 ± 0.03 


B 

1.484 ± 0.011 

0.097 ± 0.006 

17.87 ± 0.09 


G1 

1.165 ± 0.012 

0.013 ± 0.014 

17.37 ± 0.13 

17.43 ± 0.10 

G2 d 

~7.51 

-12.63 

17.46 ± 0.04 


G3 

1.781 ± 0.051 

2.803 ± 0.066 

18.47 ± 0.81 



SDSS J1254+2235 {9 = 1.561" ± 0.004" 

z 3 = 3.626, zi - 

- 0.3) 

A 

0.000 ± 0.001 

0.000 ± 0.002 

16.95 ± 0.06 


B 

-1.446 ± 0.003 

-0.588 ± 0.002 

17.82 ± 0.05 


Gla 

-0.871 ± 0.004 

-0.317 ± 0.003 

17.58 ± 0.12 

17.40 ± 0.15 

Gib 

-0.871 ± 0.004 

-0.317 ± 0.003 

17.54 ± 0.11 

17.71 ± 0.17 / 

Gla+Glb 

-0.871 ± 0.004 

-0.317 ± 0.003 

16.80 ± 0.03 


G2 

-2.456 ± 0.022 

-0.458 ± 0.025 

21.01 ± 0.24 



SDSS J1313+5151 (0 

= 1.223" ±0.001", 

z s = 1.875, zi = 

0.194) 

A 

0.000 ± 0.000 

0.000 ± 0.000 

14.62 ± 0.02 


B 

0.309 ± 0.000 

-1.184 ± 0.000 

16.07 ± 0.01 


G 

0.255 ± 0.001 

-0.977 ± 0.002 

15.39 ± 0.03 

15.69± 0.03 


SDSS J1322+1052 (0 

= 1.999" ± 0.002" 

, z 3 = 1.716, zi ~ 

0.55) 

A 

0.000 ± 0.000 

0.000 ± 0.000 

16.22 ± 0.06“ 


B 

1.070 ± 0.001 

1.689 ± 0.001 

17.80 ± 0.06“ 


G1 

0.738 ± 0.0025 

1.187 ± 0.001 

15.24 (?) 

16.70 ± 0.02 

G2 

-6.27 

~9.92 

18.81 ± 0.02 


G3 

-2.672 ± 0.051 

1.443 ± 0.020 

19.28 ± 0.14 


G4 

~5.94 

~4.31 

19.43 ± 0.08 


G5 

~1.09 

~11.94 

19.51 ± 0.02 


SDSS J1330+1810 (max. 

sep. 1.710" ± 0.003", = 1.393, z t 

= 0.373) 

A 

0.000 ± 0.001 

0.000 ± 0.000 

16.72 ± 0.02 


B 

0.414 ± 0.0015 

-0.011 ± 0.001 

17.15 ± 0.03 


C 

1.248 ± 0.002 

1.168 ± 0.0015 

17.61 ± 0.03 


D 

-0.244 ± 0.002 

1.579 ± 0.002 

18.88 ± 0.02 


G 

0.221 ± 0.008 

0.966 ± 0.0015 

15.31 ± 0.06 



SDSS J1334+3315 9 (0.838" ± 0.002", 2 

■ s = 2.426, Zi = 0.557) 

A 

0.001 ± 0.001 

0.000 ± 0.000 

18.74 ± 0.03 J 
18.31 ± 0.02 H 
17.67 ± 0.06 K' 


B 

-0.491 ± 0.001 

0.680 ± 0.001 

19.07 ± 0.04 J 
18.63 ± 0.03 H 
17.97 ± 0.07 K' 


G 

-0.134 ± 0.009 

0.220 ± 0.009 

20.53 ± 0.36 J 

19.93 ± 0.35 H 




18.88 ± 0.38 K' 


G1 

-3.171 ± 0.001 

-2.609 ± 0.001 

18.74 ± 0.01 J 
18.09 ± 0.02 H 
17.56 ± 0.01 K' 


G2 

1.151 ± 0.009 

2.049 ± 0.009 

23.11 ± 0.26 J 
20.81 ± 0.10 H 
20.31 ± 0.07 K' 



SDSS J1353+1138 {0 

= 1.392" ±0.0015” 

, z 3 = 1.629, zi r . 

- 0.25) 

A 

0.000 ± 0.000 

0.000 ± 0.000 

15.02 ± 0.01 


B 

0.245 ± 0.000 

-1.371 ± 0.000 

15.98 ± 0.00 


Gla 

0.248 ± 0.001 

-1.042 ± 0.000 

15.46 ± 0.02 

16.00 ± 0.20 

Gib 

[0.248 ± 0.001] [ 

-1.042 ± 0.000] 

17.13 ± 0.05 

16.20 ± O^ 

Gla+Glb 

[0.248 ± 0.001] 

-1.042 ± 0.000] 

15.25 ± 0.01 

15.36 ± 0.01 
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G2 

2.824 ± 0.009 

5.392 ± 0.005 

18.82 ± 0.09 


SDSS J1400+3134 (0 = 1.744" ± 0.003' 

', = 3.317, zi ~ 0.8) 

A 

0.000 ± 0.000 

0.000 ± 0.000 

16.94 ± 0.03 

B 

1.032 ± 0.001 

1.406 ± 0.001 

17.42 ± 0.03 

G1 

0.327 ± 0.005 

1.025 ± 0.004 

17.81 ± 0.14 

G2 d 

-6.37 

~5.78 

18.70 ± 0.05 


SDSS J1406+6126 (0 = 2.029" ± 0.011", 

= 2.134, zi = 0.271)' 1 

A 

0.000 ± 0.006 

0.000 ± 0.002 

17.43 ± 0.02 

B 

1.686 ± 0.006 

-1.129 ± 0.006 

17.99 ± 0.03 

Ga 

1.160 ± 0.002 

-0.750 ± 0.003 

16.53 ± 0.13 

Gb 

1.160 ± 0.002 

-0.750 ± 0.003 

16.07 ± 0.11 

Ga+Gb 

1.160 ± 0.002 

-0.750 ± 0.003 

15.52 ± 0.12 


SDSS J1455+1447 (0 = 1.753" ± 0.001", 

, 2 S = 1.424, zi ~ 0.42)" 

A 

0.000 ± 0.000 

0.000 ± 0.000 

15.81 ± 0.00 

B 

1.609 ± 0.000 

0.696 ± 0.000 

16.95 ± 0.01 

G1 

1.003 ± 0.005 

0.538 ± 0.004 

16.05 ± 0.10 16.37 ± 0.02 

G2 

~10.75 

~2.10 

18.03 ± 0.03 


SDSS J1515+1511 (0 = 2.010" ± 0.001”, 

, z 3 = 2.054, zi = 0.742) 

A 

0.000 ± 0.001 

0.000 ± 0.000 

16.01 ± 0.01 

B 

1.668 ± 0.001 

-1.122 ± 0.001 

16.35 ± 0.01 

G1 

1.357 ± 0.003 

-0.848 ± 0.003 

17.29 ± 0.02 

G2 

13.49 —8.30 

17.41 ± 0.02 d 


SDSS J1620+1203 {0 = 2.807" ± 0.004", 

z B = 1.158, zi = 0.398)" 

A 

0.000 ± 0.000 

0.000 ± 0.000 

17.02 ± 0.01 

B 

-2.107 ± 0.001 

-1.856 ± 0.001 

18.56 ± 0.02 

G 

-1.764 ± 0.001 

-1.458 ± 0.001 

15.52 ± 0.02 15.72 ± 0.01 


Magnitudes are in the Vega system. All magnitudes are corrected for atmospheric and Galactic extinction. The positive directions 
of X and Y are West and North, respectively. The objects without attached error bars on the positions are located at a relatively 
large distance from the target, where the residual distortion is expected to be large. Here n = 4 refers to the Sersic index, and zi , z s 
are the redshifts of the main lens galaxy and source quasar, respectively; where Zi is approximate, its value is a photometric redshift 
estimate. For the objects where the quasar host galaxy has been explicitly modelled, relative photometry of the images includes the 
contribution of the point source as well as the host. The results for SDSS J1405+0959, presented in a separate publication, are not 
included (see Appendix IB 191) . 

a The fluxes are estimated with Hostlens for the point source, and reflect the uncertainty in the host galaxy - quasar disentangling, 
therefore attached errors should not be used to estimate the error on the flux ratio, which are smaller. 
b The zero-point error is 0.14 mag for both SDSS J0820+0812 and SDSS J1620+1203. 

c In addition to the quoted errors on photometry, a systematic error of 0.05 mag was obtained during simulations, which should be 
added to the absolute photometry, but does not affect the relative values. 

d Aperture photometry is given, and the positions and their errors are measured with the PHOT task. 
e The astrometry for SDSS J1206+4332 should be considered less reliable (see Appendix IB 101) 
f Sersic index of Gla (not Gib) is fixed at 4. 

9 The astrometry represents the average over the H and K' bands, and the magnitudes of G1 and G2 were obtained from aperture 
photometry. 

h Absolute photometry unreliable due to non-photometric conditions. 


Table 4: Morphological parameters 


Object 

e 

PA 

r-e n 

n 

SDSS J0743+2457 

G 

0.80 ± 0.12 

23.8 ± 2.0 

0.19 ± 0.05 

4.13 ± 0.92 

G 

0.79 ± 0.12 

23.8 ± 2.0 

0.19 ± 0.05 

[4] 

SDSS J0819+5356 


© 2014 RAS, MNRAS 000.111 f58l 



















18 C.E. Rusu et al. 


G1 

0.25 ± 0.03 

-35.5 ± 2.4 

2.89 ± 0.45 

6.28 ± 0.34 

G1 

0.15 ± 0.01 

-31.8 ± 1.2 

1.22 ± 0.05 

[4] 

SDSS J0820+0812 

G“ 

0.25 ± 0.05 

-44.8 ± 5.9 

1.58 ± 0.30 

1.37 ± 0.75 

SDSS J0832+0404 

G 

0.21 ± 0.02 

9.0 ± 8.3 

0.26 ± 0.05 

3.85 ± 0.70 

SDSS J0904+1512 

G1 

[0] 

_ 

0.32 ± 0.09 

4.89 ± 0.71 

G1 

[0] 

- 

0.29 ± 0.01 

[4] 

G2 

0.21 ± 0.12 

43.0 ± 11.7 

0.46 ± 0.05 

0.56 ± 0.21 

SDSS J0946+1835 

G1 

0.09 ± 0.01 

-11.9 ± 3.9 

1.37 ± 0.09 

5.64 ± 0.30 

G1 

0.09 ± 0.00 

-11.6 ± 3.6 

0.79 ± 0.04 

[4] 

G2 

0.34 ± 0.04 

68.7 ± 4.4 

0.30 ± 0.01 

0.78 ± 0.11 

G3 

0.66 ± 0.09 

-9.6 ± 4.1 

0.21 ± 0.07 

0.27 ± - 0.5 

SDSS J1001+5027 

G1 

0.10 ± 0.08 

8.7 ± 5.8 

1.77 ± 0.36 

5.74 ± 0.56 

G1 

0.10 ± 0.01 

1.0 ± 7.31 

1.06 ± 0.04 

[4] 

G2 

0.16 ± 0.03 

-43.9 ± 10.1 

0.41 ± 0.13 

3.23 ± 1.80 

SDSS J1055+4628 

G1 

0.36 ± 0.04 

-1.7 ± 3.3 

0.48 ± 0.30 

5.83 ± 2.40 

G1 

0.41 ± 0.05 

-0.5 ± 2.2 

0.39 ± 0.05 

[4] 

SDSS J1131+1915 

G1 

0.47 ± 0.07 

87.5 ± 5.2 

0.43 ± 0.11 

6.58 ± 1.14 

G1 

0.45 ± 0.07 

87.8 ± 4.9 

0.27 ± 0.02 

[4] 

SDSS J1216+3529 

G1 

0.67 ± 0.11 

-14.60 ± 4.10 

0.32 ± 0.10 

5.02 ± 2.83 

G1 

0.68 ± 0.11 

-15.20 ± 2.20 

0.28 ± 0.05 

[4] 

SDSS J1254+2235 

Gla 

0.32 ± 0.05 

-61.3 ± 2.7 

0.14 ± 0.03 

2.6 ± 1.6 

Gib 

0.28 ± 0.05 

-26.2 ± 3.9 

1.09 ± 0.06 

0.23 ± 0.07 

Gla 

0.31 ± 0.05 

-60.8 ± 2.6 

0.72 ± 0.19 

[4] 

Gib 

0.30 ± 0.05 

-23.8 ± 3.1 

1.13 ± 0.05 

0.17 ± 0.09 

SDSS J1313+5151 

G 

0.35 ± 0.01 

-68.7 ± 1.0 

1.32 ± 0.09 

7.0 ± 0.43 

G 

0.32 ± 0.01 

-68.8 ± 1.2 

0.69 ± 0.03 

[4] 

SDSS J1322+1052 

G1 

0.07 ± 0.03 

-52.3 ± 7.8 

0.48 ± 0.02 

[4] 

G3 

0.59 ± 0.09 

72.1 ± 5.7 

0.59 ± 0.28 

0.94 ± 0.30 

SDSS J1330+1810 

G 

0.70 ± 0.03 

25.5 ± 0.4 

0.49 ± 0.03 

2.23 ± 0.12 

SDSS J1334+3315 

G (J) 

unreliable 

- -21 ± 22 

0.25 ± 0.05 

2.53 ± 1.36 

G (H) 

0.61 ± 0.29 

~ - 23 ± 5.4 

0.22 ± 0.11 

[4] 

G ( K ') 

0.43 ± 0.23 

-36.5 ± 18.7 

0.34 ± 0.17 

3.4 ± 2.1 

G ( K' host) 

0.62 ± 0.23 

-5.0 ± 3.8 

0.18 ± 0.05 

~2.6 ± 0.81 

G1 (J) 

0.44 ± 0.02 

-57.4 ± 1.4 

0.78 ± 0.03 

1.51 ± 0.12 

G1 {H) 

0.55 ± 0.03 

-60.5 ± 1.7 

0.73 ± 0.16 

2.58 ± 0.54 


© 2014 RAS, MNRAS OOO.HH581 





























Subaru Telescope adaptive optics observations of gravitationally lensed quasars in the SDSS 


19 


G1 (K') 

0.60 ± 0.02 

-59.5 ± 1.0 

0.73 ± 0.13 

3.40 ± 0.42 

SDSS J1353+1138 

Gla 

0.50 ± 0.08 

-64.8 ± 0.9 

0.47 ± 0.01 

5.77 ± 0.19 

Gib 

0.69 ± 0.03 

-66.8 ± 0.9 

1.04 ± 0.03 

0.32 ± 0.08 

Gla 

0.52 ± 0.04 

-63.9 ± 1.4 

0.20 ± 0.05 

[4] 

Gib 

0.61 ± 0.03 

-66.6 ± 0.6 

1.04 ± 0.07 

0.88 ± 0.25 

SDSS J1400+3134 

G1 

0.32 ± 0.05 

-66.6 ± 7.2 

0.31 ± 0.09 

[4] 

SDSS J1406+6126 

Gla 

[0] 

— 

0.13 ± 0.04 

[4] 

Gib 

0.69 ± 0.01 

-28.0 ± 0.60 

1.14 ± 0.06 

[1] 

SDSS J1455+1447 

G1 

0.23 ± 0.02 

82.1 ± 5.3 

0.87 ± 0.17 

8.40 ± 1.25 

G1 

0.26 ± 0.02 

86.3 ± 2.5 

0.45 ± 0.05 

[4] 

SDSS J1515+1511 

G1 

0.81 ± 0.03 

-17.1 ± 0.3 

0.66 ± 0.01 

1.65 ± 0.15 

SDSS J1620+1203 

G 

0.24 ± 0.02 

-27.2 ± 1.2 

1.30 ± 0.07 

4.88 ± 0.20 

G 

0.22 ± 0.01 

-24.7 ± 0.8 

1.02 ± 0.01 

[4] 


The position angle is positive from North towards East. 
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Table 5: Lensing models 


bO 

o 


Model 

a [km/s] (M [h L M+) / 0 Bin 

e,9 e [deg] 

7 ,© [deg] 

At [days] 



SDSS J0743+2457 (z s = 2.165, z t = 

0.381) 



SIS +7 

165.2+1/0.58+;!! 

/ 

0.063++, -31.7+2 

23.ll!;! 

q co+ 0.22 

0.90-045 

SIE 

163.3+1/0.56+;!! 

0.17+!®, -31.9+2 

/ 

22 . 0 l°! 

0 cc+0.18 
o.oo_047 

SIE +7 

187.9+ 3 3 6 /0.74++ 

(0.80+; ( j , ?,23.8+;!) 

0.496++!, —63.4+;! 

30.91!;! 

4 45+0-27 

^•^0-0.91 


SDSS J0819+5356 (z s = 2.237, z t = 

0.294) 



SIS +7 

304.5++/2.12++ 

/ 

0 f)49 ”*~ 0 ' 003 —60 0 - * -0 ' 3 
u.u^y _ 0 003 , ou.u _ 0 4 

66.41+ 

16.951!;!! 

SIE 

301.2±g;| /2.07±g;gS 

0.13+°!, —60.0+1 

/ 

63.4+;! 

16.901+! 

SIE +7 

300.2+72.06+°! 

(0.25+ °!, -35.5+;!) 

0.079+;!!!, 69.91!;! 

61 5 + 1 - 2 

O 1 .O- 0.6 

1 7 99 + + O 8 
± (.zz o gg 


SDSS J0820+0812 ( z 3 = 2.024, z t = 

0.803) 



SIS +7 

277.3±3;g/1.04±g;“ 

/ 

0 n9Q + 0 - 009 95 q + 13 - 7 
u.uzy _ 0 021 , ^o.y_ 2 o.i 

1 99 7+21-9 

100 .1 _io .6 

5.39++ 

SIE 

277.8+1/1.05+°! 

0.08+;°!, 24.2+!;° 

/ 

135.21!+ 

5 95+0-30 

O.OO-Q gl 

SIE +7 

272.3+l/1.00+;° 3 

(0.25+;°°, -44.8+;®) 

0 122 4-0 ' 009 40 4”*" 8 - 6 
u.izz _ 0 046 , 4 tu. 4_ 6 3 

121.7+ 5 3 ® 

5.56+; 3 ° 


SDSS J0832+0404 (z s = 1.115, zi = 

0.659) 



SIS +7 

332.6±+/1.08+;°! 

/ 

n m k+ 0-000 99 0 + 6.0 
u.u±o _ 0 013 , oz.o — 12 .0 

233.0+1 

9 Q 4+ 0 - 10 
o.y4t_oo2 

SIE 

333.9+;!/1.08+;gJ 

0.04+;°!, 32.8+;! 

/ 

236.4+'! 

3.951°;!° 

SIE +7 

337.2+;7l.l0+°! 

(0.21+;°!, 9.0+;l) 

0.068+;°+ 95.1+;! 

244.2+;! 

3.96l°;°° 


SDSS J0904+1512 ( z s 

= 1.826, zi r**> 0.2 — 0.5, zi = 0.3 used for modelling) 



SIS +7 

161.3±g;U/0.57+;g 

/ 

0.082+ °!!, -85.6+;! 

8 . 21 + 

7 QH + 0 - 46 

' - yu 0.39 


SDSS J0946+1835 (z„ = 4.799, zi = 

0.388) 



SIS +7 

9+2.8 /-x 0-1 +0.03 

zo/ .z_ x g/ i.oi_q 02 

/ 

0 1 52”*”°' 011 —46 Q ~*~ 1-4 
u.ioz_ 0018 , z to.y _ 1 0 

121 . 61 + 

q 9 /A+O.O 2 

O.ZU -0 06 

SIE 

249.4+;l 4 /1.45+;°! 

0.38++, -46.7+; 3 5 

/ 

143.2+;! 

9 1 9 + 0-01 
°-- lo -0.04 

SIE +7 

237.4+;®/1.31+;°? 

(0.09+;°!, -11.9+; 3 ) 

0 145+ 0 010 _ 5 i q + 10 
u.I9:0_oo23’ O1 - y -0.06 

122.3+;! 

9 1 q+ 0-02 

O.1O-006 


SDSS J1001+5027 {z s = 1.839, z t = 

0.415) 



Glonly; Power law+ 7 ; 7 *= 1 .72++ 
Glonly; SIE +7 + k; k = 0.18++ 
Glonly; SIE +7 (h free) h = 0.82+;!! 
Gl/G2=ct. Power law +7 (Gl) + SIS (G2) 

/1.44±g;g« 

225.1+;°/1.00+;° 3 

247.9+;°/1.21+;° 3 

/0.96+°® 

(0.10+;°®, 8.7+;®) 
(0.10+;!!, 8.7+;°) 
(0.10+82, 8.7+;®) 
(0.10+;^, 8.7+;°) 

0.050:!;!°!, 18.2+®;! 
0.150+;°!!, 24.6+;! 
0.181+;°+ 24.61!;° 
0.170+;!!!, 28.81+ 

(119.3+;!) 

(119.3+;!) 

(119.3+;!) 

(119.3+;!) 

6 97+ 120 
u -°' -0.04 

4 60"*” 0,38 

4t.oy _ 0 34 

9 1 q+ 0-04 
o.iy_ 0 .i3 

9 90 +O.H 

O.OO-Q so 

ry> — O 1 2 + 0,30 

... 7 —z.iz _ 0 16 

Gl/G2=ct. SIE +7 + k (Gl) + SIS (G2) 

231.8+°/1.05+;°® 

(o.io+;82,8-7+:!) 

0 127”*~ 9 ' 929 20 4 T 3 -^ 
u.iz^_o 0375 ^ y -^_3.0 

(119.3+;°) 

9 5«+0 32 

0.0 U-o.28 

... « = -0.06++ 

G2 free; SIE +7 (Gl) + SIS (G2) 

230.1+;7l.04+;° 3 

(0.10+;!!, 8.7+;!) 

n 1 91 +0.025 qo r\+2.9 

U. 131_o. 037 , z 8 .U _5 6 

(119.3+;!) 

q rro + 0.16 

' °-0.13 

■■■ crWo-gi = 0.59+H 






Glonly; Power law+ 7 ; 7 ' = 1.76jlo q 8 

/l 9 Q + 0-05 

/ 1 .00 — Q.Q4 

/ 

0 104+0- 008 16 4 - *" 2-4 

U.lU9:_o 0335 iD - 4 -3.1 

(119.3+;!) 

5 68"*" 1,26 

o.oo_o 00 

Glonly; SIS +7 ( h free) h = 0.79+'!? 

245.7+;7l.l9+;°7 

/ 

0 210"*" 0 ' 016 22 6”*" 10 

U.Z1U —0 020 ’ zz -°- 1.2 

(119.3+;!) 

q qq + 0.06 

O.ZO-Q 12 

Glonly; SIS +7 + k; k = 0.15±g g| 

226.9+2/1.01+33 

/ 

0 170”*" 0 016 22 6 + 0.8 
u.i/y_0 .029’ zz -°-i.o 

(119.3+;!) 

4 44 ”*" 9,32 
^•^-0.28 

Gl/G2=ct. Power law +7 (Gl) + SIS (G2) 

/o.9i±g;gJ 

/ 

n 99/1 +0.022 ry-r q + 0.7 
U.Z09 :—0 085 ’ Z ' ,Z - 1.0 

(119.3+1) 

q qq + 0.00 

O.OO-Q 64 

— 9 24+ 0 10 
... 7 — z.z 4_ 0 09 


Gl/G2=ct. 2 SIS +7 + re; k = -0.09+;!? 

233.8+;°/1.07+;°? 

/ 

Q-tcr+0.0i9 26 1+ 9 ' 9 

U.1OO-0 023’ ZD i -l.l 

(H9.3+1) 

9 9c;+0-19 

O.OO-q 24 
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G2 free; 2 SIS+ 7 ; o^/cgi = 0.54l?'8| 

231.2±| ; ^/1.05±g;S 

/ 

0.1631? ???, 24.61?-? 

(119.3! 2 ;?) 

3.69!?;?? 

Gl/G2=ct. 2 SIS +7 (h free) h = 0.62±g;g| 

/0.98tg;« 

/ 

0 142~*”°' 010 26 I 4 " 1-1 

U.14tz_ 0 029’ ZD * A —1.1 

(119.31?-°) 

3.99!?;?? 


SDSS J1055+4628 {z„ = 1.249, z t = 

0.388) 



SIS +7 

180.3±g-g /0.58±g-g? 

/ 

0 020 ~*~ 0 ' 005 —35 4+ 16 - 1 
u.uzu_Q 003 ’ °°*^_9.8 

99 O+0.3 
ZZ-Z -0.3 

4 16"*"°' 10 
^• ±D -o.io 

SIE 

180.1±0;®/0.58l°;“ 

0.067°;, -35.0±ii;g 

/ 

99 1 + 6-2 
ZZ-± - 0.1 

4 15+ 0 ' 08 
^• lo _ 0.12 

SIE +7 

173.8!?;?/0.54!g;?? 

(0-357;°°, -1-717) 

0 144+°- uls —87 

0.035 ’ °' - °-3.6 

18.9!?;® 

4 24”*”°' 12 

SDSS J1131+1915 {z s = 2.915, zi ~ 0.3) 

SIS +7 

178.ll?;?/0.74l°;?! 

/ 

0 . 0661 ° °°?, 80 . 81 ?;? 

1 cc r + 0.2 
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The position angle is positive from North towards East. Values inside parentheses are free to vary, but with a gaussian prior given by the observed uncertainties. For the lenses where a spectroscopic lens 
redshift is not available (SDSS J0904+1512, SDSS J1131+1915, SDSS J1216+3529, SDSS J1254+2235, SDSS J1322+1052, SDSS J1353+1138 and SDSS J1400+3134), the effect of the redshift uncertainty on 
the estimated velocity dispersion and time delays is not included; this can be calculated by assuming typical redshift uncertainties of ~ 0.1 in the discovery papers. Unless otherwise specified, all models have 
(x^/d.o.f. = 0/0). The nearby galaxy companions are ignored in these models. G1/G2 refers to the velocity dispersion ratio. Notation p, represents the total magnification factor; S and -y / represent the third order 
perturbation and the density slope. The results for SDSS J1405+0959, presented in a separate publication, are not included (see Appendix |B19l ). For SDSS J1330+1810, as well as all other systems with nearby 
galaxies, these models do not consider explicitly the nearby galaxies. For SDSS J0904+1512, SDSS J1206+4332, SDSS J1254+2235, SDSS J1334+3315 and SDSS J1406+6126, where the ellipticity of the lensing 
galaxy is either negligible or unreliable, the SIE+'y model was omitted. The systems where a close companion is identified can be inferred from Tableful and those with a noticeable field overdensity in SDSS from 
Figure [B3l 
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6 DISCUSSION OF THE TECHNIQUE 

Here we take a look back on the morphological modelling 
technique we employed, by summarising several points: 

• For the majority of the systems, most of the error bars 
on the fitted parameters come from the simulations with 
either the PSF built on the bright quasar image or a sepa¬ 
rately observed PSF star. This is expected since these sim¬ 
ulations are designed to test for systematic errors. These 
PSFs include components that cannot be modelled analyti¬ 
cally, the noise frames used in the simulations may include 
faint background sources, and in the case that the PSF is 
built on the bright image, more noise is introduced, boost¬ 
ing the error estimate. In the cases where the systems are 
bright, it can be argued that the use of noise frame is not 
realistic because they do not contain Poisson noise, which is 
expected to be dominant. In regards to this, the use of noise 
maps produced by Hostlens/Galfit includes an estimate of 
Poisson noise, which was therefore taken into account. 

• While the separately observed PSF star that is most 
similar to the PSF of each system in terms of ana¬ 
lytic parameters was used for simulations, we tested for 
SDSS ,11620+1203 that the use of a second PSF star (from 
SDSS J0832+0404) would produce similar error bars, in case 
we model this system with an analytical PSF. Therefore the 
fact that we simulate our systems with different but similar 
AO PSFs does not appear to have a large effect. 

• The analytical PSF method can lead to spurious detec¬ 
tions and mistaken interpretations, if it is not cross-checked. 
This has been seen in the case of SDSS J1002+4449 and 
SDSS J1353+1138. For SDSS J1002+4449, a comparison 
of the \ 2 values indicates that a lensing galaxy is present, 
whereas simulations indicate that noise is most likely mis¬ 
taken for a galaxy. In the case of SDSS J1353+1138, mod¬ 
elling with an analytical PSF suggests the detection of an 
additional stellar component, whereas modelling with a hy¬ 
brid PSF shows that, once again, noise was mistaken for 
another component. 

• Inconsistencies in astrometry are found in the multi¬ 
ple observations of extended objects in SDSS J1405+0959 
and SDSS J1330+1810. While the scatter between observa¬ 
tions has been taken into account in order to infer realistic 
error bars for these particular systems, the inconsistencies 
do not cast doubt on the astrometric precision estimated 
for the other objects in the current sample, because there 
are known sources of systematics for the two objects. For 
SDSS J1405+0959, the low S/N and the uncertain mor¬ 
phology of GX, located unusually close to G2, is a source 
of systematics. For SDSS J1330+1810, the residuals show 
that there are systematics in the modelling of the PSF (20 
mas mode). However the astrometry of point sources in these 
two systems, as well as for SDSS J1334+3315, with typical 
error bars of ~ 1 mas, is consistent between observations, 
supporting the astrometric precision estimation. 

• As is expected because the final science frame has larger 
S/N and the residual astrometric distortion estimate was 
likely dominated by measurement errors, the systematic dis¬ 
tortion estimate typically dominates the astrometry errors. 
These errors would be further reduced by using a more pre¬ 
cise distortion map and applying corrections, not just error 
estimates, for differential refraction. 

• Most of the galaxies are not consistent with a Sersic 


index n = 4 (de Vaucouleurs profile). They show more pro¬ 
nounced residuals when the de Vaucouleurs profile is used. 
Instead, most galaxies are fitted with a higher value of 
n = 5 — 6. This is in agreem ent w ith the mean index value 
for the lens galaxies found in lBolton et al.l (l2012ll . 

• For SDSS J1515+1511, the only object in the current 
sample for which a suitable PSF star exists in the FOV, 
masking the central 2x2 pixels (in order to reject pixels 
in both the PSF and the target that may be dominated by 
Poisson noise) has introduced significant changes in the pa¬ 
rameters of the host in the sense of increasing the effective 
radius and Sersic index, and reducing its flux; on the other 
hand, for SDSS J0819+5356 (modelled with an analytical 
PSF) the changes are insignificant. The effect could not be 
checked for SDSS J0904+1512, which was modelled itera¬ 
tively with a hybrid PSF built on the target. 

• A matter of concern is inferring the lens galaxy / host 
galaxy parameters in shallow images where the outer re¬ 
gions of the galaxies may be poorly detected. However, the 
simulations performed should in principle account for this 
effect. A more robust method is to compare deep and shal¬ 
low exposures of the same object, which was done for the 
lens galaxies in SDSS J0946+1835 and SDSS J1330+1810, 
revealing that the morphological parameters match within 
~ 1 — 2 a. This however was not done for the host galaxies: 
for SDSS J1330+1810, where both deep and shallow ob¬ 
servations are availa ble, the ho st is not fit satisfactorily in 
either observation. IPeng et al.l (l2006h undertook this study 
and found differences in the estimated host magnitude and 
effective radius of ~ 0.15 mag and 10% - 20%, respectively, 
comparable to the effects they found from using different 
PSFs. 


7 QUASAR HOST GALAXIES IN THE AO 

SAMPLE 

Here we discuss the properties of the detected quasar host 
galaxies, and carry out a study on whether their lumi¬ 
nosities follow the relations found in the literature for 
quasars at similar redshifts. As shown in Table [S] we mea¬ 
sured the quasar host galaxy luminosity in four systems 
(SDSS J0819+5356, SDSS J0832+0404, SDSS J0904+1512 
and SDSS J1515+1511), and found upper limits found for 
an additional system (SDSS J1322+1052). 

When modelling the host galaxies, we considered a sin¬ 
gular isothermal profile for the lensing galaxy. Since the real 
lenses may not be isothermal , this affects th e infer red size 
and luminosity of the hosts. iMarshall et ali (l2007l l trans¬ 
late the uncertainty in the mass profile slope cry into addi¬ 
tional uncertainty on the size and luminosity as cry/r oc cry, 
where r stands for the galaxy radius, and cr mAB oc 2.2cry, 
respectively. Here cry can be taken to be the intrinsic sprea d 
of power-law induced measured bv lKoopmans et al.l (|2009h . 
cry ~ 0.20. This dominates the error budget for the majority 
of the systems. 

All four hosts with available measurements have 
small effective radii of 2.2, 3.2, 0.4 and 2.1 kpc for 
SDSS J0819+5356, SDSS J0832+0404, SDSS J0904+1512 
and SDSS J1515+1511, respectively. It is a known fact that 
quiescent galaxies at high redshift z ~ 2 are m ore compact 
than galaxies of comparable mass at z ~ 0 fe.g.. lDaddi et all 
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120051 : FlVuiillo et ahll2007l ). Also. lPeng et al.l (120061 ) measured 
for their sample of quasar hosts effective radii <3 — 5 kpc, 
which are not much larger than obtained here (with the ex¬ 
ception of SDSS J0904+1512). 


All four hosts have small Sersic indices indicative of disk 
galaxies. At low redshifts, bulge-dominated hosts are more 
numerous, especially for the more luminous quasa r s (e.g ., 
iDunlop et all 120031 : iGuvon. Sanders, fe Stocktonl 1200a) . 
However, at z — 1.5 — 2.5, iKocevski et al.1 ( 20121 1 find a 
high f ractio n of d isk h osts for moderately luminous AGN, 
ISchawinski et af] (120121 ') find heavily obscu red quasar hosts 
to be mostly disks, and IPeng et al.l (j2006l l find signs that 
the hosts are of less concentrated morphology. On the other 
hand, it is plausible that, given the overall small effective 
radii, the even more compact and concentrated bulge com¬ 
ponent cannot be reliably distinguished and is integrated 
into the PSF by our modelling technique. As a result, we also 
performed fitting with de Vaucouleurs profile for the host 
galaxies. In particular, this needs to be used when check¬ 
ing if the same correlation between the luminosity of the 
host bulge and Mbh found in the literature holds in our 
AO sample. We also note that the S/N of all hosts is low, 
and visually the residuals after modelling with Sersic or de 
Vaucouleurs profiles look very similar. 

As the relations between the quasars and their host 
galaxies need to be determined in the quasar rest frame, 
we converted the observed magnitudes of the de-magnified 
hosts and quasars, recreated in the source plane, to absolute 
magnitudes in the rest-frame 77-band. We used the equation 


Mr = rriK — K(z) — 51ogdz, — 25 , (2) 

where K [z] is the Tv-co r recti on between the K s and R bands 
found bv lFalomo et al.l (120081 1. and cLl is the luminosity dis¬ 
tance (in Mpc). The choice of spectral template impacts the 
conversion at most by ~ 0.3 mag for SDSS J0832+0404. We 
did not apply an evolution correction. In addition, some ex¬ 
tinction may be caused by light being deflected close to the 
lensing galaxy. However, no differential extinction is found 
for any of the five systems in the foll ow-up spect r oscop y pre¬ 
sented in the discovery papers, and IPeng et al.l (|200d l con¬ 
clude that such extinction is likely negligible. We present 
the results in Tabic [7] 

We assumed that all five systems are radio quiet, as 
they are not found in any available catalogue of radio 
sources. Comparing the nuclear luminositi es to t h e samp le 
of hig h- and low-luminosity quasars from iKotilainen et al] 
(1200711 . all systems contain low-luminosity quasars, except 
for SDSS J1515+1511, which contains a high-luminosity 
quasar. 

In order to study the relation between the host galaxies 
and their central black h oles, Mb h needs to be estimated. 
An estimate exists in IShen et alj (1201 lh . where Mbh was 
computed for all SDSS quasars in Data Release 7. This es¬ 
timate however cannot be used, as it does not account for 
the lensing nature of the systems, the multiple images and 
the spectral contamination from the lensing galaxies. 

In order to estimate Mbh, as well as to facilitate 
comparison , we e mployed the same virial relation used by 
iPeng et akl (l2006l ): 


Mbh = An n e 


FWHM(line) 

2 

XL\ (Aline) 

1000 km s -1 


10 44 erg s' 1 


(3) 


where L = A L\ is the continuum luminosity at 
Aiine{CIV,MgII} = {1350,3000} A, Ai ine {CIV, Mgll} = 
14.5,6.1} x 10 6 , and 7 ii ne{CIV , Mgll} = {0.53,0.47} 

( Vestergaard fe Petersonl (120061 4 and iMcLure fe Jarvisl 
( 2002i) , respectively). At least one of the two emission lines, 
CIV and Mgll, is found in the original SDSS spectrum of 
each system. _ 

The measured line FWHMs in IShen et all (l201lh are 
not modified by lensing, and can therefore be employed. We 
determin ed the co ntinuum luminosities using the method of 
IPeng et al.l (l2006l f. First, by combining the discovery papers 
with the new AO data, all systems have individual resolved 
magnitudes for the quasar images and lensing galaxies in at 
least three bands. We converted the measured magnitudes 
of the bright quasar image in each system into flux densi¬ 
ties, which we then fitted with a power law F\ oc A“ 2 
and extrapolated to the observed wavelength Aiine (1 + z). 
We then divided the extrapolated flux densities by the mag¬ 
nification ratio of the bright quasar image, resulting in the 
de-lensed Fa, which we converted to the desired luminosi¬ 
ties L = 47rd|Aiine(l + s)Fa- Although this approach is an 
approximation, and does not account, for example, for un¬ 
reliable magnitudes in the low-resolution data or the ef¬ 
fect of strong emission lines on the magnitudes, Mbh is 
only weakly dependent on the estimated continuum lumi¬ 
nosity, roughly through its square root. We used the fitted 
power law to estimate the quasar nuclei Mr in Table [7] 
with K(z) = 2.5 log(a — 1), and the final values in the Vega 
system. 

We report the resulting Mbh, together with the 
Eddington ratios, in Table [7] SDSS ,10819+535 and 
SDSS J1322+1052, where the CIV line is affected by par¬ 
tial absorption, show differences in Mbh estimated from the 
Mgll and CIV lines, and therefore we discarded the result 
based on the CIV lines. For the four objects with reliable 
results, Mbh — 10 s — 10 9 M 0 . 

In Figure [5] we show the 5 objects in the curr e nt A O 
sample compared to the sample from IPeng et al.l (120061 ). 
They follow closely the correlation found between Mbh and 
the host galaxy luminosity. 

Figure [TD] plots the absolute 77-band magnitudes of the 
quasar hosts versus those of their nuclei, as well as the 17- 
band magnitudes as a function of redshift, fortlie AO sample 
as w ell as other qu a sars c ompiled bvl Kotilainen et al.l (l2007t) 
and iFalomo et akl |2008). In both plots, the four systems 
with robust host luminosities occupy the loci found for radio 
quiet quasars. 


8 LENS GALAXY LIGHT-MASS PROFILE 
AND ENVIRONMENT 

Although a dedicated study of the environment of the SQLS 
sample used in this paper is not possible, due to the lack 


8 The conversions were performed using the tool at 

http: //www. gemini. edu/?q=node/11119 
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Table 6. Quasar host galaxy properties 


Object 

mj^i host 

mj^/ nucleus 

e 

ReS H 

Reff [kpc] 

Sersic index 

host/qso 

SDSS J0819+5356 

20.09 ± 0.22 

19.26 + 0.11 

0.74+0.12 

0.26 + 0.07 

2.2 

0.57+0.66 

0.46 

SDSS J0819+5356 

19.49 ± 0.58 

19.13 + 0.13 

0.74+0.11 

0.33 + 0.30 

~ 2.8 

[4] 

~ 0.77 

SDSS J0832+0404 

19.21 ± 0.10 

17.95 + 0.08 

0.67+0.15 

0.38 + 0.10 

3.2 

0.25 + 0.32 

0.31 

SDSS J0832+0404 

18.46 ± 0.14 

18.22 + 0.12 

0.65 + 0.13 

0.39 + 0.19 

3.2 

[4] 

0.80 

SDSS J0904+1512 

19.39 + 0.05 

17.74 + 0.02 

0.25 + 0.16 

0.04 ± 0.00 

0.4 

2.11 ± 0.24 

0.22 

SDSS J0904+1512 

19.10 + 0.06 

17.73 + 0.06 

0.15 + 0.23 

0.04 + 0.01 

0.4 

[4] 

0.28 

SDSS J1322+1052 

£ 17.15 

~ 18.02 

~ 0.3 

- 

- 

- 

> 1.2 

SDSS J1515+1511 

19.18 + 0.10 

16.68 + 0.01 

0.30+0.05 

0.25 + 0.01 

2.1 

1.71 ± 0.80 

0.10 

SDSS J1515+1511 

18.85 ± 0.04 

16.72 + 0.01 

0.33 ± 0.04 

0.18 + 0.02 

1.5 

[4] 

0.14 


Error bars include those resulting from the simulations with a real PSF and noise, but not from the uncertainty in the 
radial slope of the lens mass profile. Here host/qso is the flux ratio in the source plane, and e is the ellipticity of the 
host. 


Table 7. Quasar host luminosity and black hole related properties 


Object 


Mr 

host 

L* 

a — 2 

Mr 

nucleus 

log Z/cont 
[log erg/s] 

FWHM 

[km/s] 


^bol 

-t'Edd 

SDSS J0819+5356 

2.237 

-23.7 

10.9 

-0.71 

-25.9 

44 . 6 , 44.7 

1114, 3096 

7.1, 8.1 

1.04, 0.16 

SDSS J0832+0404 

1.115 

-22.75 

4.6 

-1.75 

-23.95 

-, 45.1 

-, 5863 

-, 8.8 

-, 0.07 

SDSS J0904+1512 

1.826 

-23.4 

8.6 

-1.65 

-25.65 

45 . 7 , 45.5 

5754, 4364 

9.1, 8.8 

0.13, 0.22 

SDSS J1322+1052 

1.716 

> -25.2 

< 44.5 

-1.68 

-25.2 

45 . 6 , 45.4 

1723, 4553 

0° 

0 

0° 

bo 

1.26, 0.17 

SDSS J1515+1511 

2.054 

-24.1 

15.9 

-1.88 

-26.65 

46.4, 46.1 

3860, 4267 

9 . 1 , 9.0 

0.60, 0.47 


For all systems, the magnitudes and luminosities are calculated assuming Sersic index n = 4 for the host galaxies. Here 
a is the spectral power-law index (see text). The two values in columns 7 (8) refer to 1350 A(CIV) , and 3000 A(MgII), 
respectively; Corresponding values are given in columns 9 and 10. L* (Mb +=21.2; lGardner et al.l d 19971) ’) is the characteristic 
luminosity of the Schechter luminosity function for elliptical galaxies. L CO nt is the continuum luminosity in the proximity 
of the considered emission line. L^ 0 \/ Z/Edd is the Eddington ratio, where L^ 0 \ is the bolometric lu minosity computed from 
T 3000 and T 1350 with bolometric corrections 5.15 and 3.81, respectively, from the composite SED in Richards et aT] (120061 ’): 
Z/Edd — 47T GMBumpc/aT — 1.3 x 1O 38 Mbh/-^ 0 erg/s is the Eddington luminosity, in the usual units. 


of necessary spec troscopic data (e.g. jMomcheva et al.l l2006; 

I Wong et al.ll201ll ), in Appendix[B]we present considerations, 
whenever possible, based on the available AO data and the 
SDSS multi-band data. As mentioned in Section [lj a moti¬ 
vating goal is the high-precision determination of the Hub¬ 
ble constant. This is possible in the case of strongly lensed 
quasars with measured time delays, however biases are in¬ 
troduced by degeneracies between the lens galaxy mass pro¬ 
file and the environment. It is therefore useful to identify 
and prioritise the modelling of lenses with a comparatively 
better understood environmental effect when planning time- 
consuming monitoring, for the purpose of measuring time 
delays. In this section, the goal is to study based on our AO 
sample whether or not the light ellipticity can be used as a 
proxy for the mass ellipticity of the lens galaxies, and then 
to obtain an overall classification of the present AO sample 
in terms of the degree in which the environmental influence 
is understood. 

Using an appro a ch similar to th at _ of 

[Chantry. Sluse. fc Magainl (j2010l i and ISIuse et al.l (|2012l f . 
we modelled each object (except for the quad lens 
SDSS J1330+1810), whenever possible, with two indepen¬ 
dent lens models: SIS +7 and SIE+ 7 , plus the SIE model 
(generally equivalent with SIS+ 7 ), in order to estimate the 
effects of the lens ellipticity on the direction and strength 
of the predicted shear. Since for doubles the SIE +7 models 
have a negative number of degrees of freedom, we used 
a prior, which assumes that the observed luminous lens 


ellipticity and orientation also hold for the mass profile. The 
considerations in this section are based on the assumption 
that we measure accurately the ellipticity of the luminous 
profile of the lensing galaxies. The simulations we performed 
for each system, as described in Section ^. 61 were indeed de¬ 
signed to provide an accurate estimate of the uncertainties 
in these measurements, even when we approximate with an 
analytical or hybrid PSF the true, non-analytical PSF. In 
addition, many of our lenses have large ellipticity, which is 
less likely to be affected by inaccuracies in the PSF. 

8.1 Lens galaxy light-mass profile 

W e first pe rform an insightful comparison to the sample of 
1 Sluse etah llml). They model a sample of 10 doubles and 
14 quads from the CASTLES survey, further processed with 
image deconvolution. The only quad in the current AO sam¬ 
ple, SDSS J1330+1810, is likely a disk-like galaxy with a 
significantly larger observed ellipticity, ~ 0.7. In addition, 
the model ellipticity is highly dependent on the power-law 
model, with a value of ~ 0.35 for slope 7 ' = 2, and ~ 0.6 for 
slo pe 7 ' ~ 2 .4 (t he best-fit model assuming no perturbers); 
ISIuse et al.l (120121) assume a uniform slope 7 ' = 2 (isother¬ 
mal) . _ _ 

Among the 10 1 Sluse et al.1 (12012! ') doubles, 7 have shear 
orientations in agreement between the SIS +7 and SIE+ 7 . 
All of these have observed light cllipticiy e < 0.22 and in ad¬ 
dition, the SIS +7 model shear values are such that e < 37 
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M r (host) 


Figure 9. The correlation between the black hole mass and the 
absolute R— band magnitude of th e quasar hosts. The black points 
are from the sample of iPeng et al.l 11200a ). The line shows the best- 
fit rela tion found for local AGNs (not plotted here; see[Peng_et_a]J 
ll2006l ') for details). The error cross from IPeng et al.l 1 120061 ). 
shown here, is fairly representative for the AO sample (magenta 
points) as well, with the exception of SDSS J0819+5356. For 
SDSS J0819+5356 and SDSS J1322-t-1052, the Mbh estimates 
based on CIV are ignored. A colour version of this figure is avail¬ 
able in the online version of the paper. 


dKeeton. Kochanek. fc Seliaklll997 ). therefore the lens ellip- 
ticity is not the main source of shear. There are two excep¬ 
tions, in which case the light orientation and shear are nearly 
orthogonal, such that the intrinsic and extrinsic sources of 
shear are hard to disentangle. For the remaining objects, 
the SIS +7 and SIE +7 shear orientations are different, the 
observed ellipticity is larger, and the shear in the SIE +7 
model is nearly orthogonal to the lens. For these doubles (HE 
2149-2745, CLASS B1600+4 34 an d FBQS 0951+2635), the 
interpretation attempted in Isiuse et al. (2012 1 of shear in 
terms of environment is not secure, if the assumption that 
the light and mass ellipticity are the same. On the other 
hand, all three objects are consistent with the observed ori¬ 
entation and SIS +7 shear being in agreement to within ~ 10 
deg. 

In the subsequent analysis of our sample of confirmed 
doubles with robust modelling, we remove 9 systems which 
may bias our analysis. These systems have unreliable light 
cllipticity/orientation, or are too complicated to be reli¬ 
ably modelled with simple single SIE or SIE +7 models. 
They are SDSS J0904+1512 (ellipticity measurement not 
available), SDSS J1001+5027 (very close galaxy compan¬ 
ion), SDSS J1055+4628 (morphology uncertain due to pos¬ 
sible quasar host contamination), SDSS J1206+4332 (com¬ 
plicated system of lenses), SDSS J1216+3529 (very close 
galaxy companion), SDSS J1254+2235 (multiple Sersic corn¬ 
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Figure 10. Upper plot: The absolute R— band magnitudes of the 
quasar hosts, assuming Sersic index n = 4, and their nuclei, for 
the 5 objects with host galaxy estimates in the present AO sample 
(magenta) and various samples of radio quiet quasars at 1 < z < 
2. Lower plot: The absolute R— band magnitudes of quasar hosts 
and their redshift. The bl ack points are sampl e averages of radio 
quiet quasars, compiled in I Falomo et al.1 (|2008h . A colour version 
of this figure is available in the online version of the paper. 
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ponents), SDSS J1334+3315 (close galaxy companion re¬ 
sponsible for most of the shear), SDSS J1405+0959 (binary 
lens system) and SDSS J1406+6126 (multiple Sersic com¬ 
ponents). We keep however SDSS J1353+1138, which was 
modelled with two Sersic components of similar orientation 
and ellipticity. 

From the remaining sample of 13 doubles with ro¬ 
bust SIS +7 (SIE) and SIE +7 mass models, two objects, 
SDSS J0946+1835 and SDSS J1322+1052, have the shear 
orientation in the SIS +7 and SIE +7 models consistent with 
each other. These have observed light e < 0.10 consistent 
with e < 37 , resulting in the conclusion that the match be¬ 
tween shear model orientations is due to the fact that the 
internal shear created by the lens ellipticity is too small to 
account, and in fact significantly modify, the model shear. 
Qual itative ly, t hese d oubles are counterparts of the 7 objects 
in ISluse et all (l 2012 f) . 

Eleven objects in the current sample 
(SDSS J0743+2457, SDSS J0819+5356, SDSS J0820+0812, 
SDSS J0832+0404, SDSS J1131+1915, SDSS J1313+5151, 
SDSS J1353+1138, SDSS J1400+3134, SDSS ,11455+1447, 
SDSS ,11515+1511 and SDSS J1620+1203) have observed 
orientations and SIE +7 shear that are nearly orthogonal. 
This means that the 11/13 systems belong to this cate- 
gory are qu alitative counterparts of the three systems in 
ISluse et all J2OI2I ). Most of these have ellipticity > 0.25. 
Lens modelling shows that once the ellipticity in the SIE +7 
model is larger than e < 37 , where 7 is the shear strength 
of the SIS +7 model, the models require a nearly orthogonal 
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shear in order to reproduce the observed configuration. 
As this ellipti c ity ra nge lies outside the correlation range 
in ISluse et al.1 (l2012l f . the most likely explanation is that 
for these systems the overall mass ellipticity, dominated 
by the dark matter halo, is less elliptical than observed. 
Statistically, this is more likely than for such a large number 
of systems to be serendipitously located in environments 
where the main source of external shear acts orthogonal to 
the observed light orientation. 

There is nonetheless a theoretically predicted tendency 
for doubles (as opposed to quads), to have shears preferen¬ 
tially oriented perpendicular to the major axis of the lensing 
galaxies, but this is at a level of ~ 15% (Figure [lT] and ref¬ 
erence therein), certainly not enough to explain the large 
fraction of objects in the AO sample. Statistically, the ob¬ 
served and theoretical histograms in Figure [lT] have only 
a ~ 3% probability of representing the same distribution, 
following a Kolmogorov-Smirnov test. 

Therefore, this statistical analysis shows that in our 
sample of lensed quasars with an intermediate range of 
REin/Refc a few times larger than ~ 0.5 (with mean 
REin/Refi ~ 1.5, larger than the mean value for the SLACS 
and SL2S lenses), and at large observed light ellipticity 
> 0.25, the light profile tends to be less ellip tical than 
the m ass. A similar conclusion was suggested bv lDve et al.l 
ll20l4) from a much smaller sample of individual lenses with 
extended source modelling. 

In Figure [T^] we compile the mass versus light orien¬ 
tation and ellipticity for the robust doubles in our sam¬ 
ple, using the alternative SIE (n.b. not the SIE+ 7 ) model. 
The error bars on the mass ellipticity and orientation are 
boosted to include devia tions from isothermality consistent 
with the lens sample of iKoopmans et al.l (120091 1 . i.e. with 
slope 1.8 < 7 ' < 2.2. That is, the values of the mass el¬ 
lipticity, an d hence orientation , depen d on the radial mass 
profile fe.g.. iKeeton. Kochanek, fc Falcdll998l h and we add 
in quadrature the error bars in Table0for an isothermal pro¬ 
file ( 7 ' = 2 . 0 ) to the variation in ellipticity/orientation for 
the corresponding models with 7 ' = 1.8 and 7 ' = 2.2. The 
ellipticity plot shows that the great majority of lenses with 
observed ellipticity > 0.25 require mass models with smaller 
ellipticity than the light, even after we allow for freedom in 
the radial slope. This supports the conclusion above. The 
orientation (PA) plot shows scatter due to the influence of 
environm ent. The sc atter i s larger than for the sample of 
quads in ISluse et al.l (l 2012 f ), where the additional observa¬ 
tional constraints can more reliably separate external and 
internal sources of shear. 


8.2 Sample classification based on environment 

We have seen above that the use of the SIE +7 model with a 
prior on the observed light ellipticity/orientation is mislead¬ 
ing in breaking the degeneracy between the external shear 
and lens shape for our sample, as in most cases the light 
shape is not a good proxy for mass shape. Below we classify 
our sample based solely on the SIS +7 model and clues from 
the environment. We identify three types of objects: 

• For 6 objects (SDSS J1131+1915, SDSS ,11254+2235, 
SDSS ,11313+5151, SDSS J1353+1138, SDSS J1406+6126 
and SDSS J1515+1511), the measured light orientation of 
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Figure 11. Histogram of the difference in orientation between 
the observed lens and the shear (assuming the SIE +7 model) 
for the doubles in the AO sample. As noted in the text, 10 ob¬ 
jects which may bias our analyses were removed. Also plotted is 
the t heoretically p redicted alignment for doubles (solid line) from 
lOguri fc Marshall j201C)h ; the non-biased (non-lens) population is 
flat. 


the main lensing galaxy and the direction of the shear (in 
the SIS +7 model), or alternatively the orientation of the SIE 
mass model (without shear), match to within ~10 deg. As 
no prior on the ellipticity or orientation was used, the align¬ 
ment of the orientation predicted by the mass models and 
the observed light orientation is therefore real. As expected, 
the SIS +7 and SIE models are equivalent, the intrinsic lens 
ellipticity acting as internal shear, with e ~ 37 . For these ob¬ 
jects, the environmental influence is not necessary in order 
to model the system, and is likely unimportant. 

One caveat needs mentioning. SDSS J1254+2235 is mod¬ 
elled by a superposition of two concentric Sersic profiles, 
whose orientations differ by ~ 35 deg. The SIE mass model 
orientation corresponds to the orientation of the brighter 
profile. If the additional newly-discovered object G2 were 
close to the redshift of the lens, its effect on the model would 
be drastic. Also, for SDSS J1406+6126, which was modelled 
with a superposition of two Sersic profiles, the orientation of 
the extended profile was considered, the other being spher¬ 
ical. 

• For 12 objects (SDSS J0819+5356, SDSS J0820+0812, 
SDSS J0832+0404, SDSS J0946+1835, SDSS J1216+3529, 
SDSS J1322+1052, SDSS J1330+1810, SDSS ,11334+3315, 
SDSS ,11400+3134, SDSS J1405+0959, SDSS J1455+1447 
and SDSS J1620+1203), the environment is likely to play 
an important role, as the observed light orientation and 
required shear (the SIS +7 model) are different (see below 
for two exceptions), and there is some evidence as to what 
causes the external shear. 

For SDSS ,11322+1052, the difference in orientations 
is just ~17 deg, and there are several nearby galax¬ 
ies, at least one of which could affect the system. For 
SDSS J0819+5356, SDSS J0832+0404, SDSS J0946+1835 
and SDSS J1620+1203, the shears point toward nearby 
groups of galaxies, the likely source of influence. 
SDSS ,11330+1810, the only 4-image lens in the present sam¬ 
ple, has enough observational constraints to disentangle the 
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internal and external shear, the latter of which is due to 
one of two nearby galaxy groups. SDSS J0820+0812 is lo¬ 
cated behind an elongated galaxy group and a larger struc¬ 
ture, whose orientation corresponds to the required shear. 
For SDSS J1400+3134, the two orientations differ by ~ 14 
deg, but are brought in agreement if the effect of a nearby 
galaxy is considered. For SDSS J1334+3315, a bright nearby 
galaxy has been sho wn to be r espon sible for the major¬ 
ity of the shear (see Rusu et al.ll201lh . SDSS J1405+0959 
(see iRusu et all l2014f T is a complicated system with a bi¬ 
nary lens, and likely an additional contribution from the 
observed nearby galaxies located in the foreground. 

There are two exceptions, where the observed light orien¬ 
tation and required shear are in fact consistent at ~10 deg, 
but there is likely an environmental contribution. For both 
SDSS J1216+3529 and SDSS J1455+1447, a single nearby 
galaxy, which may affect each system, is also located close 
to the direction of the shear. For the latter system, explic¬ 
itly accounting for the companion galaxy in the lens model 
by scaling the strength of the two galaxies further improves 
the agreement between the shear orientation and the light 
orientation. 

• For 5 objects (SDSS J0743+2457, SDSS J0904+1512, 
SDSS J1001+5027, SDSS J1055+4628, and 
SDSS J1206+4332), the environment is likely to play 
an important role, but there is no clue as to what causes 
the measured external shear. For SDSS J0743+2457 and 
SDSS J1055+4628, the shear direction is significantly 
offset from the light orientation, and there is no obvious 
galaxy or galaxy group to account for the discrepancy. 
However the observed light orientation may be biased 
for SDSS J1055+4628 due to possible unaccounted 
host contribution during the morphological modelling. 
SDSS J0904+1512 has an ellipticity consistent with zero, 
therefore the source of the required shear is not known. 
SDSS J1206+4332 has several nearby galaxies, and it is 
unclear if they are physically associated, and what their 
exact contributions to the model are. SDSS J1001+5027 
requires a large external shear whose direction does not 
match the nearby foreground group of galaxies. 

The doubles that are most promising in terms of re¬ 
liably disentangling internal and external sources of shear, 
and would therefore be most useful for a determination of 
the Hubble constant based on time delay measurements, are 
SDSS J1131+1915, SDSS J1313+5151, SDSS J1334+3315, 
SDSS J1353+1138, SDSS J1400+3134, SDSS J1406+6126, 
SDSS J1455+1447 and SDSS J1515+1511. 


9 CONCLUSIONS AND FUTURE PROSPECTS 

In this paper we described the results of a large imaging 
campaign with the Subaru Telescope adaptive optics, aimed 
at obtaining high-resolution observations of gravitationally 
lensed quasars from SQLS. The purpose of the high reso¬ 
lution observations was multifold. First, it was to provide 
accurate astrometry and relative photometry in the near- 
infrared, which we used to obtain accurate lensing models 
and masses inside the Einstein radii. Second, it was to derive 
robust morphology for the lensing galaxies and compare it 
with the mass profile. Third, it was to study (if detected) 


the quasar host galaxies at high redshift ~ 2, as well as other 
potential features previously hidden in low-resolution data. 

In this work we introduced a new method for fitting the 
morphology of gravitationally lensed quasars in AO data. 
This was necessary, as the separately observed PSF stars 
proved to be ill-suited representations of the actual PSF, 
due to varying atmospheric turbulence characteristics. The 
new method consists of approximating the PSF analytically, 
and fitting it simultaneously to the quasar images and to 
the lensing galaxy (in the latter case, convolved with a Ser- 
sic profile). For the low Strehl ratios (< 0.10) typically ob¬ 
tained, we determined empirically that the best simple an¬ 
alytical PSF approximation is a superposition of 2 Moffat 
profiles. For the systems with larger signal-to-noise, and a 
relatively isolated brighter image, we created a hybrid PSF 
by using the observed brightness profile of the core, and re¬ 
placing the wings with analytical approximations, due to 
the high noise level. We used a variety of methods to obtain 
realistic error bars on the observed parameters, including 
simulations which employed the separately observed stars. 

In addition, using the natural stretch of the quasar host 
galaxy into arc-like features, we characterized the quasar 
host galaxy, or placed limits on its brightness, for 5 sys¬ 
tems. The hosts follow relations found in the literature with 
redshift, nuclear luminosity and black hole mass. This is 
the first time such study was attempted without an a-priori 
known PSF. While the small number of host detections in 
the present AO sample cannot at the present significantly 
improve the results of previous studies, by showing that 
these hosts follow a known correlation determined for HST 
data, and have similar error bars, we have shown that stud¬ 
ies of hosts with AO are possible to perform in the future 
on larger samples, in a way that is competitive and comple¬ 
mentary to HST -based studies. 

The typical achieved relative astrometry precision was 
~ 1 — 2 mas for the quasar images, but less for the lens¬ 
ing galaxies. This is comparable to the HST results. Based 
on the new astrometry, we have shown that for many sys¬ 
tems the original low-resolution astrometry was not accurate 
within the given error bars, in many cases providing unreli¬ 
able gravitational lensing mass models. Due to the necessity 
of modelling the PSF on the target, the morphological pa¬ 
rameters of the lensing galaxy were typically less accurate 
than those provided by HST. 

By using the SIE mass model together with the ob¬ 
served light distribution in the lensing galaxies, and using 
the known correlation between the orientation of light and 
mass, we have classified our sample according to the degree 
to which we expect the lens environment to influence the 
observed configuration. We have pointed out systems which 
seem well understood and would benefit most from time de¬ 
lay measurements in order to infer the Hubble constant. On 
the other hand, we have shown that the tight correlation 
found for the elliptical galaxy lenses in the SLACS Survey 
between the mass and light ellipticity does not hold statis¬ 
tically for the present sample at ellipticity e > 0.25, and the 
mass profile is less elliptical than the light profile. This is 
likely due to the larger mean ratio of the Einstein radius to 
the effective radius encountered in the present sample, the 
mass profile therefore being more dark-matter dominated 
and rounder. 

As a result of the high-resolution observations, we ob- 
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Light ellipticity 


Figure 12. Light and mass orientation (left), as well as ellipticity (right) for 14 objects without obvious sources of bias or systematics 
in the AO sample. For the observed (light) profiles, the values and error bars are taken from Table [I] However for the SIE (n.b. not the 
SIE+ 7 ) mass models, the error bars from Table[5]are boosted to account for variations in the power-law slope, by considering variations 
in mass orientation and ellipticity for 1.8 < 7 ' < 2.2. For the quad SDSS J1330+1810, only the power law +7 + <5 model is considered, 
without boosting the errors. For SDSS J1353+H38, the light profile of the flatter Sersic component is considered. 


tained the following noteworthy results on particular sys¬ 
tems (see Appendix iBl for details): 

• SDSS J0819+5356 and SDSS J0832+0404: we charac¬ 
terised the host galaxy with the use of an analytical PSF. 
SDSS J0904+1512 and SDSS J1322+1052: we used an iter¬ 
ative approach with a hybrid PSF to model the quasar host 
galaxy. SDSS J1515+1511: this is the only system where 
a star serendipitously located in the field has provided a 
suitable PSF, resulting in the detection of the quasar host 
galaxy. 

• SDSS J1001+5027 and SDSS J1206+4332: the sys¬ 
tems have available time delay measurements and con¬ 
sist of multiple lensing galaxies. The secondary galaxy in 
SDSS J1001+5027 is either located at different redshift, or 
is undergoing tidal stripping. The system also suggests the 
existence of a massive nearby galaxy cluster, which future 
field spectroscopy would help confirm. In SDSS J1206+4332, 
we discovered a previously undetected galaxy in the non-AO 
data. In order to reliably model each system however, spec¬ 
troscopy is required to estimate the redshift of the individual 
lensing galaxies. 

• SDSS J1254+2235, SDSS J1353+1138 and 
SDSS J1406+6126: the two galaxies are not fitted well by 
a single Sersic profile, but instead require a combination 
of two Sersic profiles each, revealing more morphological 
detail than in the previously available low-resolution data. 

• SDSS J1330+1810: in opposition to the discovery pa¬ 
per, we showed that the lensing galaxy is a late-type galaxy. 
Also, we discovered astrometric anomalies when applying 
the customary SIE +7 mass model. Although substructure 
in the form of lens galaxy satellites is not detected but can¬ 
not be ruled out either, a power law ( 7 ' ~ 2.3) model with a 


contribution from the complicated environment in the form 
of a third order perturbation fits the data well. 

• SDSS J1 405+0959: we p resented this system in a sepa¬ 
rate paper dRusu et al.l 120141 ) : this is an interesting system 
where the AO observations have revealed dramatical new de¬ 
tail, with a host galaxy arc produced by two lensing galaxies. 

Through the observations and subsequent data anal¬ 
ysis presented here we have shown that a large campaign 
to observe gravitationally lensed quasars with AO is fea¬ 
sible even with the lack of an a-priory known PSF. We 
used the data to obtain new results on individual systems, 
as well statistically, in combination with other samples of 
objects. Upcoming wide-held imaging surveys such as the 
Large Synodic Sky Survey and Subaru Hyper Suprime-Cam 
Survey will discover thou sands of additional lensed quasars 
llOguri &; Marshalill2010l) . and LGS/AO is a viable alterna¬ 
tive to HST for the follow-up of some of these new systems 
with high resolution observations. 

The technique employed can be further enhanced by the 
availability of a more accurate distortion correction map, 
and may be further enhanced by deconvolution techniques. 
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APPENDIX A: TABLES OF TIP-TILT STARS, 
PSF STARS, ANALYTICAL PSF PARAMETERS, 
AND ASTROMETRIC DISTORTION 

APPENDIX B: DESCRIPTION OF RESULTS 
ON INDIVIDUAL OBJECTS 

A discussion is given below for the objects observed dur¬ 
ing the present campaign, summarising the available data, 
analysis method used, as well as conclusions revealed from 
the lens mass models. The relative astrometry and photom¬ 
etry of each system, the morphological parameters of the 
galaxies, as well as the gravitational lens mass models are 
compiled in Tables 0 5] and [5] respectively. The objects for 
which the interpretation of the data is problematic are dis¬ 
cussed in Appendix 0 
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Table Al. Tip-tilt and PSF stars 


Object 

Description 

Magnitude 

Separation [ /; ] 

SDSS J074355.71+245745.62 

TT J0743+2457 

R=15.00 

~42.0 

SDSS J082006.11+535620.42 

TT J0819+5356 

R=11.18 

~55.5 

SDSS J083219.00+040505.03 
SDSS J083427.71+031632.82 
SDSS J083425.13+031533.49 

TT J0832+0404 

TT (PSF) J0832+0404 
PSF J0832+0404 

R=13.84 

R=13.9 

K=12.4 

~66.8 

71.2 

SDSS J090406.54+151229.94 
SDSS J090409.04+151214.85 

TT J0904+1512 

PSF J0904+1512 

R=16.36 

K=13.87 

~42.4; 39.2 

SDSS J092635.53+310037.66 

TT J0926+3059 

R=12.4 

~ 53.2 

SDSS J094605.79+183553.34 
SDSS J094333.27+182720.20 
SDSS J094332.33+182735.79 

TT J0946+1835 

TT (PSF) J0946+1835 
PSF J0946+1835 

R=14.43 

R=14.60 

K=12.36 

~20 

20.6 

SDSS J 100135.04+502740.84 
SDSS J100216.17+501059.07 
SDSS J100222.74+501117.49 

TT J1001+5027 

TT (PSF) J1001+5027 
PSF J1001+5027 

R=18.02 

R=18.04 

K=12.24 

~63.3 

65.7 

SDSS J 100234.48+445036.47 
SDSS J 100229.61+445002.77 

TT J1002+4449 

PSF J1002+4449 

R=13.47 

K=15.16 

~75.6 

61.8 

SDSS J105440.61+273222.93 

TT J1054+2733 

R=16.89 

~43.4 

SDSS J105538.55+462804.93 
SDSS J105208.39+463816.93 
SDSS J105200.79+463830.8 

TT J1055+4628 

TT (PSF) J1055+4628 
PSF J1055+4628 

R=13.7 

R=13.9 

K=11.85 

~79.2 

79.3 

SDSS J112821.12 +240253.1 

TT J1128+2402 

R=16.84 

~50.9 

SDSS J113201.76+191518.63 
SDSS J112825.17+194704.17 
SDSS J112827.97+194747.46 

TT J1131+1915 

TT (PSF) J1131+1915 
PSF J1131+1915 

R=15.33 

R=14.73 

K=10.90 

~58.3 

58.6 

SDSS J121640+352949.51 
SDSS J121505.47+330437.83 
SDSS J121507.13+330544.54 

TT J1216+3529 

TT (PSF) J1216+3529 
PSF J1216+3529 

R=15.6 

R=15.5 

K=11.7 

~72.9 

69.5 

SDSS J 125414.9+223530.99 
SDSS J125144.95+222445.79 
SDSS J125141.47+222521.42 

TT J1254+2235 

TT (PSF) J1254+2235 
PSF J1254+2235 

R=17.59 

R=17.37 

K=10.48 

~56.0 

59.7 

SDSS J131343.88+515214.61 
SDSS J131618.12+512544.78 
SDSS J131619.66+512449.66 

TT J1313+5151 

TT (PSF) J1313+5151 
PSF J1313+5151 

R=16.32 

R=16.16 

K=11.51 

~58.8 

57.1 

SDSS J132058.29+164459.07 

TT J1320+1644 

R=16.7 

~58.0 

SDSS J132238.43+105237.71 

TT J1322+1052 

R=16.53 

~29.7 

SDSS J133023.35 +180957.3 

TT J1330+1810 

R=14.77 

~75.5 

SDSS J133402.02+331437.75 

TT J1334+3315 

R=13.64 

~57.0 

SDSS J135303.52+113752.67 

TT J1353+1138 

R=15.05 

~42.8 

SDSS J140017.13+313511.05 

TT J1400+3134 

R=14.1 

~58.8 

SDSS J140514.96+095820.10 
SDSS J140517.09+103007.16 
SDSS J140521.31+102933.31 

TT J1405+0958 

TT (PSF) J1405+0958 
PSF J1405+0958 

R=13.35 

R=13.1 

K=12.0 

~71.5 

70.8 

SDSS J140626.66+612556.94 

TT J1406+6126 

R=14.37 

~46.0 

SDSS J145459.41+144838.21 

TT J1455+1447 

R=16.54 

~73.0 

SDSS J151538.33+151232.87 
SDSS J151537.14+151114.41 

TT J1515+1511 

PSF J1515+1511 

R=13.37 

K=14.79 

~57.2 

80.3 

SDSS J162022.49+120343.81 

TT J1620+1203 

R=13.98 

~52.0 


The magni tudes are taken from the the Naval Observatory Merged Astrometric Dataset 
fNOMAD; IZacharias et al]|2005h catalogue. “Separation” specifies the distancePetween 

the tip-tilt star and the object, and between the tip-tilt and the PSF star, respectively. ’ — 
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Table A2. Analytical parameters of the best-fitted PSFs 


Object 

FWHMl 

el 

PA1 

n i 

FWHM2 

e2 

PA2 

/32 

fluxl/ 

(fluxl + flux2) 

Strehl 
ratio (%) 

SDSS J0743+2457 

0.29 

0.30 

76.5 

3.1 

0.71 

0.10 

75.0 

2.0 

0.34 

5 

SDSS J0743+2457 star 

0.21 

0.20 

-88.5 

2.8 

0.69 

0.04 

69.3 

2.3 

0.44 


SDSS J0819+5356 

0.25 

0.17 

-70.8 

1.5 

0.92 

0.09 

61.0 

1.9 

0.36 

2 

SDSS J0819+5356 star 

0.23 

0.25 

19.8 

2.6 

0.80 

0.04 

18.6 

1.8 

0.21 


SDSS J0832+0404 

0.24 

0.16 

-70.3 

1.4 

1.20 

0.23 

4.2 

2.8 

0.38 


SDSS J0832+0404 PSF 

0.19 

0.04 

47.0 

9.0 

0.52 

0.12 

-49.2 

2.2 

0.22 

5 

SDSS J0904+1512 

0.21 

0.31 

81.6 

4.5 

0.77 

0.04 

35.5 

2.8 

0.32 

5 

SDSS J0904+1512 PSF 

0.21 

0.16 

33.2 

3.3 

0.72 

0.11 

-63.6 

3.2 

0.32 

4 

SDSS J0926+3059 

0.21 

0.29 

78.8 

21.1 

0.62 

0.10 

-81 

1.8 

0.09 

2 

SDSS J0946+1835 (2011) 

0.26 

0.23 

78.4 

3.8 

0.81 

0.12 

-47.7 

2.9 

0.57 

5 

SDSS J0946+1835 (2012) 

0.16 

0.14 

-86.1 

2.7 

0.98 

0.02 

82.1 

5.5 

0.55 

9 

SDSS J0946+1835 PSF 

0.14 

0.07 

-35.9 

3.0 

0.90 

0.04 

-28.4 

4.5 

0.57 

11 

SDSS J1001+5027 

0.20 

0.11 

-73.9 

5.3 

0.67 

0.02 

31.1 

2.0 

0.22 

3 

SDSS J1001+5027 PSF 

0.28 

0.18 

66.8 

7.3 

0.63 

0.04 

79.8 

2.3 

0.26 

2 

SDSS J1002+4449 

0.21 

0.22 

29.4 

12.5 

0.48 

0.27 

19.9 

2.6 

0.20 

3 

SDSS J1002+4449 PSF 

0.21 

0.11 

-11.3 

4.4 

0.58 

0.14 

2.9 

2.5 

0.30 

5 

SDSS J1054+2733 

0.34 

0.13 

69.9 

8.6 

0.83 

0.05 

53.4 

2.3 

0.09 


SDSS J1055+4628 

0.17 

0.12 

74.4 

3.8 

0.84 

0.25 

-89.0 

3.2 

0.26 

5 

SDSS J1131+1915 

0.15 

0.09 

-84.2 

3.7 

0.65 

0.05 

-28.1 

2.4 

0.39 

12 

SDSS J1131+1915 PSF 

0.14 

0.08 

-88.5 

4.0 

0.58 

0.22 

5.2 

2.4 

0.39 

12 

SDSS J1216+3529 

0.46 

0.20 

21.3 

[100] 

0.83 

0.13 

-56.0 

2.3 

0.48 


SDSS J1254+2235 

0.19 

0.20 

88.2 

8.1 

0.56 

0.15 

-52.7 

2.0 

0.31 

6 

SDSS J1254+2235 PSF 

0.21 

0.09 

-86.6 

18.2 

0.54 

0.09 

-85.3 

2.1 

0.41 

6 

SDSS J1313+5151 

0.18 

0.15 

72.4 

3.9 

0.63 

0.12 

-79.1 

2.0 

0.25 

5 

SDSS J1313+5151 PSF 

0.20 

0.11 

-84.5 

2.1 

0.72 

0.05 

32.0 

2.3 

0.50 

6 

SDSS J1322+1052 

0.22 

0.24 

60.1 

7.0 

0.61 

0.03 

48.0 

2.3 

0.35 

5 

SDSS J1330+1810 (2012) 

0.16 

0.06 

-52.8 

5.6 

0.61 

0.10 

-56.4 

2.1 

0.35 

8 

SDSS J1330+1810 (2013) 

0.18 

0.37 

72.9 

3.9 

0.68 

0.02 

-83.9 

1.7 

0.30 


SDSS J1334+3315 J 

0.14 

0.08 

-48.1 

15.3 

0.32 

0.07 

-67.3 

2.1 

0.29 

8 

SDSS J1334+3315 H 

0.14 

0.10 

-34.5 

4.7 

0.50 

0.08 

-66.5 

3.2 

0.51 

11 

SDSS J1334+3315 K' 

0.15 

0.07 

-35.8 

7.1 

0.42 

0.19 

-81.2 

2.5 

0.42 

12 

SDSS J1353+1138 

0.21 

0.20 

62.7 

7.5 

0.61 

0.08 

-30.8 

2.4 

0.37 

5 

SDSS J1400+3134 

0.24 

0.28 

34.5 

6.6 

0.66 

0.07 

78.4 

3.0 

0.35 

4 

SDSS J1405+0959 J 

0.18 

0.07 

21.7 

6.0 

0.55 

0.16 

-75.9 

2.5 

0.15 

2 

SDSS J1405+0959 H 

0.17 

0.07 

-48.0 

8.3 

0.51 

0.27 

-63.2 

2.2 

0.16 

4 

SDSS J1405+0959 H PSF 

0.17 

0.09 

12.8 

43.3 

0.46 

0.14 

-42.5 

2.3 

0.27 

6 

SDSS J1405+0959 K' 

0.18 

0.11 

27.9 

6.3 

0.54 

0.07 

-60.4 

1.7 

0.21 

7 

SDSS J1406+6126 

0.69 

0.12 

-70.3 

2.0 

- 

- 

- 

- 

- 

- 

SDSS J1455+1477 

0.19 

0.26 

-63.0 

2.4 

0.67 

0.07 

-52.0 

2.1 

0.22 

3 

SDSS J1515+1511 

0.21 

0.08 

-40.2 

8.2 

0.64 

0.03 

-46.7 

2.0 

0.17 

3 

SDSS J1515+1511 PSF 

0.21 

0.07 

-40.0 

10.3 

0.63 

0.04 

-45.4 

2.0 

0.19 

6 

SDSS J1620+1203 (2011) 

0.16 

0.16 

-42.3 

3.0 

0.55 

0.07 

41.4 

2.0 

0.50 

9 

SDSS J1620+1203 (2013) 

0.17 

0.13 

50.4 

4.0 

0.50 

0.01 

-63.8 

1.8 

0.33 

7 

Analytical parameters for all best-fitted object 

and star PSFs, 

determined 

with 

Host lens. 

The fitting was performed on the 


coadded, final frames. Affix 1 refers to the core, and affix 2 refers to the wings, both modelled as Moffat profiles. Here FWHM 
("), e, PA (deg) and /3 are the profile full width at half maximum, ellipticity, position angle (positive from North towards East), 
and (3 parameter, respectively. Values inside square brackets are kept fixed. The Strehl ratios were computed by comparing the 
peak flux of the PSF with that of a diffraction limited PSF model. 
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Table A3. Residual distortion error 


Object (dith. pos.) 

X CT pa i r sep. 

Y <7pair sep. 

^atm. refract. 

X 0"pix. scale 

Y <3"pix. scale 

X Ototal 

Y <^total 

SDSS J0743+2457 

[1.5] 

[1.5] 

<1 

0.2 

0.2 

1.1 

1.1 

SDSS J0819+5356 

[1.5] 

[1.5] 

~1.5 

0.3 

0.3 

1.5 

1.5 

SDSS J0832+0404 (5/6) 

5.5 

5.5 

<1.0 

0.1 

0.1 

3.9 

3.9 

SDSS J0904+1512 (5/5) 

1.6 

0.7 

< 1 

0.1 

0.2 

1.1 

0.5 

SDSS J0946+1835 (52mas) 

[1.5] 

[1.5] 

~1.0 

0.3 

0.2 

1.3 

1.3 

SDSS J0946+1835 (20mas) (9/9) 

1.7 

1.4 

~1.0 

0.4 

0.4 

1.4 

1.2 

SDSS J1001+5027 (7/9) 

0.8 

0.6 

~1.0 

0.2 

0.3 

0.9 

0.9 

SDSS J1055+4628 

[1.5] 

[1.5] 

< 1 

0.2 

0 

1.1 

1.1 

SDSS J1131+1915 (7/7) 

0.9 

2.1 

< 1 

0.1 

0.2 

0.6 

1.5 

SDSS J1206+4332 

[1.5] 

[1.5] 

< 1 

0.3 

0 

1.1 

1.1 

SDSS J1216+3529 

[1.5] 

[1.5] 

< 1 

0.2 

0.1 

1.1 

1.1 

SDSS J1254+2235 (9/9) 

1.6 

2.5 

< 1 

0.1 

0.2 

1.1 

1.8 

SDSS J1313+5151 (9/9) 

1.0 

0.7 

< 1 

0.2 

0.1 

0.7 

0.5 

SDSS J1322+1052 (5/5) 

2.1 

1.7 

< 1 

0.2 

0.2 

1.5 

1.2 

SDSS J1330+1810 (52mas) 

[1.5] 

[1.5] 

< 1 

< 0.2 

< 0.2 

1.1 

1.1 

SDSS J1330+1810 (20mas) 

[1.5] 

[1.5] 

< 1 

< 0.3 

< 0.4 

1.2 

1.0 

SDSS J1334+3315 (4/5) (J) 

2.2 

1.7 

<1 

0.1 

0.1 

1.6 

1.2 

SDSS J1353+1138 (5/5) 

1.2 

1.6 

<1 

0.2 

0.1 

0.9 

1.1 

SDSS J1400+3134 (4/5) 

1.6 

2.0 

~1.5 

0.2 

0.2 

1.6 

1.8 

SDSS J1405+0959 (J, H, K') 

[1.5] 

[1.5] 

< 1 

0.2 

0.1 

1.1 

1.1 

SDSS J1455+1447 (8/8) 

2.2 

2.1 

<1.0 

< 0.1 

0.1 

1.5 

1.5 

SDSS J1515+1511 (8/9) 

1.7 

1.0 

< 1 

0.2 

0.2 

1.2 

0.7 

SDSS 11620+1203 (2013) (4/5) 

3.8 

3.6 

~1.0 

0.2 

0.3 

2.8 

2.6 


The detector X axis is typically aligned to the north-south direction, and therefore corresponds to the Y axis in Table l3l In the 
Objects column, the values in parenthesis refer to the number of dither positions used to calculate cr pair sep Here cr pa i r sep> is 
the standard deviation of the quasar image separation between different dither positions. cr atm . re fract. is the expected stretch 
towards zenith due to the atmospheric distortion. cr p j x sca i e is the uncertainty in the image separation due to the uncertainty 
in the pixel scale (\/No.pixels x average pixel size error). cr tota i is calculated for the coordinates of an individual image on each 
axis as y'CT^ air sep + cr^ tm refract + scale /V^- Values inside square brackets are for objects for which <7 pair sep . could not be 
calculated, and the average value from the objects with available calculations was used. 


B1 SDSS J0743+2457 

SDSS J0743+2457 is a small-separation (6 ~ 1") dou- 
ble discovered inde pendently in the MUSCLES s urvey 
lljackson et al.ll2012h and in SQLS lllnada et al.ll2014h . The 
source redshift is z s = 2.165, an d the spectroscopic lens 
redshift is z s — 0.381 lllnada et alJl20l3 ). Follow-up low- 
resolu tion o bserva tions were obtained in the VRI bands 
lllnada et al.ll 20 lT identifying the lens galaxy. 

The AO observations were performed in NGS mode, 
and there is a bright star in the FOV, closer to the NGS in 
relation to the system (Figure IB II) . This star was used as a 
PSF to perform simulations of the system. The system itself 
was modelled using an analytical PSF, which resulted in 
clean residuals. Image B is particularly close to the lensing 
galaxy, but the two can nonetheless be disentangled. 

As the orientation of shear in the SIS +7 model (Table 
EJ is significantly different from the light orientation of the 
lensing galaxy (Tablc[4]), and there are no nearby perturbing 
galaxies in the FOV, the origin of the shear, although rela¬ 
tively small, cannot be determined. The shear orientation is 
vastly different in the discovery paper (by ~ 60 deg), where 
no error bars are provided. The SIE +7 model requires an 
unreasonably large amount of shear. We rule out that an 
inaccurate estimate of the galaxy morphology could bring 
the galaxy and shear orientation to match, as the ~ 55 deg 
difference is very large, and the very large lens ellipticity 
would also facilitate the measurement of its orientation. 


B2 SDSS J0819+5356 

The_ discove ry o f SDSS J0819+5356 is reported in 
llnada et al.l (120091 '). where it was followed-up at low reso¬ 
lution in the VRI bands. It is a large-separation (~ 4") 
galaxy-scale lensed quasar at z B = 2.237, with spectroscopi¬ 
cally _estimated__Zi = 0.294. Even in the low-resolution data, 
llnada et al.l (120091 ) notice an extended component C near A, 
and surmise it may be the quasar host galaxy. 

Due to non-photometric observing conditions, the zero- 
point was calculated in the new NGS AO data using the stars 
in the FOV. Two of the stars have catalogue magnitudes 
with which the observed relative photometry is in agree¬ 
ment. The resulting photometric accuracy is dominated by 
the catalogue error bar of 0.1 mag. The system was mod¬ 
elled with an entirely analytical PSF, including the quasar 
host galaxy, which resulted in clean residuals. Simulations 
were performed using the star towards south in Figure [B2l 
which is located at about the same distance from the NGS. 
The large error bars on the magnitudes of A and B (Table 
H reflect the uncertainty in the host-quasar decomposition. 
Systematics due to large Poisson noise were not found when 
the system was remodelled after masking the central 4 pixels 
of image A. The inferred quasar host properties are shown in 
Table [ 6 ] where the fit was done for a host galaxy of uncon¬ 
strained Sersic index, as well as n = 4. Residuals are shown 
in Figure IB2I As a note, if the system is modelled without 
a host, the morphology of the lensing galaxy changes dra¬ 
matically (effective radius ~ 1.2”, Sersic index ~ 5). The 
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Figure Bl. Modelling of SDSS J0743+2457. North is up and East is to the left, left: original frame (5" X 5”), left center: lens 
galaxy revealed after the subtraction of the quasar images. Right center: residuals after the subtraction of all objects. Right: FOV for 
SDSS J0743+2457 (66" x 66"). 


companion galaxy G2 (Figure TB2I) is disk-like, according to 
the SDSS classification, and requires two Sersic components 
to model. 

The morphological parameters of the len s galaxy G1 
are ge nerally similar to the ones reported in llnada et alJ 
(120091) . with the light orientation differing by 5 deg and the 
Sersic index slightly smaller. However the effective radius 
is half the value in the discovery paper. A possible reason 
could be that for this particularly extended object the sky 
background estimate is less reliable in the AO data. The 
wings of a steep Sersic profile are known to be particularly 
sensitive to the sky level. 

The shear direction in the S 1 S +7 model is offset by ~ 
25 deg from the light orientation, suggesting external influ¬ 
ence. The model parameters are very similar to the ones in 
the discovery paper, due to the large image separation which 
allows more robust astrometry. If accounting for G2 in the 
lens modelling, the shear direction changes by ~ 16 deg, 
making it more aligned with the observed light orientation. 
Although the SDSS photometric redshift of G2 is consis¬ 
tent the redshift of Gl, it is a disk-like galaxy, therefore its 
effect is expected to be smaller than inferred from the Faber- 
Jackson law. The value of the shear is not very large, partic¬ 
ularly in the SIS +7 model, and Gl is the brightest among 
the nearby galaxies, suggesting that it may be located close 
to the centre of a group. There is a group of galaxies (Figure 
IB3I) consistent with the redshift of Gl, surrounding the sys¬ 
tem but concentrated towards west. This is consistent with 
the rough direction of the shear in the mass models. In ad¬ 
dition, the large image separation/velocity dispersion of the 
lens may be boosted by the environment. 

B3 SDSS J0820+0812 

SDSS J0820+0812 was observed without AO. The object 
was originally reported in ljackson, Ofek. fe QgurU (I2009T) , as 
a two-image lensed quasar with separation ~ 2.3", with z s = 
2.024 and zi = 0.803. The observations in that paper were so 
unreliable that no astrometry/lens modelling is published. 

In our new data, with the exception of the target, there 
are two galaxies and a star in the FOV. Due to the pres¬ 
ence of the star in the FOV, the morphology of the system 


was fitted with Galfit in the traditional way. Because the 
original PSF FWHM (~ 30 pixels = 1.35") is extremely 
oversampled, 2x and 4x pixel binning was performed. The 
photometry of all components and the effective radius of the 
galaxy are the same, but the positions, Sersic index, ellip- 
ticity and position angle (PA) differ for different binning. 
The results from 4x pixel binning (final pixel scale 0.210”, 
Figure lB4ll were used, because the PSF for 2x pixels bin¬ 
ning is still largely oversampled, the S/N is lower, and the 
best-fit Sersic index ~ 0.25 is unreasonably small. The sky 
was estimated in four different empty regions surrounding 
the target, the dispersion sigma of the values was calculated, 
and the sky was held fix at ±lcr when using Galfit. For each 
parameter, the statistical error reported by Galfit and the 
dispersion due to the sky was added in quadrature. The lat¬ 
ter is dominant in most cases, including for the astrometry. 

The 2MASS catalogue was used to calibrate the pho¬ 
tometry. The K— band magnitude of the star is 15.63±0.21, 
and 15.53 ± 0.14 for the galaxy. Since the error bar on the 
galaxy is smaller, and there is also the possibility that the 
star is variable, the galaxy aperture photometry was used for 
calibration. The magnitude of the star then becomes 15.87, 
marginally inconsistent with the 2MASS value. 

A small Sersic value ~ 1 is preferred for the lensing 
galaxy. However, if the Sersic index is fixed to 4, the dif¬ 
ference is statistically small (Ax 2 /d.o.f. = 0 . 002 ), and the 
residuals look the same. The flux ratio 2.03 ± 0.13 is sig¬ 
nificantly lower than reported in the discovery paper (~ 6 ). 
Even if the Sersic index 4 is used, it is no larger than ~ 3. 
This is likely due to the poor quality of the data used in 
the discovery article. In addition, the flux ratio determined 
here is obtained in a redder band, and therefore less prone 
to quasar intrinsic variability, microlensing and extinction 
effects. 

The large lens velocity dispersion ~ 270 km/s may be 
boosted by the lens environment, and indeed there are seven 
nearby bright galaxies in SDSS photometric redshift consis¬ 
tent with ~ 0.19 (therefore in the foreground of the lens), 
surrounding the system. Most of the group members are 
located on the NE-SW direction. In addition, many more 
galaxies consistent with this redshift are found, particularly 
towards SW (Figure fB3l) . This is consistent with the direc¬ 
tion of the shear in both SIS +7 and SIE +7 models. Since 
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Figure B2. Hostlens modelling for SDSS J0819+5356. North is up and East is to the left. Top left: original frame (10" X 10"), top 
centre: residuals after modelling with an analytical PSF, without considering a host. Top right: residuals after modelling with host 
galaxy, which is not subtracted. Bottom left: residuals after modelling with a host galaxy and subtracting all components. Bottom right: 
original 66" X 66 " FOV. 


the differences in shear position angle are small in the two 
models (however the error bars are very large), this suggests 
that the environmental contribution is the main source of 
shear, dominating over the lens internal cllipticity and ex¬ 
act mass profile. 

B4 SDSS J0832+0404 

SDSS J0832+0404 was reported in lOguri et al.1 (i2008al ) as 
a double with z a = 1.115 and zi = 0.659, with available 
V, I, K band images. 

In the AO data, the presence of a faint arc-like fea¬ 
ture towards East (Figure IB5I) . which is significantly re¬ 
duced if the quasar host galaxy is considered, indicates a 
host detection. As few residuals remain if the PSF is mod¬ 
elled analytically, a hybrid PSF was not considered, and 
therefore the system was modelled in a manner similar to 
SDSS J0819+5356. In addition, since the overall AO cor¬ 
rection is similar to that of the bright star in the FOV of 
SDSS J0819+5356 (Table lA2l) . that star was used to perform 
simulations. The parameters of the host galaxy are given in 
Table |6] 

This object is also available in archival HST obser¬ 
vations (see Section E2I). The 3 x 400 s observations 
were obtained with t he Wide-Field Planetary Camera 2 
llHoltzman et al.lll995l) on 2009 July 8. We performed stan¬ 
dard reduction using Multidrizzl^fl for distortion correction, 

9 Multidrizzle, part of the STSDAS package, is a product of the 


cosmic ray removal, frame alignment and rotation. The re¬ 
duced science frame is shown in Figure IB6I We used two 
stars in the FOV as PSF templates (Figure fB6l) . and mod¬ 
elled the system using Galfit. Compact residuals remain af¬ 
ter the subtraction of all objects, particularly at the location 
of the lens and of image A. These may be due either to vari¬ 
ations in the PSF across the field, or due to Poisson noise. 
We checked that compact residuals remain if, for example, 
the star PSF1 is modelled with the PSF2 template, in sup¬ 
port of the first hypothesis. The fact that the residuals are 
compact and are not visible in image B suggest that they 
are not due to the quasar host galaxy, which appears to be 
undetected. This is not surprising even though the HST ob¬ 
servations are deeper, since the F606W filter falls shortward 
of the 4000 A rest-frame break of the host (assuming it is 
an early-type galaxy). 

In addition to A, B and Gl, an extended object GX is 
visible very close to image A, stretched tangentially to the 
location of the lens. This object, if indeed detected at all in 
the AO data, is much fainter, therefore has a blue colour. Its 
elongated morphology and position close to A suggest it is 
a lensed background object, although it cannot be ruled out 
that it is associated with the lens, without additional colour 
information. In the SIS+7 model, the shear direction does 
not match the position of this object, making it unlikely to 
be associated with the lens. If the redshift of GX it smaller 
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Figure B3. SDSS environment (left) and redshift cuts (right) for selected s ystems. A H g alaxie s kept after the cut have i < 21 and 
photometric redshift error A z <0.2, ensuring accurate photometric redshifts IlCsabai et al.|[2003h . The cuts satisfy at least one of the 
two redshift estimates in the SDSS database, PHOTOZ1 and PHOTOZ2 (see http://www.sdss3.org/drl0/algorithms/photo-z.php). 



Figure B4. Galfit/Hostlens modelling for SDSS J0820+0812. Left: original frame (12.6" X 12.6", 4x4 binning); centre: lens galaxy 
revealed after the subtraction of the quasar images. Right: residuals after the subtraction of all objects. 


than ~ 1.2, it is singly imaged, whereas if not, a counter 
image should be produced, given its proximity to the lens, 
of flux ~ 1/3. Such image should likely be observable in the 
HST data in the proximity of B, but is undetected. 

There is a companion galaxy G2 at ~ 4.1" from the 
main lens, towards south-east. Assuming the Faber-Jackson 
law and that G1 and G2 are at the same redshift, it would 
be expected that G2 produces a shear of ~ 0.05, however the 
measured shear assuming the SIS +7 model is just < 0.02, 
directed towards the foreground group of galaxies detected 


in SDSS and the HST data towards north-east (Section l5.2D . 
For this double, the most plausible explanation is that the 
there is an interplay between the effect of the intrinsic lens 
ellipticity, as well as the shear from the nearby companion 
and the nearby group, with the later being dominant. 

B5 SDSS J0904+1512 

The discov ery o f SDSS J0904-I-1512 was reported in 
iKavo et al.l d2010l) . It is a small-separation (~ 1.1") double 
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Figure B5. Hostlens modelling for the AO observation of SDSS J0832+0404. North is up and East is to the left. Top left-, original frame 
cut (10" X 10"). Top centre-, residuals after modelling the system, but without subtracting the lensing galaxy. Bottom left: residuals 
after modelling with an analytical PSF without fitting the host galaxy. Bottom centre: residuals after subtracting the quasar images and 
host galaxy, using an analytical PSF. Bottom right: residuals after subtracting the lens galaxy and quasar point sources, but without 
subtracting the host galaxy. 


at z s = 1.826. The redshift of the lensing galaxy has not been 
spectroscop i cally determined, and estimates are discrepant. 
iKavo et al.l (120 101) find zi ~ 0.2 based on the Faber-Jackson 
method ( Rusin et al.ll2003l ). but zi ~ 0.5 based on the V — I 
colour. 


Figure [B7] shows asymmetric residuals in the form of an 
arc (over-subtraction on one side, under-subtraction on an¬ 
other) at the location of A, oriented tangentially to the crit¬ 
ical curves. These are accounted for if a quasar host galaxy 
is considered. Since Figure IB7I shows non-analytical PSF 
residuals at the location of image A, these may affect the 
properties of the host galaxy fitted with Hostlens. Therefore 
an iterative approach was taken to construct a hybrid PSF 
which accounts for the non-analytical components, as de¬ 
scribed in Section ESI The hybrid PSF core consists of the 
observed pixel values (after subtracting the host contribu¬ 
tion) in a circle of 10 pixel radius, centred on image A. This 
was chosen to be large enough to include the non-analytical 
PSF components, but as small as possible to leave outside 
a visible contribution from the host galaxy. The iteration 
was stopped once the \' 2 stopped decreasing, in the third 
iteration. A similar iteration was also performed by fixing 
the Sersic index of the host galaxy at n = 4. The resulting 
host parameters are written in Table [6] Although it is very 


compact (~ 1 pixel effective radius in the source plane, but 
larger in the observed, magnified image plane), it is consid¬ 
ered a true detection because of the visually confirmed arcs. 
In order to perform simulations, the image A cut produced 
in the last step of the iteration was used, with the host pre¬ 
viously subtracted. This introduces additional noise in the 
simulations, however the image is quite bright. The same 
iterative procedure was applied to the simulations. The el- 
lipticity of the lens galaxy converges towards zero and was 
therefore ignored. The nearby companion G2 (Figure IB7I 
top left) was modelled with Galfit using the hybrid PSF. 

This object shows the largest differences compared to 
the discovery paper. It is also one of the smallest separation 
systems in the AO sample, which might explain the discrep¬ 
ancies. I 11 the discovery paper, A and B are only ~ 4 and 
~ 2.5 pixels apart, on each axis, with FWHM ~ 3.5 pix¬ 
els. The lens galaxy astrometry differs by ~ 1.5 pixels on 
each axis from the AO data. In the AO data, the lens G1 
appears clearly to be located closer to image B, as opposed 
to the discovery paper. While the effective radius and Sersic 
index are compatible, the discovery paper fits an ellipticity 
of ~ 0.5, while in the AO data it is consistent with 0. 

Lens modelling was performed by assuming zi = 0.3. 
Converting the SIE model result in the discovery paper with 
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Figure B6. Galfit/Hostlens modelling of the IIST observation for SDSS J0832+0404. North is up and East is to the left. Top left: 
original frame cut (75" X 75"), showing the location of the two PSF stars, and a group of nearby galaxies (circle). Top centre: original 
frame cut (10" X 10"). Top right: residuals after modelling the system, but without subtracting the lensing galaxy. Bottom left & bottom 
centre: residuals after modelling the system with PSF1 & PSF2, respectively. Bottom right: residuals after subtracting G2 and GX as 

well. 


the e ~ 37 relation gives a shear of ~ 0.24, the largest in 
this sample for a SIS +7 model. The value using the AO 
astrometry is ~ 0.08, more in line with the other objects; 
the orientation of the shear is different from that in the 
discovery paper at ~ 60 deg. This further supports the fact 
that the astrometry in the discovery paper is inaccurate. As 
the ellipticity was fixed at 0 , mass models incorporating it 
are not used. The nearby galaxy G2, assuming it lies at the 
redshift of Gl, would produce a shear of ~ 0.03, or less than 
half of the value required by the SIS +7 model. In addition, 
the shear direction is orthogonal to the direction towards 
G2. As a result, no conclusion can be drawn on the source 
of the shear, although the SDSS data shows that the system 
lies in a crowded environment. 


B6 SDSS J0946+1835 

The discovery of SPSS J0946+1835 is reported in 
iMcGreer et all J20ld) . At z s = 4.799, it is the highest red- 
shift lensed quasar known to date. The lens redshift has been 
determined spectroscopically to be zi = 0.388. The lens is 
a massive elliptical galaxy, causing a relatively large image 
separation ~ 3.1". Two faint components were discovered in 


the close proximity of the system, the bright one inconsistent 
with the colours of the images. 

The object was observed twice in the K' — band with 
AO, with the observation in 2012 being much deeper. This 
latter observation was performed in the 20 mas detector 
mode, with 9-point dither of separation 5". In order to in¬ 
crease S/N even more in the hope of the quasar host galaxy 
detection, 2x2 binning was performed. Due to the long ex¬ 
posure time, a significant number of cosmic rays were visible 
in each frame. These were removed automatically using the 
rejection-combine, however those in the immediate proxim¬ 
ity of the target were manually removed from each individual 
frame in order to prevent any contamination. 

Figure [2] shows that the separately observed PSF star 
cannot be used as a suitable PSF. Therefore, a hybrid PSF 
was used to model the system. The core component contains 
the observed light distribution within a circle of radius 12 
pixels centered on image A. A larger hybrid core size was not 
considered, because it might be affected by inaccurate sky to 
the North of image A, which was located close to the edge of 
the field in several dither sequences. Simulations were done 
using the separate PSF star. The system was simulated with 
this PSF, but all objects were simulated 0.3 mag fainter, in 
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Figure B7. Hostlens modelling for SDSS J0904+1512. North is up and East is to the left. Top left: original frame cut (13" X 13"), 
showing the location of G2 to the South. Top center left: original frame cut (5.1" X 5.1"). Top center right: residual after modelling 
with an analytical PSF, without fitting a host galaxy. Top right: residual after modelling with an analytical PSF and fitting the host 
galaxy. Bottom left: residuals after subtracting the quasar images and host galaxy, using a hybrid PSF. Bottom center: residuals after 
subtracting the lens galaxy and quasar point sources. Bottom right: residuals after subtracting all components. 


order to produce about the same counts as the real system, 
since the separate PSF has slightly higher Strehl ratio. 

The faint circular residuals seen in Figure IB 8 I around 
image B are due to masking when creating the sky frame, 
and not due to quasar host galaxy arcs. There are no visual 
signs of arc detection. A quasar host was fit using an ana¬ 
lytical PSF to model the system, which resulted in a bright 
compact host ~l/3 the flux of the point source, and the lack 
of visible arcs. It is likely that part of the quasar light is in¬ 
corporated in the host, whose magnitude it therefore overes¬ 
timated. This was tested by creating 100 simulated systems, 
using the separate PSF star, without a host, and fitted with 
an analytical PSF, both with and without a host. The differ¬ 
ence between the reduced °f the models with a host and 
without is —0.0197 ± 0.0014 (lcr of the scatter from the 100 
values), whereas in the two corresponding models of the real 
system it is x7 d -°- f - = 63425.86/74147-64450.97/74152 = 
—0.0138. It is therefore very likely that the host detection 
is spurious. Since this is a very high-redshift quasar, the 
study of its host galaxy would be interesting. The fact that 
there is no clear detection is due to the cosmological surface 
brightness dimming and the low intrinsic magnification of 
this system. 

Figure lB 8 l shows that the two faint components noted 
in the discovery paper are also detected in the AO data 
(G2 and G3). Both have elongated morphology, indicative 
of galaxies. Galaxy G2, which the discovery paper suggests 
to be similar in colour with Gl, is 5 mag fainter than G1 
in the new observations, and photometric redshift suggests 
that it is a faint local galaxy. For G3, there is no colour 
information to estimate a redshift. 

The shear direction is consistent in the SIS +7 and 
SIE +7 models, but different at ~ 35 deg from the orienta¬ 


tion of the galaxy. The effect of G2 and G3 on the lens model 
is small, modifying the direction of the estimated shear by 
~ 3 deg and ~ 1 deg, respectively. SDSS imaging (Figure 
IB3I) shows that there is a group of galaxies in the proximity 
of the system, towards NW, that are consistent with the lens 
galaxy redshift. The direction of the shear is consistent with 
this direction. In addition, the estimated shear values are 
quite large, consistent with the existence of nearby struc¬ 
ture. In conclusion, the lensing models are well explained by 
the observed environment, which dominates over the contri¬ 
bution of the intrinsic ellipticity of the lensing galaxy. 

B7 SDSS J1001+5027 

SDSS ,11001+5027 was reported in lOguri et al.l (120051) . 
based on observations in the V. R and I bands, as a lensed 
quasar with two images of intermediate separation ~ 2.9", 
and with a second galaxy in very close proximity to the main 
lens. The redshift of the main lens is determined spectro¬ 
scopically in llnada et al.1 (12012! ) to be 0.415, and the source 
redshift is z s = 1.839. The _time_delay betw een the two 
images is measured by iRathna Kumar et al. (12013! ) to be 
119.3 + 3.3 days. 

For this object, image A is isolated enough in the AO 
data to be used as a PSF. However, this leaves more residu¬ 
als than usual at the location of B and Gl, indicating that 
image A is not a very good approximation of the PSF at 
the location these components. Therefore the decision was 
made to model the system with the analytical PSF (Figure 
IB91) . 

For the simulations, the separate PSF star observed 
with this system was a natural choice, however it has a 
smaller FWHM than the real system, producing objects that 
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Figure B8. Galfit/Hostlens modelling for SDSS J0946+1835. North is up and East is to the left. Top left: original frame (12" X 12"), 
top centre-, residuals after modelling with a free Sersic index and a hybrid PSF. Top right: residuals after modelling with Sersic index 4. 


are too sharp. Therefore the system was simulated using an 
actual cut around image A as PSF. This likely overestimates 
the error bars, a s noise is added to an already noisy PSF. 

lOguri et al.l (120051 ) noted that the two lensing galaxies 
of SDSS J1001+5027 have colours consistent with those of 
early-type galaxies at 0.2 ^ z ^ 0.5. However, in the AO K' 
band G2 appears significantly fainter compared to Gl, sig¬ 
nifying that the two galaxies (both ellipticals, as suggested 
by the value of the best-fit Sersic index) may not be located 
at the same redshift. 

As time delays are measured for this system, the mass 
models can have an extra parameter compared to the others 
used in the present paper (for d.o.f. = 0). Due to its prox¬ 
imity to the system, G2 is expected to have an important 
effect on the observed configuration. Therefore it was taken 
into consideration as an SIS when modelling the system. For 
Gl, it was proceeded as follows: 

• The fiducial SIS+ 7 , SIE, SIE +7 models were used for 
Gl, but the Hubble constant h was allowed to vary as a free 
parameter. 

• As a b ove, bu t h was f ixed a t the value measured by 
IPlanck Collaboration et alJ J2013I ). and instead the relative 
masses of Gl and G2 were no longer constrained by the 
Faber-Jackson relation. 

• h and the relative strengths of Gl and G2 were fixed, 
but models were considered with one of the following two 
changes: instead of the singular isothermal profile, a power 
law profile of free slope ( 7 ' = 2 in the case of the singu¬ 
lar isothermal profile) was considered for Gl; in addition 
to shear, external convergence was also considered at the 
location of Gl. 

The critical curves and caustics for two of the models 
above are depicted in Figure [B9l (bottom right). Of course, 
the real model would be a combination of the models consid¬ 
ered above, but such a model would be significantly under¬ 
constrained. Therefore these models were considered sepa¬ 
rately. 

The one-dimensional probability distributions for the 
Gl density slope 7 ' and for h, in Figure IB 101 show that 
models which account for the effect of G2 require a rela¬ 
tively larger h and less steep slope. Also, from the two di¬ 


mensional probability contours, among others, the following 
conclusions can be drawn, where all other parameters have 
been marginalised over: 

• models with steeper mass profiles require less mass and 
more external shear; 

• models with more external convergence require less 
mass and less external shear; 

• models with more massive G2 require less massive Gl 
and less external shear; 

The most interesting result is that all spherical mass 
models considered have similar shear orientation, regardless 
of whether G2 is modelled or not. The shear direction in 
the SIS +7 model differs from that of the observed light pro¬ 
file of Gl by only ~ 20 deg. However the light profile of 
this galaxy is not elliptical enough for the quadrupole to 
account for the large measured shear. This suggests strong 
environmental influence, dominant over both the internal 
shear generated by Gl and the influence of G2. lOguri et al.l 
(2005) d oes report the exi stence of a group of galaxies (Fig¬ 
ure 8 in lOguri et al.l 120051 ) consistent with z ~ 0.2 towards 
north-west, however this is roughly perpendicular to the di¬ 
rection of the shear. 

Another fact suggesting that the influence of G2 is not 
strong is that, if it is considered, not just small convergence, 
but negative convergence is required to reproduce the sys¬ 
tem in the SIS/SIE (Gl) + 7 + k + SIS (G2) case. This is 
peculiar, since the SIS profile produces 7 = k, and conver¬ 
gence is a positive scalar quantity, whereas shear is a head¬ 
less vector (invariant under rotation by 180 deg); therefore it 
is expected that 7 < k, assuming the shear is caused by the 
added contributions from multiple SIS perturbers. On the 
other hand, negative convergence does exist for regions of 
space that are relatively void. This however is not the case, 
as is shown by the existence of the nearby group of galaxies 
in the foreground. Comparing the one-dimensional probabil¬ 
ity distributions for convergence, it is seen that 7 ~ k if the 
contribution of G2 is ignored. Also, a weaker G2 reduces the 
mass slope of Gl, making it consistent with an isothermal 
profile (Figure |B1 111 . 

Therefore, the models and the observed environment 
strongly suggest that G2 should have a small effect on the 
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system, corresponding to either a different redshift from Gl, 
or a mass much smaller than suggested from the relative 
brightnes s of Gl and G2, a nd the Faber-Jackson relation. 
Note that lOguri et al.l (l2005l l measure even closer brightness 
for the two galaxies at shorter wavelengths, and therefore 
without the AO data the discrepancy would be even larger. 
There are two explanations for this result: 

• The redshift of G2 is indeed different, as suggested from 
the different colours of the two galaxies in the visible and at 
I\' band. Therefore, its gravitational effect on the system is 
reduced. 

• Alternatively, the proximity of Gl and G2 could sug¬ 
gest that G2 is tidally stripped by the larger galaxy, which 
reduces its mass without ch angin g the overall luminosity of 
the galaxy ('e.g.. iTreu eTaLlboogi'l . 

This appears to be a complicated system, where the 
measured time delay cannot be straightforwardly converted 
into a measurement of h. A spectroscopic redshift of G2 is 
required for additional insight, as well as field spectroscopy 
to understand the origin of the large required shear. 


B8 SDSS J1055+4628 


This object was reported in iKavo et al.l ( 201dll as a small- 
separation (~ 1.1”) double at z s = 1.249. Ilnada et a ll (1201 2|) 
spectroscopic lens redshift of zi = 0.388. IKavo et all 


2010 h also identifies two nearby galaxy companions and 
notes that the one to the North has similar colours to the 
main lens, but the one to the South is much bluer. 

In the AO data, an entirely analytical PSF was em¬ 
ployed due to the low S/N, with only 8 frames of 60 s expo¬ 
sure obtained for this objects. Simulations were performed 
with the PSF star observed for SDSS J1405+0959. There 
are features around image A that may be interpreted as arcs 
due to the lensed host galaxy (Figure [B121 bottom right), 
however due to the low S/N and lack of improvement in 
the residuals, a host galaxy fit is not reported. The possible 
presence of arcs, which are unaccounted for, makes this a 
less reliable object to explore with the present technique. 
For comparison, the morphological parameters of the lens- 
ing galaxy are all consistent within ~ la with the discovery 
paper, except for the tens ellipticity, which is significantly 
larger in the AO data. 

The companion galaxies from the discovery paper, al¬ 
though visible, appear very faint in these shallow observa¬ 
tions. The magnitude of G2 (Figure [B12] top left) was mea¬ 
sured with aperture photometry, since it is too distant from 
the system to have a reliable PSF. 

In the SIS +7 model, the required shear is small, and its 
orientation is different at ~ 34 deg from the light orienta¬ 
tion. The shear increases significantly for the SIE +7 model, 
with a large change in the position angle. The effect of con¬ 
sidering G2 is to change the position angle of the shear by 
~ 3 deg. This effect is expected to be even smaller, since the 
SDSS photometric redshift of G2, z ~ 0.3, is lower than the 
spectroscopic redshift of Gl. The effect of the two nearby 
galaxies identified in the discovery paper is not considered. 
In addition to being very faint, they are significantly offset 
from the direction of the shear in all three models consid¬ 
ered, therefore they are unlikely to be important in the mass 


models. No environmental clues for the origin of the shear 
are found in the larger-scale SDSS environment. 


B9 SDSS J1131 + 1915 

This object was reported in lKavo et al.l d2010lf . It is a small- 
separation (~ 1.5”) double at z B = 2.915. A spectroscopic 
tens redshift is n ot availa ble, but the use of the Faber- 
Jackson relations ijRusin et al provides an estimate 

of zi ~ 0.3. 

Using the AO data, a hybrid PSF was created from 
image A. Two different core radii, of 15 and 30 pixels, were 
used. As this is a bright system, simulations were produced 
with a cut around image A itself. For the nearby galaxy 
G2 (Figure rB13ll . astrometry and aperture photometry were 
measured with the IRAF PHOT task. 

The astrometry is highly incompatible with the one in 
the discovery paper. In particular, the tens galaxy is closer 
to the brighter image in the original data, and to the fainter 
image in the AO data The problem lies with the original low- 
resolution, large seeing observations and not the AO data, 
because components B and Gl are clearly resolved in the 
AO imaging. The direction of shear in the SIS +7 model is 
consistent with the observed light orientation to within 7 
deg, suggesting that the intrinsic tens ellipticity is the main 
source of shear. The tens model in the discovery paper is 
highly inconsistent, due to the vastly different astrometry. 
Scaling its relative velocity dispersion by the Faber-Jackson 
relation, nearby galaxy G2 is expected to change the direc¬ 
tion of the shear in the SIS +7 model by 18 deg. However 
this effect would be smaller if G2 is located at a different 
redshift. 


BIO SDSS J1206+4332 

The discov ery o f SDSS J1206+4332 is reported in 
lOeuri et al.l (120051 '). as an intermediate-separation (~ 2.9”) 
doub le at z s = 1.789. Time delay was measured in 
lEulaers et alj (|2013l l to be 111.3 ± 3 days. The system ap¬ 
pears to have three tensing galaxies where the main tens, 
Gl, has an associated absorber at redshift zi = 0.748. The 
second galaxy, G2, is identified as a high redshift galaxy at 
z ^ 0.7 and, the third galaxy G3, is blue and may be a 
chance superposition of a foreground galaxy. 

Since the PSF is severely oversampled in the non-AO 
observations, 2 x 2 pixel binning was performed. As there is 
no star in the FOV, the PSF was built analytically. When 
modelling the system with Hostlens, only A, B, Gl and G3 
(Figure IB14I) were considered, and the other objects were 
masked. Despite the poor seeing, a cut around image A, af¬ 
ter analytically removing the contribution of other objects, 
was successful. Using this cut as PSF however produces sig¬ 
nificantly more residuals at the location of image B, and 
subsequently it is discarded. Nonetheless, in the absence of 
other suitable non-AO objects to use for simulation, the cut 
around image A itself was used to simulate the system, and 
the resulting simulations were fitted in the same way as the 
real system. 

I 11 addition to the known galaxies, a new one (G4) 
was identified to the South of G2, which explains the pe¬ 
culiar morphology of G2 in the deconvolved images from 
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Figure B9. Galfit/Hostlens modelling and lens model for SDSS J1001+5027. North is up and East is to the left. Left: original frame 
(10” X 10”), centre: residuals after modelling with an analytical PSF. Right: Lens models: critical curves / caustics are shown in cyan 
and magenta for the 2SIS+7, green and black for a 2Sersic+7 model (which assumes mass fallows light), respectively. Squares mark the 
positions of the images, crosses mark the position of the lenses, and triangles mark the position of the source. A colour version of this 
figure is available in the online version of the paper. 










Figure BIO. ID and 2D selected lens model plots for SDSS J1001+5027, created with Glafic and MCMC. The 2D contours represent 
the la and 2 a regions. 


lEulaers et al.l d2013l ). In addition, another galaxy G5 was 
modelled, further away from the system. Several of the mor¬ 
phological parameters of the model galaxies were marked 
by Galht as unreliable, and therefore no morphological pa¬ 
rameters are reported. However, G1 appears to have a large 
Sersic index ~ 10, and G3 is also modelled with a large in¬ 


dex, whereas the other galaxies appear to require low values. 
Also, G1 seems to have a large ellipticity of ~ 0.7, with the 
major axis at ~ — 70 deg. Given the quality of the data, we 
do not consider mass models incorporating the observed el¬ 
lipticity, and we also omit error bars on the estimated mass 
models. 
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Figure Bll. Degeneracy for SDSS J1001+5027 in the case con¬ 
vergence k = 0 (black) and k = 0.15 (red). An elliptical power low 
profile +7+K (Gl) + SIS (G2) mass model was used, with the rel¬ 
ative contribution of G2 free to vary, but with fix n (d.o.f.= — 1). 
The contours represent the 1 a and 2tr regions. The horizontal line 
shows the relative ratio of the Einstein radii as inferred from the 
Faber-Jackson relation. A colour version of this figure is available 
in the online version of the paper. 


This system is similar to SDSS J1001+5027, in the sense 
that a time delay measurement is available, but it is even 
more complicated, due to the presence of more galaxies. If 
the morphology is to be trusted, Gl and G2 are early-type 
galaxies, and it is likely that they are physically associated, 
since they have similar colours. They are therefore e xpected 
to dominate the mass models. lOeuri et al.l (120051 ) obtain 
that G2 is indeed required in order to reproduce the observed 
configuration, but based on the new data Gl is enough to 
produce a perfect-fit SIS +7 or SIE model. The counterpart 
model of the one in the discovery paper, based on the new 
data (Gl as SIE with G2 as SIS), produces similar results. 

According to the models summarised in Table [5j 
whether or not G2 is considered (scaled using the Faber- 
Jackson relation), the shear position angles are in relative 
agreement with the tentatively measured light orientation 
(within ~ 20 deg), suggesting that the intrinsic mass ellip- 
ticity is likely to dominate the quadrupole. The strength of 
G3 is constrained by the number of observed images, and 
modelling it with SIS at the same redshift as Gl puts an 
upper limit of G3/G1 < 0.3 (compared to < 0.7 from the 
Faber-Jackson law) on its velocity dispersion in relation to 
Gl (in a model where G2 is also considered). This is in 
agreement with the different colours of Gl and G3, in the 
discovery paper. 

Analogous to the models considered for 
SDSS J1001+5027, we obtain that Gl should have a 
slope of ~ 1.8 assuming a power law mass model, or 
alternatively, for an SIS profile, convergence k ~ 0.15 or 
Hubble constant h ~ 0.8. Models with or without G2 are 
more similar to each other in terms of these results than 
in the case of SDSS J1001+5027, and in particular the 


required convergence is positive. The critical curves and 
caustics for two representative models are shown in Figure 
IB14I (bottom right). 

The observed flu x ratio in iv'Mjand is ~ 1, whereas 
it is ~ 1.5 both in lOeuri et al.l (120051 ) (at smaller wave- 
le ngth), and usin g the method which estimated a time delay 
in lEulaers et al.l (12013 ). This value has a small effect on the 
results in Tabled 

As a final remark, during the review process of this pa¬ 
per, we have become aware of the analysis presented for 
this system bv lAgnello et al.l d2015l) . Based on higher qual¬ 
ity adaptive optics imaging and additional spectroscopy, 
they obtained improved morphology results and were able 
to show that G3 is in fact a lensed image of the quasar host 
galaxy. 

Bll SDSS J1216+3529 

SDSS J1216+3529 was reported in iQguri et all d2008al ) as 
a small-separation (~ 1.5”) double at z s = 2.012. Based on 
the galaxy magnitude and colours, a redshift zi ~ 0.55 was 
estimated. A secondary galaxy was identified close to the 
system. 

This object has a very poor AO correction (PSF FWHM 
~ 0.46”), and only 9 frames of 60 s exposure were obtained. 
As a result, 2 x 2 pixel binning was used. The system was 
modelled with an analytical PSF, and simulated with the 
separate PSF star from SDSS J1001+5027, appropriately 
binned and convolved to match the observed FWHM of the 
system PSF. The magnitude and position of the more re¬ 
mote, brighter galaxy G2 (Figure fBlSI top left) was obtained 
with the IRAF PHOT task. For the closer, fainter compo¬ 
nent G3 measurements were made using the same analytical 
PSF. 

Here, similar to the discovery paper, G3, whose magni¬ 
tude is unreliable and photometric redshift is uncertain, was 
ignored in the lens models. If physically associated with Gl, 
its effect on the lens model would be drastic however, mod¬ 
elling only Gl, the position angle of the shear is in agreement 
within 7 deg with the light orientation of Gl. In addition, 
the SIS +7 shear direction is also close to the direction to G3, 
within 10 deg. It is therefore likely that there is a degener¬ 
ate contribution of the effect of G3 and intrinsic ellipticity 
of Gl for this system. As a note, G2, whose photometric 
redshift is consistent with Gl, would change the estimated 
shear direction by 6 deg. 

B12 SDSS J1254+2235 

This is a small-separation ~ 1.6” double (ilnada et al.ll2009l ) 
at high redshift (z s = 3.626), with a disk-like lens at zi ~ 0.3, 
based on its R — I colour. 

In order to eliminate the residuals in the AO data, and 
obtain a suitable fit, two concentric Sersic profiles were re¬ 
quired for the lensing galaxy. These are likely representative 
of the bulge and disc components of late-type galaxies. The 
orientations of these components are relatively misaligned 
by ~ 35 deg. modelling was performed without using a hy¬ 
brid PSF, because the proximity of the objects and the use 
of two Sersic profiles are likely to affect the construction of 
the hybrid PSF. The system was simulated using the sepa¬ 
rately observed PSF star. 
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Figure B12. Galfit/Hostlens modelling for SDSS J1055+4628. North is up and East is to the left. Left: original frame (66.5” X 66.5”). 
Centre left: original frame (10” X 10”). Centre right: residuals after modelling with a free Sersic index and an analytical PSF; galaxy is 
not subtracted. Right: residuals after subtracting all components. 



Figure B13. Galfit/Hostlens modelling for SDSS J1131+1915. North is up and East is to the left. Left: complete FOV (66” X 66”). 
Centre left: original frame (10” X 10”), centre right: residuals after modelling with an analytical PSF. Right: the lens galaxy is kept, 
while the quasar images are subtracted. 


A new, faint object well-modelled by a Sersic profile 
was found just East of image B (Figure IB 161 top right). 
Since colour information is unavailable, it is impossible to 
tell whether it is located in the foreground, background, or 
physically associated with the lensing galaxy. 

An SIE +7 model is not employed, since the position 
angle to be used with this model is not obvious. It is worth 
noticing that in the SIS +7 model, the position angle cor¬ 
responds within 5 deg to that of the brighter Sersic profile 
needed to fit the morphology. It is therefore likely that the 
dark matter halo of the galaxy is aligned to this brighter 
Sersic profile. 

It must be noted that taking into account the compo¬ 
nent G2 significantly changes the mass models (7 ~ 0.044 
at —12.4 deg for the SIS +7 with G2). This is because of its 
proximity to the system, despite its faintness. However, since 
the redshift of this component is unknown, its real effect on 
the system cannot be estimated. 

B13 SDSS J1313+5151 

The system is reported in lOfek et ~ahl (120071 ) as a small- 
separation ~ 1.2” double at z s = 1.875, with spectroscopic 
lens redshift zi = 0.194. 

For this bright object, the standard deviation of pixels 
in a region containing B and G when fitting the object with 
a hybrid PSF and an analytical PSF is smaller in the case of 
the hybrid PSF, regardless of the PSF cut size. This means 
that the contribution from the non-analytical PSF structure 


is more dominant than the noise in the PSF wing. The ob¬ 
ject was modelled with a hybrid PSF built on the bright 
image A. The scatter between two different PSF core radii 
was included in the error budget: a small radius of 25 pixels, 
and another including the whole image A cut. The simula¬ 
tions were done with a cut of image A, after all other con¬ 
tributions were removed. This object is particularly bright 
enough that even accounting for additional noise introduced 
in the simulations, the resulting error bars are small. 

As a caveat however, this object has more significant 
residuals compared to the other in the current AO sample, 
amounting to as much as 10% of the flux peak of G, at its lo¬ 
cation (Figure [B17ll . Adding additional Sersic profiles would 
not improve the residuals, which are concentrated over a 
small area around the peaks of G and B. The reason for the 
residuals is likely the fact that Poisson noise is large at the 
peaks. 

Compared to lOfek et ~ai1 (l2007h , the image brightness 
rank has changed. The new brightness ranking, with the 
brightest image more distant from the lens, was predicted 
in that paper, with microlensing and extinction brought for¬ 
ward as explanations for the anomalous flux ratio and its 
observed strong dependence on wavelength. Indeed, both of 
these effects should be weaker at longer wavelengths. 

The shear value in the SIS +7 model is very small, and 
the m odel parameters are similar to the ones in lOfek et al.l 
£2007]). The observed ellipticity is ~ 3 times larger than the 
equivalent ellipticity e ~ 37 , however the position angles are 
consistent to within ~ 8 deg. 
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Figure B14. Galfit/Hostlens modelling and lens model for SDSS J1206+4332. North is up and East is to the left. Left: original frame 
(17" X 17"); centre: residuals after modelling with an analytical PSF. Right: Lens models: critical curves are shown in blue, and caustics 
in red for the 1SIS+7 model, cyan and magenta for the 2SIS+7 model, respectively. Squares mark the positions of the images, a cross 
marks the position of the lens, and triangles mark the position of the source. A colour version of this figure is available in the online 
version of the paper. 



Figure B15. Galfit/Hostlens modelling for SDSS J1216+3529. North is up and East is to the left. Left: complete FOV (65" X 65"). 
Centre left.: original frame (10" X 10"), centre right: residuals after modelling with an analytical PSF; the lens galaxy is not subtracted. 
Right: residuals after modelling with an analytical PSF and subtracting all components. 


B14 SDSS J1322+1052 

SDSS J1322+1052 was reported in lOguri et al.l d2008al) . It 
is a small-separation ~ 2.0” double at z 3 = 1.716, with the 
lens redshift estimate zi ~ 0.55 based on the Faber-Jackson 
relation. 

I 11 the AO data, this object is similar to 
SDSS J0904+1512, in the sense that a quasar host 
galaxy was found and modelled iteratively with a hybrid 
PSF (with a 10 pixel core radius). It is the brightest and 
most extended host galaxy in the current sample (Figure 
IB18I) . with an almost complete Einstein ring. 

It became apparent during the modelling that the fit 
is problematic, with a very large Sersic index > 10 for the 
lensing galaxy. This is the largest value encountered in the 
current sample. Therefore, another model was constructed 
in the same way, where the Sersic index is fixed at n = 4. 
The residuals look remarkably similar in both models, and 
the final x 2 /d.o.f values (obtained after the forth iteration) 
are similar. It is likely that the Einstein ring from the host 
galaxy biases the modelling of the lensing galaxy, introduc¬ 
ing systematics. In view of the difficulty in estimating the 
correct morphology of the two galaxies, the host galaxy mag¬ 
nitude in Table O obtained when the Sersic index of the lens 
is n = 4, serves as an upper brightness estimate. 

There are four galaxies in the close proximity of the sys¬ 


tem (Figure fB18l bottom right). The direction of the shear 
in the SIS +7 model is fairly close, within 17 deg, to the 
light orientation. The measured lens ellipticity is small, and 
including it in the SIE +7 model has the effect of lower¬ 
ing the amount of required shear. Both models have similar 
shear orientations. Among the nearby galaxies, G3 the clos¬ 
est changes the estimated shear direction by 10 deg, making 
it less compatible with the orientation of Gl, if modelled as 
an SIS scaled to the main lens through the Faber-Jackson 
law; the effect of the other galaxies is negligible. 

B15 SDSS J1330+1810 

B15.1 Background and morphological modelling 

This object was reported and studied in lOguri et al.l d200sh . 
based on visible and near-infrared observations. It is the 
only four-image lensed quasar in the current AO sample. It 
is a small-separation (~ 1 . 8 ”) system, with source redshift 
z s = 1.3 93 and spe c trosco pically measured lens redshift zi = 
0.373. lOguri et alJ ( 2008 ) reports no flux ratio anomaly in 
their modelling. 

This system was observed with AO twice. The deeper 
AO data from 2013, obtained with 20 mas pixel scale, was 
binned (2x2 pixels), since the PSF was oversampled. The 
imaging data from both 2012 and 2013 was fitted with an an- 
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Figure B16. Galfit/Hostlens modelling for SDSS J !254+2235. North is up and East is to the left. Left: original frame (10” X 10”), 
centre: residuals after modelling with 2 Sersic profiles, using an analytical PSF. Right: residuals after modelling with one Sersic profile. 



Figure B17. Galfit/Hostlens modelling for SDSS J ! 3 1 3+5 151. North is up and East is to the left. Left: original frame (10” X 10”), 
centre: residuals after modelling with a free Sersic index and a hybrid PSF. Right: residuals after subtracting only the quasar images. 


alytical PSF, as four point-sources and a galaxy. The galaxy 
was modelled as two concentric Sersic profiles, resulting in 
a good overall subtraction (Figure (Bl^). In the deeper 2013 
data, two Sersic profiles were not found to be a significant 
improvement over a single profile. As seen in Figure IB19I 
the residuals look worse in the 2013 data, as the objects 
seem to be overestimated on one side, and underestimated 
on another. The shallower 2012 data shows less signs of a 
host galaxy, which becomes slightly more prominent, in the 
form of arcs around images A, B and C, in the 2013 data. 
Unfortunately, due to the peculiar PSF shape in the data 
and the faintness of the arcs, the host galaxy could not be 
fitted convincingly, and it was ignored in the subsequent 
morphological modelling, with the exceptions of the simula¬ 
tions described below. 

When simulating the 2012 data, the separate PSF star 
in SDSS J0946+1835 was used, as it is the most similar 
in terms of the analytical approximation. The simulations 
were only used to infer error bars for the astrometry and 
photometry. The morphology of the lens in Table [4] and 
photometry of all objects in Tabled is taken from the deeper 
observations in 2013. When simulating the 2013 data, the 
separate PSF star from SDSS J0946+1835 was used, with 
a Gaussian convolution performed so that it matches the 


observed ellipticity and FWHM of the PSF estimated in 
the observed data. In order to check if the faint host galaxy 
introduces any biases, a model was created where the system 
was fitted with Hostlens as a lensed system with a lensing 
galaxy (SIE+ 7 ) model and four lensed quasar images. Since 
this does not produce the observed flux ratios, a dark lens 
was used near image D in order to produce the observed ratio 
(see the section below on lens modelling). A host galaxy was 
then fitted, and its brightness was varied until it matched 
the brightness of the observed arcs near the three brighter 
images. This resulting brightness map of the host was added 
on top of each of the 100 simulations of the system, which 
were refit as five unrelated objects, without a lensing mass 
model. 

The morphology of the lensing galaxy in the 2012 data, 
if a single Sersic component is fitted, is in good overall agree¬ 
ment with the morphology obtained in the 2013 data. The 
astrometry is quite consistent within error bars for A, B 
and C in the two data sets, in spite of the large residuals. 
Larger errors exist for the lensing galaxy, likely due to the 
less accurate PSF fitting of this more extended object. As a 
result, the average of the astrometry in the two observations 
was used, and the scatter between the two was used as the 
error bar on astrometry, calculated as the sample standard 
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Figure B18. Hostlens modelling for SDSS J1322+1052. North is up and East is to the left. Top left: original frame (8" X 8"), top 
centre: residuals after modelling with an analytical PSF without considering the host galaxy. Top right: residuals after modelling with an 
analytical PSF, and accounting for the host galaxy. Bottom left: residuals after modelling the system with a hybrid PSF and subtracting 
all components except the host galaxy. Bottom centre: residuals after modelling the system with a hybrid PSF and subtracting all 
components. Residuals when fixing the Sersic index of the lens at 4 or leaving it as a free parameter are visually indistinguishable. 
Bottom right: 20" X 20" region around the system. 



Figure B19. SDSS J1330+1810 imaging observations in 52 mas mode (2012; left, centerleft) and 20 mas mode (2x2 binning, 2013; 
centerright, right). Shown are original 5” X 5” frames and residuals after Hostlens/Galfit modelling. The lens is modelled as two concentric 
Sersic profiles in the 52 mas mode, and one Sersic profile in the 20 mas mode. 


deviation. This method accounts for any possible systemat- 
ics, as the data is of course independent between the two 
observations, conducted with different pixel size, with a dis¬ 
tinctively different PSF, and different profiles for the galaxy 
(two Sersic and one Sersic component, respectively). 

The photometry of the lensing galaxy is c onsistent be¬ 
tween the AO data and the AVband value in lOguri et al.1 
(l2008ll . Images A, B and C appear to have become brighter 
by about 0.15 — 0.3 mag, while D is ~ 0.25 mag fainter. 
Assuming that the error bars in lOguri et al.l (l2008h . which 
do not account for systematic errors coming from uncertain¬ 
ties in the PSF and the galaxy profile, are not significantly 
underestimated, this effect must be due to quasar intrinsic 


variab ility and microlensing. As remarked by lOguri et al.l 
(l2008ll . microlensing likely plays the main role: while intrin¬ 
sic variability would not change the observed flux ratios, 
particularly in view of the short estimated time delays (Ta¬ 
ble 0, the flux ratios A/B and A/C change (Figure [B20I) . 
particularly in the shortest wavelength band, whereas A/C 
does not. In addition, B and D are located closer to the 
galactic disk, and therefore more likely to encounter stars. 

The relative astrometry is consistent with the larger 
error bars in the discovery paper. The main difference in 
morphology is the lower Sersic index (n = 2.2 compared 
to 3.6). This is closer to t he fiducial value of n = 1 for a 
disk-like galaxy. lOguri et al.1 (120081 ) notice that the galaxy 
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is brighter than expected from the Faber-Jackson relation 
given its velocity dispersion (or alternatively should have a 
lower redshift), and surmise that this could be due to tidal 
stripping. However, they find no obvious companion galaxy 
near G, and no obvious signature of interaction. No such 
sign is observed in the AO data either. However, the new, 
lower Sersic index offers a more direct alternative. As shown 
in Figure IB22I an Sbc spectral template produces a photo¬ 
metric redshift that is in agreement with the spectroscopic 
redshift (x 2 /d.o.f.=1.7 for 3 d.o.f.). 

B15.2 Lens modelling 

Instead of the 10% error bars on the flux ratios, employed 
throughout this paper, since information is available in mul¬ 
tiple bands, the average flux ratios from Figure IB20I (ex¬ 
cept for those in E-band) and their scatter were considered: 
B/A=0.75+0.05, C/A=0.46+0.02, D/A=0.19+0.05. Incase 
of image D, this error is larger than 10%. 

The SDSS en vironment is presented in Figure 7 from 
lOguri et al.1 J2008I') . They identify two galaxy groups, one 
towards south-west, and a second to north-east, closers to 
the system ; in w hich G3 (using the original notation in 
lOguri et al.l 120081 1 is the brightest galaxy. G3 is located at 
a = 0.311, however for the majority of the other galaxies the 
photometric redshift is not accurate enough to distinguish 
which galaxies are associated with the redshift of G, or of 
G3. 

Here we consider two classes of mass models, one which 
accounts for the nearby galaxies (only those visible in the 
iC-band FOV), by modelling them as SIS lenses scaled rel¬ 
ative to the main lens through the Faber-Jackson relation, 
and one which does not account for the environment. In both 
of these cases, models that account, and do not account for 
the measured flux ratios are considered separately. A sum¬ 
mary of the models used, their corresponding \ 2 and d.o.f. 
are shown in Table nm Taking account of the environment 
is not a rigorous procedure, and is meant only for reference. 
This is mainly because it uses the untested assumption that 
the mass distribution is discrete, following individual galax¬ 
ies rather than tracing the potential of a galaxy group. In 
addition, galaxies outside the A _, -band FOV are not consid¬ 
ered, and uncertainties in the scaling based on the Faber- 
Jackson relation are also ignored. We also assume that all 
galaxies are at the redshift of G . How ever this effect is not 
large: according to lOguri et alj d2008h . it is estimated that 
a weak external shear 7 at z = 0.31 is equivalent to an 
external shear of ~ O .87 at z = 0.373. 

As this system is a quad, there are enough d.o.f. to 
model both the external ellipticity and shear. The first 
model considered here, SIE+ 7 , was also considered in the 
disc over y pa per. Accou nting or not for the measured flux ra¬ 
tios, [OgurTetlil] (hoosh finds x 2 /d.o.f.=1.14/4 and 0.33/1, 
respectively. They therefore conclude that the object ex¬ 
hibits no sign of flux anomaly, which is typically encountered 
in the case of quads. Based on the order-of-magnitude more 
accurate astrometry, the updated values of 46/4 and 17/1 
are found, respectively. The fits become relatively better it 
the environment is considered (20/4 and 14/1, respectively). 
These model are therefore significantly worse fits than sug¬ 
gested by the low resolution data, in terms of astrometry 
and flux ratios. 



Figure B20. Flux ratios of the SDSS J1330+1810 quasar im¬ 
ages from follow-up observations. Plotted are magnitude differ¬ 
ences between images B-D and the brightest image A, A 7774 .x = 
m. 4 — nix (V = B, C, D). Dotted horizontal lines indicate average 
values of individual ratios (excluding VjjHgg), which are adopted 
in the mass modelling. 


A possibility is that the center of mass does not coincide 
with that of the light. However, the minimum offset from the 
observed lens position, necessary to obtain a good fit (when 
the flux ratios are discarded) is ~ 30 mas, or ~ 150 pc at z = 
0.37. This is five times the offset found by l|Yoo et al.ll2006l f 
based on the analysis of four lensed systems with Einstein 
rings, and is therefore unlikely. 

Next, more complicated but physically motivated mass 
models were fitted to explain the observed astrometry and 
photometry. These models are SIE + 7 + 5 (where 5 refers 
to the third order perturbation component), Power law el¬ 
lipsoid + 7 , NFW + Sersic + 7 , SIE + 7 + substructure 
and Power law ellipsoid + 7 + 5. Only two of the models, 
SIE + 7 +substructure and Power law ellipsoid + 7 + 5, 
produce a formally good fit. The introduction of the 5 com¬ 
ponent (Tables l5l and IBU Figure [B21I) is suggested in terms 
of the lens being located in a crowded environment, which 
contributes more than a shear term. But the compatibility 
(best-fit offset ~ 16 deg and ~ 14 deg w/o and w/ mod¬ 
elling the nearby galaxies, respectively) between the orien¬ 
tation of this component and the orientation of the lensing 
galaxy disk makes it likely that the component is physi¬ 
cally associated with mass in the disk. The fitted mass slope 
~ 2.4 is larger than the average value of the SLACS lenses 
(IKoopmans et al.l 120091 ). but consistent with their intrinsic 
scatter. As a caveat however, the SLACS lenses are elliptical 
galaxies. 

For the SIE + 7 + substructure model, the position that 
a single perturber (satellite galaxy) would occupy in order to 
explain the observables was studied. There are two possible 
locations, in the vicinity of image A and D, respectively 
(Figure IB23I) . The question is whether such substructure 
would be detected in the deep AO image. For comparison, 
if the satellite is compact, so as to be fit by a point-source, 
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Figure B21. SIE +7 (left) and power law ellipsoid +7 + <5 (right) mass models for SDSS J1330+1810. Critical lines (blue), caustics 
(red) and time delay contours (gray) are shown. The cross shows the position of the lensing galaxy, and the squares show the positions 
of the images. A colour version of this figure is available in the online version of the paper. 


Table Bl. y 2 table for mass models of SDSS J1330+1810 


Model 

w/ flux 
w/o field 

w/o flux 
w/o field 

w/ flux 
w/ field 

w/o flux 
w/ field 

Sersic +7 + 5 

20/4 

15/4 

20/4 

15/4 

Sersic +7 

43/6 

29/3 


29/3 

NFW+Sersic +7 

11/2 


11/2 


Sersic+ 7 +substruct. 

43/3 

0/0 



SIE+ 7 +substruct. 

0 . 9/1 

o/-i 

1/1 

o/-i 

SIE +7 

46/4 

17/1 

20/4 

14/1 

SIE +7 (gal. pos. free) 

2/2 

o/-i 

0.6/2 


SIE +7 + 5. 

15/2 

0/0 

29/2 


Pow. law ellip. +7 

11/3 

0/0 

7/3 

0/0 

Pow. law ellip. +7 + 5 

1/1 


0.4/1 



The values are reported as x 2 /d.o.f. 


and has the same brightness as image D, it would have a 
velocity dispersion of ~ 95 km/s. The most likely velocity 
dispersion of 60 km/s, in case the perturber is close to image 
D, translates to a virial mass of 6.5 x 10 8 Mq, and ~ 21 
mag (~2 mag fainter than image D). As the PSF was built 
from the images themselves, and the satellite galaxy could 
be located very close to image D, it is difficult to infer a 
detection magnitude. On the other hand, there might be 
multiple smaller pertur bers causing the observed anomaly. 
iFadelv fe Keeton! (l2012l l have considered models with a few 
as well as a population of substructure for HE 0435+1223, 
a quad with flux anomaly, and found that both are favoured 
by the data. 

As a last remark, Figure lB24l plots the orientation of the 
external shear for the SIE + 7 and the Power law ellipsoid 
+ 7 + 5 model, in order to check which of G1 and G3 most 
affects the system. While the SIE + 7 clearly favours a shear 
direction towards G3, the Power law ellipsoid + 7 + 5 model 
does not distinguish b etween the t wo. 

To conclude, iQguri et al.l ^OOsI l , based on low- 
resolution data, notes that a standard SIE + 7 model can 
well reproduce both the astrometry and flux ratios in the 


system. Based on the new data, it was concluded that this 
is no longer the case. There is both astrometric and flux 
anomaly, but these can be solved with a smooth mass model, 
using a steep Power law ellipsoid + 7 + 5 model, where 5 
likely results from mass in the lens galaxy disk. A satellite 
galaxy perturber cannot be ruled out, but it is not required 
by the data in this new model. Taking into account (albeit 
incompletely) the nearby lens environment does not modify 
these conclusions, but it does attenuate the discrepancies in 
the SIE + 7 mod el, particula rly pertaining to the flux ratios. 
For comparison, ISluse et al. fl 201 i/ :l End 3 of the 14 quads 
with high-precision astrometry to have astrometric anomaly, 
when modelled with SIE + 7 and including the contribution 
of the nearest perturber. They do not consider more com¬ 
plicated models for the lensing galaxy. Spectroscopy of the 
field (which consists of at least two groups of galaxies at 
different redshifts) , such as performed for other lenses by 
IWong et~ah] 1201 li), is necessary to better account for the 
rich environment. 


B16 SDSS J1334+3315 

The discovery of this do uble-image s mall-s eparation lensed 
quasar was published in iR.usu et al.l ( 2011 1 . based on spec¬ 
troscopy, low-resolution imaging in the zIJHK s —bands and 
LGSAO observations in the JHKs— bands. This was the 
first object observed in the present AO imaging campaign. 
It is one of only two sub-arcsecond systems in SQLS. 

We noticed that the error bars on the fitted ph otome¬ 
try in the non-AO JHKs data also reported in IR.usu et al.l 
J 201 l|) are more uncertain than the published error bars, 
which included just statistical errors reported by Galfit. 
This broadens the photometric redshift estimate of the 
lensing galaxy to 0.5 — 1.7 ( 68 % confidence interval), for 
X 2 /d.o.f. ~ 1 , assuming the spectral template of an early- 
type galaxy. This is nonetheless lower than the source quasar 
redshift, in agreement with the gravitational lensing nature 
of the system. Although the Iz— bands are particularly un- 
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Figure B22. Photometric redshift for the SDSS J1330+1810 
lensing galaxy, obtained with Hyperz for an Sbc template. The 
vertical line marks the spectroscopic redshift of the lens. The ver¬ 
tical axis represents the probability of exceeding the x' 2 value 
obtained during the fit by chance, as a function of the numbers 
of d.o.f. (= 3). 

reliable given the galaxy faintness, including them in the fit 
lowers the redshift to 0.53 — 0.96 (68% confidence interval), 
best -fit value 0. 6 7 (x 2 /d.o.f. ~ 0.6). 

iRusu et al.l (1201 ll i give an analysis of the AO data with 
a Gaussian + Moffat PSF built from image A, and quote 
only statistical errors. Here, the analysis was redone using 
Hostlens and an analytical PSF fit simultaneously on all 
components (A, B and G). The updated astrometry, pho¬ 
tometry and morphology are summarised in Tables [3] and 

SI 

This is the only system in the current sample where the 
lensing galaxy is closer to the brighter image, which may 
make the PSF estimation more difficult. Also, the K'— band 
shows residuals around image B that may or may not be 
due to the quasar host galaxy (Figure IB25I) . Whether or 
not the (possible) host galaxy is incorporated in the models 
introduces significant changes in the estimated morphology 
and photometry (~ 0.6 mag) of the lensing galaxy. 

In addition, the photometry/morphology of the lensing 
galaxy is not considered trustworthy in the AO data, due 
to the following caveat. There is a significant number of 
frames, particularly in the HK'~ bands, that show a dark 
sky pattern. These were kept in the reduced data in order to 
increase the S/N, since overall a small number of frames were 
obtained, but may affect the photometric and morphological 
results. As a result, the SIE+7 model was not used in Table 
[5] In particular, there are vast differences in the photometry 
of the lensing galaxy com pared to the no n-AO data, a fact 
which was also noted bv lRusu et al.l (1201 il l. The conclusions 
obtained in that paper regarding this system do not change. 

B17 SDSS J1353+1138 

SDSS J1353+1138 was reported in llnada et al.l d2006i l as a 
small-separation ~ 1.4” double at z s = 1.629, with lens pho¬ 
tometric redshift estimate of zi ~ 0.25. llnada et al.l (|2006i f 
identify an additional component C in the vicinity of image 
A. They rule out several explanations for this component, 


concluding most likely that it is a foreground star, in agree¬ 
ment with its bluer colours. 

Since for this system image A is distinctly separated 
from B and G, the hybrid PSF method was employed. In 
the initial analytical modelling of all components, two case 
were considered, according to whether or not component C 
from the discovery paper is detected. In the case that compo¬ 
nent C is considered at the position inferred in llnada et al.l 
J2006h . the X' 2 value is lower than in the case that the com¬ 
ponent is not considered. However, comparing the results 
in both cases with subsequently generated hybrid PSFs re¬ 
veals that in fact there are noticeably more residuals close 
to image B (Figure IB26I bottom right), in the case that 
the hybrid PSF was built by first subtracting component C. 
Therefore, component C is considered undetected in the AO 
data, consistent with its known blue colours, in support of 
the foreground star hypothesis. 

For this object, as was the case with SDSS J1313+5151, 
analytical PSF modelling produces a larger pixel standard 
deviation in the area containing B and G than the hybrid 
PSF does, therefore showing that the actual PSF is highly 
non-analytical. The scatter between the results when using 
a hybrid PSF core of 10 pixel radius (i.e. the region that was 
cut around image A, while the rest was replaced with an an¬ 
alytical wing), and another much larger radius, was included 
in the error budget. In the noise simulations, the system was 
created using a cut of image A as PSF, and remodelled using 
the hybrid PSF. As in the case of SDSS J1254+2235, two 
concentric Sersic profiles were found to be necessary for a 
suitable modelling of the lensing galaxy. The galaxy has a 
large ellipticity, and the two Sersic components have orien¬ 
tations consistent with each other. 

The shear orientation in the SIS+7 model differs from 
the orientation of G1 by only 10 deg, suggesting that the 
lens ellipticity is the main source of the quadrupole. However 
in the discover y paper, the t wo or ientations differ by ~ 30 
deg, which led llnada et al.1 (|2006l f to the conclusion that 
they are inconsistent. For the SIE+7 model, the ellipticity 
of Gla was used, this being the galaxy bulge component, 
with the smallest ellipticity. The effect of a nearby galaxy 
G2, assuming it is located at the same redshift, is small, 
changing the orientation of the shear by just 4 deg. 

B18 SDSS J1400+3134 

This system is reported in llnada et al.l (l2009i f as a small- 
separation (~ 1.7”) double at z s = 3.317, with lens redshift 
estimate ~ 0.8. 

Although images A, B, and lensing galaxy G1 are suf¬ 
ficiently separated, the S/N is small, and as a result there 
are no significant residuals when the PSF is modelled an¬ 
alytically ('Figure [B27I) . If the Sersic index of G1 is a free 
parameter, the preferred value is very large, ~ 10. Therefore, 
for this object, the Sersic index has been fixed at the fidu¬ 
cial value of 4, typical for an elliptical galaxy. This produces 
virtually no visual increase in the residuals at the location 
of the galaxy. Simulations were performed with the star in 
the SDSS J0743+2457 FOV and its subsequent analytical 
fit. For the nearby galaxy G2, its position and magnitude 
were determined using the IRAF PHOT task. 

The orientation of the shear in the SIS+7 model is dif¬ 
ferent at only ~ 14 deg from the orientation of Gl. While 
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Figure B23. SIE +7 (left) model of SDSS J1330+1810 with the location of a possible single substructure superposed. Contours are 
drawn from MCMC at la and 2a. All other parameters are marginalised over. The two smaller plots represent velocity dispersion 
probabilities for the perturber. Vertical bars contain 68 % of the probability. A colour version of this figure is available in the online 
version of the paper. 



Figure B25. Hostlens modelling of the SDSS J1334+3315 AO data. North is up and East is to the left. Top left: residuals after 
subtracting all components in K' band (4" X 4"), without considering a host galaxy, top center: residuals after modelling with an 
analytical PSF and considering the host galaxy. Top right: residuals after modelling with an analytical PSF, but not subtracting the 
host galaxy. 


this suggests that the intrinsic ellipticity of G1 is the main 
source of the quadrupole, it is likely that the nearby galaxy 
G2 also affects the model. Indeed, by adding G2 as an SIS 
perturber to the SIS +7 model, the shear position angle is 
in perfect agreement with the orientation of Gl. 

B19 SDSS J1405+0959 

This system was discovered independently as a two - image 
lens ed quasar candidat e by IJackson, Ofek. fe Qguril (12009! ) 
and llnada et alj feoidh The AO observations, subsequent 
analysis and res ults are presented in a separate paper, 
iRusu et al.l d2014l ~) . There it is shown that the system is likely 
to be a triple, produced by two galaxy lenses. 


B20 SDSS J1406+6126 

This system is reported in llnada et al.l (120071') . as a double 
with z s = 2.13, zi = 0.27 and small image separation 1.08”. 
The morphology of the lensing galaxy is not presented in 
that paper. 

The AO correction in the present campaign was very 
poor, with PSF ~ 0.75”. The PSF at image A is well approx¬ 
imated by as single Moffat profile. The system was simulated 
using this analytical PSF. The lens requires 2 concentric 
Sersic profiles (n = 1 and 4, respectively; Figure [B28I) . The 
flux ratio is in agreement with the discovery paper (V and I 
bands). The orientation of shear in the SIS+shear model is 
in agreement with the observed light orientation. 
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Figure B26. Galfit/Hostlens modelling for SDSS J1353+1138. North is up and East is to the left. Top left : complete FOV (62" X 62"), 
top centre-, original frame (9" X 9"); the possible component C is marked with a cross. Top right-, residuals after modelling with 2 Sersic 
profiles and a hybrid PSF; the lens galaxy is not subtracted. Middle left : residuals after modelling with a 1 Sersic profile. Middle centre: 
residuals after modelling with 2 Sersic profiles and subtracting all components. Middle right: residuals after accounting for component 
C. A colour version of this figure is available in the online version of the paper. 



Figure B27. Galfit/Hostlens modelling for SDSS J1400+3134. North is up and East is to the left. Left: 62" X 62" FOV. centre left: 
original frame (10" X 10"), centre right: residuals after modelling with Sersic index 4 and an analytical PSF; lens galaxy is not subtracted. 
Right: residuals after subtracting all components. 


B21 SDSS J1455+1447 

The system was reported in iKavo et al.l (l2010l l as a dou¬ 
ble at z s = 1.424, with image separation 1.73”. The AO 
data is suitable to construct a hybrid PSF from image A. 
Simulations were performed with the separate PSF star 
from SDSS J0832+0404. The inferred relative astrometry 
is severely incompatible with the discovery paper, with all 
components being clearly resolved in the AO data. This also 
translates to large discrepancies in the inferred mass models. 
The orientation of the shear in the SIS +7 model is close to 
both the orientation of the luminous profile of the lensing 
galaxy (at ~ 7 deg) and the direction to the nearby com¬ 
ponent G2 (at ~ 17 deg; Figure [B29h . and an even better 


agreement (at ~ 1 deg) is obtained with the light orientation 
if G2 is explicitly accounted for in the mass model. 

B22 SDSS J1515+1511 

The d iscovery of this object is reported in llnada et al.l 
d2014h . It is a small-separation (~ 2 . 0 ”) double at z 3 = 
2.054, with spectroscopically determined lens redshift zi = 
0.742. 

This is the only object for which a suitable star to use 
as PSF exists in the FOV. The LGS was centred between 
the system and the nearby bright star, to insure as much 
as possible the same AO correction. To discard pixels that 
might be unreliable due to Poisson noise, the 2x2 pixels at 
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Figure B28. Galfit/Hostlens modelling for SDSS J1406+6126. North is up and East is to the left. Left: original frame cut (10” X 10”), 
showing the location of G2 to the South. Centre left: Residuals after subtracting images A and B, modelled with an analytical PSF. 
Centre right: residuals after subtracting all components and modelling the lens with a Sersic index of 4. Right: residuals after subtracting 
all components and modelling the lens with a free Sersic index. 



Figure B29. Galfit/Hostlens modelling SDSS J1455+1447. North is up and East is to the left. Top left: original frame cut (20” X 20”), 
showing the location of G2 to the West. Top centre: original frame cut (10” X 10”). Top right: residual after modelling with a hybrid 
PSF, without subtracting the lens galaxy. Bottom left: residual after subtracting the lens galaxy as well, modelled with a free Sersic 
index. Bottom centre: residual after subtracting the lens galaxy as well, modelled with Sersic index 4. 


the peaks of A and B were masked, resulting in a quasar host 
galaxy that is three times as extended but half as luminous 
compared to the case in which the objects are not masked. 
100 simulations with noise were also performed and mod¬ 
elled in the same way. The parameters of the host galaxy 
are shown in Table [S] 

A question that needs to be addressed is to what extent 
the star is a suitable PSF, especially considering the claim 
that the host galaxy is discovered. First, the analytical pa¬ 
rameters of the star and the fitted system PSF in Table [A~2l 
are nearly identical. Second, the extended residuals seen in 
case that the host galaxy is not considered cannot be caused 
by anisokineticism, since the tip-tilt star is located closer to 


the system, and not the star; therefore any effect due to 
anisokineticism would be expected to broaden the star, and 
not the system. Third, the host model produces a good fit 
which results in improved residuals, and the morphological 
parameters of the host are reasonable. 

A nearby bright galaxy G2 is located in the frame, and 
its position and magnitude was measured with the IRAF 
PHOT task. The same was done for a faint galaxy located 
closer to the system (20.18 ± 0.07 mag). 

The lens appears to be an edge-on disk-like galaxy, with 
very large ellipticity. The shear position angle in the SIS +7 
model is consistent with the galaxy orientation at about 
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Figure B24. SDSS J1330+1810 la and 2 a contours of the 
position angles of the ellipticity and external shear. All other 
model parameters are marginalised over. Flux ratios are included 
as c onstraints . T he len s models used are SIE +7 from the orig¬ 
inal [OgurTliFld] d200§l) observations (black), as well as SIE +7 
(red) and power law ellipsoid +7 + <5 (blue) respectively, using 
the IRCS+AOI 88 observations. Contours are estimated using 
Ax 2 = 2.3 and 6.2 for d.o.f. = 2 in case of the original data, 
and MCMC for the IRCS+AOI 88 observations. The crosses in¬ 
dicates the best-fit values. The vertical dotted line indicates the 
observed light orientation. Horizontal dotted lines are expected 
values when the external perturbations are caused by the centre 
of the group ( G1/G2) or the nea rby bright galaxy (G3; following 
the notation in lQguri et al.1120081) . A colour version of this figure 
is available in the online version of the paper. 

~ 10 deg, therefore the shear is likely caused by the galaxy 
intrinsic ellipticity. 

We checked that the expected influence of G2 and an 
even closer galaxy to the West are minimal, as including 
them scaled by the Faber-Jackson relation changes the shear 
orientation by only and ~2 deg, respectively. The SDSS 
photometric redshift of G2 is consistent with the redshift 
of Gl, while the nearby perturber appears to be a disk-like 
galaxy. 

B23 SDSS J1620+1203 

The discov ery of SDSS J1620+1203 is reported in 
iKavo et al.l |2Q10l l. It is an intermediate-separation (~ 2.8”) 
double at z B = 1.158, with spectroscopic lens redshift zi = 
0.398. 

As the 2011 data consisted of only two frames, and 
therefore is significantly noisy, in addition to containing dark 
detector pattern, it was discarded entirely, except to check 
that the A/B flux ratio is consistent with the newer obser¬ 
vations. 

Since the system was observed both times in non¬ 
photometric conditions, we proceeded as follows. The sys¬ 
tem is detected in 2MASS with K = 15.203 + 0.136; the 
image separation is 2.8”, but the 2MASS pixel size is 2”, 
and the seeing was 2.7”, so the system is not resolved in 
2MASS. The photometric the zero-point was determined by 
equating the integrated magnitude of the objects, measured 


with Hostlens and Galfit, with the 2MASS value. Of course, 
this does not take into account any changes in the inte¬ 
grated magnitude of the quasar images since the 2MASS 
observations on 24 March 1999. Such effects are unlikely 
to be large, since the galaxy is ~ 2 mag brighter than the 
brightest quasar image. In addition, the lens is an elliptical 
galaxy with known redshift. From the Faber-Jackson rela¬ 
tion, with the velocity dispersion i nferred from the im age 
separation, and calibrated with the iRusin et al.l (l2003h for¬ 
mula, the predicted magnitude of the lens is K = 15.05. 
The lens magnitude estimated with the previous method is 
K = 15.50 ± ~ 0.04. Adding to this the uncertainty in the 
2MASS observation as zero-point, and the scatter of ~ 0.5 
mag in the Rusin formula, these estimates are in agreement. 

The morphological modelling of this system was per¬ 
form ed by creating a hybrid PSF with the core radius of 15 
and 30 pixels. Simulations were made using the PSF star of 
SDSS J1405+0959. 

The observed light orientation is different from the 
shear direction in the SIS +7 model. There are several galax¬ 
ies around the system in the IRCS+A0188 FOV, and in ad¬ 
dition there are several galaxies in the SDSS in the vicinity 
of the system (Figure IB3I) . towards south and south-west, 
roughly the direction of the estimated shear. 


APPENDIX C: SUMMARY OF BINARY 
QUASARS AND SYSTEMS WITH 
UNRELIABLE DATA 

Cl SDSS J0926+3100 

SDSS J092634.560+305945.90 has been selected by SQLS as 
a gravitationally lensed quasar candidate from SDSS DR(0 
The system does not appear to have a published spectrum in 
SDSS, and no follow-up spectrum was obtained. AO imag¬ 
ing of this system is presented in Figure Eu The system 
appears to consists of two well-resolved components sepa¬ 
rated by ~ 2.00”. Faint extended emission remains around 
the faint component after subtraction of an analytical PSF 
built on the bright object. Whereas the extended emission 
can be accounted for with an additional Sersic profile, it 
cannot be interpreted as a lensing galaxy, as its preferred 
location is just ~ 1 pixel offset from the point source, in 
the direction opposite to the bright component. Therefore 
we interpret this system to consist of a galaxy (the faint 
component to south-west) and either a quasar or a star (the 
bright component). 

C2 SDSS J1002+4449 

This system has been reported in llnada et al.l (|2008| J as a 
two-image lensed quasar candidate. Several follow-up obser¬ 
vations have been inconclusive. Spectroscopy in the SDSS is 
available for only one of the two quasars, which is located at 
2 = 1.998. The separation between the two components is 
~ 0.7”, which would make it the smallest-separation system 
in SQLS. 

The observations were conducted with the FOV (and 
therefore the LGS) centered in between the system and the 

10 https://www.sdss3.org/dr9/ 
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Figure B30. Hostlens modelling for SDSS J1515+1511. North is up and East is to the left. Left: 66 " X 66 " FOV. Top centre: original 
frame (10" X 10"), top right: residuals after modelling with a free Sersic index and the nearby PSF star. Bottom centre: residuals after 
modelling with a host galaxy. Bottom right: residuals after modelling without a host galaxy. 



Figure B31. Galfit/hostlens modelling for SDSS J1620+1203. North is up and East is to the left. Left: SDSS J1620+1203 FOV 
(64" X 64"). Centre left: original frame (10" X 10"), centre right: residuals after modelling with a free Sersic index and an analytical 
PSF. Right: residuals after modelling with a Sersic index 4 and a hybrid PSF. 


bright star to the North. However after eliminating bad- 
quality frames the configuration (PSF star, LGS, system) is 
no longer axially symmetric with respect to the LGS. As a 
result, the bright star has proven to be a bad PSF match 
(Figure lG2l top right). Analytical PSF fitting with Hostlens 
produces much improved results (Figure fC2l bottom). 

For this unconfirmed lensed candidate, the goal was to 
check if there is a lensing galaxy in between the quasars. 
Therefore the system was modelled both as two point-like 
objects, and as two point-like objects with a Sersic profile 
in between. When modelling the morphology of the system 
with a lens galaxy, x 2 /d.o.f. = 102359/90579 = 1.1300 was 
obtained, and x 2 /d.o.f. = 102612/90585 = 1.1328 with just 
two quasars, for a difference of —0.0028. The system was sim¬ 
ulated 100 times as two point-like objects, using the bright 
star as a PSF, and then remodelled with analytical PSFs, 


both with and without a galaxy. This produced an average 
difference of —0.0027±0.005 (la of the distribution from the 
100 values) between the models with and without a lensing 
galaxy. This strongly suggests that the galaxy detection is 
spurious. 

In addition, if the system is modelled with an SIS +7 
model, on the assumption that the galaxy detection is real 
and its derived relative astrometry is accurate, a shear of 
~ 0.27 is derived. This is the largest value obtained in the 
present sample for a system with a single lensing galaxy. 
There is no environmental clues in the FOV to suggest an 
explanation for this large shear. 

We conclude that the detection of a lensing galaxy can¬ 
not be claimed, making this object most likely a binary 
quasar. In principle, if this is a lens system, an estimate of 
the upper limit of the galaxy luminosity may be compared 
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Figure Cl. Galfit/Hostlens modelling for the gravitationally lensed quasar candidate SDSS J0926+3100. North is up and East is to 
the left. Left-, original frame cut (10” X 10”). Centre: residuals after modelling with an analytical PSF, without fitting a lensing galaxy. 
Right: residual after modelling with an analytical PSF, with a lensing galaxy. 


to its expected luminosity, derived from the observed quasar 
separation through the Faber-Jackson law, for a reasonable 
redshift range. However, the low S/N, the unavailability of 
an accurate PSF and the small separation does not permit 
a reliable estimate, and the best-fit estimate for the galaxy 
magnitude, ~ 17.15 mag, is brighter than some of the other 
lenses in this paper, which means that the lens hypothesis 
cannot be entirely disregarded. 


C3 SDSS J1054+2733 

The system was reported in iKavo et al.1 (2010:) as a two- 
image lensed quasar at z = 1.452, with an image separation 
~ 1.27”. The new AO data is shown in Figure [(ATI There 
are prominent residuals left at the centre of images A after 
the subtraction of an analytical PSF, as well as a flux ex¬ 
cess in the vicinity of image B, in between A and B. Since 
a suitable independent PSF is unavailable, the hypothesis 
that the flux excess is due to non-analytical PSF residuals 
at image B cannot be discarded. If on the other hand this 
additional component is real (component C in Fi gurelC3 [l. it 
has no t been detected at shorter wavelengths in lKavo et alj 
(l2010j) . The new component cannot still remains if the lens¬ 
ing galaxy is modelled with two sersic profiles, and it requires 
a point source to model satisfactorily (Figure IC31 bottom 
right). On the other hand, the lensing galaxy appears to be 
offset towards south. 

The astrometry of this system, if accounting for com¬ 
ponent C, is given in Table EU Due to the alignment of A, 
B, C, and the offset of the lensing galaxy, we were not able 
to fit a reasonable lensing model that can account for com¬ 
ponent C either as an extra quasar image, a central third 
image (with a non-singular isothermal sphere), or as part 
of a four-image system. If C was a chance superposition of 
a star, it would be expected to be brighter at visible wave¬ 
lengths, although its non-detection in the discovery paper 
may be attributed to the poor resolution and small image 
separation. We conclude that the AO data of this system is 
not well characterised. 


Table Cl. Inferred Relative Astrometry and Absolute Pho¬ 
tometry for SDSS J1054+2733 


Object 

AX ["] 

Ay ["] 

K' [mag] 

A 

0.000 

0.000 

15.54 

B 

-1.226 

0.144 

17.01 

C 

-0.869 

0.086 

18.25 

G 

-0.717 

-0.389 

17.65 


C4 SDSS J1128+2402 

This system is reported in llnada et al.l (|2014) ) as a two- 
image lensed quasar at z = 1.608, with an image separation 
~ 0.84”. In the new AO data, the inferred image separa¬ 
tion is 0.75”, making this system the smallest-separation 
lensed quasar in the SQLS. However, the lensing galaxy is 
undetected in the new observations (Figure IC4I) . with the 
residuals being identical whether a lens galaxy is modelled 
or not. Without visual confirmation of improvement in the 
residuals, the extracted physical parameters of the lens are 
unreliable, as it is likely that noise is being fitted. Con¬ 
sidering the very small separation of this system and the 
fact that there are non-analytical residuals which cannot be 
w ell fit, we conclude that the faint lensing galaxy, visible 
in llnada et al.l (120141 '). cannot be reliably resolved from the 
bright quasar images with the current technique. 

C5 SDSS J1320+1644 

The discovery of SDSS J1320+1644 and subs equent data 
analysis based on non-AO data is reported in IRusu et al.l 
ll2013l) . This is a large-separation quasar pair, likely to be 
images of a lensed quasar source at z s = 1.502, by a lens¬ 
ing galaxy group/cluster at zi = 0.89 9. However the lensing 
nature of this object is not secure. IRusu et ~ahl (120131 ) noted 
that the most efficient way to test the lensing hypothesis is 
by performing AO observations sensitive enough to detect 
the quasar host galaxy, and search for evidence of lensed 
arcs. 

The acquired AO data was affected by cloud coverage, 
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Figure C2. Galfit/Hostlens modelling for SDSS J1002+4449. North is up and East is to the left. Left: (64" X 64") FOV. Top centre: 
original frame (10" X 10"), top right: residuals after modelling the two quasars with the PSF star. Bottom centre: residuals after modelling 
the system with a galaxy of Sersic index 4 and two analytical PSFs. Bottom right: residuals after modelling with an analytical PSF, 
without a lens galaxy. 
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Figure C3. Galfit/Hostlens modelling for SDSS J1054+2733. North is up and East is to the left. Left: original frame cut (10" X 10"). 
Centre right: residuals after modelling with an analytical PSF, without fitting the lens galaxy. Centre right: residual after including the 
lens galaxy in the fit. Right: residual after modelling A, B, C, G simultaneously. Crosses mark the relative positions of all four objects. 
A colour version of this figure is available in the online version of the paper. 


and as such did not reach the intended depth, being in fact 
shallower than the original non-AO data. Even after 2x2 
pixel binning to increase S/N, there are no visible signs of a 
host galaxy detection (Figure fC5l) . Both an analytical PSF 
and a PSF built on image A were used, leading to clean 
residuals for image B. 

The only new insight brought by the new AO data is 
on the quasar flux ratio. This is A/B ~ 1.8, and is signif- 
icantly larg e r tha n the values of 1.4 and 1.0, considered in 
iRusu et ali (|2013l ). While this has little impact on their in¬ 
ferred lens models, since the main model constraints were 
astrometric, it does impac t their discussion on the intrin¬ 
sic flux ratio. iRusu et alJ J2013I ) argued that, since typi¬ 
cal changes due to intrinsic variability and microlensing are 
small in the near-infrared, the flux ratio ~ 1.4 observed in 
the spectra at long wavelengths and in the JHK S — bands is 


representative of the true flux ratio. An alternative hypoth¬ 
esis was that the flux ratio ~ 1.0 found at longer infrared 
wavelengths is more accurate, and the changes due to mi¬ 
crolensing at shorter wavelengths occur on time scales longer 
than the ~ 1 year spanned by the observations available at 
that time. The new flux ratio supports this second hypoth¬ 
esis in the sense that larger flux variations do indeed occur 
on time scales of ~ 3 years, assuming, of course, that the 
gravitational lens hypothesis is correct for this system. 

This paper has been typeset from a T^jX/ lATjjjX file prepared 
by the author. 
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Figure C4. Galfit/Hostlens modelling for SDSS J1128+2402. North is up and East is to the left. Left : original frame cut (10" x 10"). 
Centre : residuals after modelling without a lens galaxy. Right: residuals after modelling with a lens galaxy. 



Figure C5. 16 // x 16" cut around SDSS J1320+1644. Left: original frame. Object names correspond to those in iRusu et alJ (|2013tl . 
Right: residuals after modelling A, B, G1 and G2 using a hybrid PSF built on image A. 
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